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FOREWORD 


This Final Report is submitted in 
accordance with the Statement of 
Work, Exhibit "A" for Contract 
NAS8-32660. The study was dir- 
ected from the Guidance Systems 
Division (GSD) of The Bendix 
Corporation. The program manager 
at this location for most of the 
contractural period was Mr. James 
Jennings. Contributors from GSD 
were Mr. Raymond Kaczynski (Sections 
1 through 4) and Dr. Frederick 
Chichester (Section 7). Other 
tasks were completed by personnel 
from the Bendix Research Laboratories 
Division (BRL) and the Bendix Energy, 
Environment and Technology Office 
(BEETO). Material in Sections 5 and 
6 was contributed by Dr. Kelvin Smith, 
Don Johnson and Dave Sidlosky. Mr. 

Art Cornell of BRL generated the 
dynamic mathematical models for the 
Space Construction Base and Mr. 

Donald Lipski of B.IL wrote the 
programs for digital simulation. Mr. 
Calvin Rybak of BEETO wrote a tech- 
nical paper which provides the basis 
for the material presented in Section 
8. The guidance of Dr. Michael Borelli 
of MSFC throughout the study is 
gratefully achnowledged. 



ABSTRACT 


Many aspects ot an attitude control system are studied and de- 
veloped for a large space base that is structurally flexible and 
whose mass properties change rather dramatically during its 
orbital lifetime. Topics of discussion Include the following: 

a. Space base orbital pointing and maneuvering 

b. Angular momentum sizing of actuators 

c. Momentum desaturation selection and sizing 

d. Multilevel control technique applied to 
configuration 1 

e. One-dimensional model simulation 

f. N-body discrete coordinate simulation 

g. Structural analysis math model formulation 

h. Discussion of control problems and control 
methods 
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SECTION 1 


1 . 0 INTRODUCTION 

This report is submitted in compliance with the 
Statement of Work contract NAS8-32660 Exhibit "A." 

The period of performance covered by the report is 
the calendar year of July 27, 1977 to July 27, 1978. 

The submission and approval of this report constitute 
the successful completion of the Exhibit "A" portion 
of the contract. Additional studies, basically a 
continuation of the present work effort, will be 
conducted under an Exhibit "B" Scope of Work. 

1 . 1 OBJECTIVES 

The sections that follow summarize the effort expended 
on the Space Construction Base (SCB) Control System 
Study contract. The topics discussed are diverse in 
content but all represent a portion of the overall 
study. The primary objective of the study was to 
develop a control system and flexible control tech- 
niques that will stabilize a large and growing space 
station of the future. 

1 . 2 SCOPE 

The thrust of the program occurred in several unique 
directions. A hypothetical space mission was developed 
U'"ing the basic information generated by NASA, and 
requirements and timeline were formulated. From this 
initial baseline the effort branched into several 
categories : 
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1. Defining a mathematical model to be used in an 
all up and limited stability model. 

2. Investigation of modern control techniques as an 
Integral part of the control system. 

3. Development of a mathematical model that would 
describe the flexible structure characteristics of 
a large vehicle assuming the model characteristics 
of the individual pieces are known. 

4. Determine sizing requirements for a momentum 
storage system taking into account gravity gradient, 
magnetic, aerodynamic and radiation torques. 
Requirements were determined for: 

a. All configurations, various orientations 
(X-POP, X-LV, XP-POP, etc.) 

b. Inertial pointing modes anywhere in the 
celestial sphere (configurations 4 through 7). 

c. Earth pointing modes (configuration 12). 

The results of these studies have become the bas s for 
the development of the control system. At preseiic, 
the simulation effort is in progress and is expected 
to continue into the extension effort. The structural 
math model effort will also be continued. An attempt 
will be made to validate the model by using known 
spacecraft data. Other tasks such as the mission 
profile, multilevel analysis and angular momentum 
sizing ari! essentially complete in themselves, and will 
be used only as reference material during the conti- 
nuation . 
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1 . 3 GENERAL 

Several appendices containing detailed equations and 
modeling are Included at the end of the report. When 
necessary, sequences of equations are included in the 
body of the report. 

The RFQ requested that the International System of 
units (designated as SI) be used in the program and 
in any reporting. Expression in customary units would 
be acceptable if it is useful to the primary recipients, 
but SI should be stated first with the customary units 
afterwards, in parenthesis. In the report, SI units 
are often used (magnetic and radiation); torques, 
angular momentum, moments of inertia and distances, 
however, are stated in English units since this was 
the method used in presenting all o^ the data in the 
RFQ. 


A liberal approach was also taken in the numbering of 
SCB configurations. The RFQ initially used Roman 
numerals, but the text here uses Arabic numerals in 
most cases — except where material is reprinted from 
the RFQ. 
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SECTION 2 

2 0 MISSION DEFINITION 

A mission profile is defined which identifies those 
parameters which are the necessary inputs to the con- 
trol system definition. Those parameters are defined 
in reference 2-1 using information obtained from 

a. the RFP data package (reference 2-2) 

b. conversation with MSFC and others 

c. best estimates. 

2.1 PHYSICAL CONFIGURATION 

The Space Construction Base concept consists of a 
baseline configuration shown in Figure 2-1. This 
configuration, which is later expanded, includes a 
Habitability Module (crew quarters), Subsystem Module 
(control center). Turret Assembly (rotating member 
and structural support for the solar wings), Solar 
Wings A and B, and a Docking Module. The docking 
mo iule plays a major role in the buildup to a larger 
structure since, in many configurations, it becomes 
the attachment point for the next module. 

The space base modules are transported to orbit by the 
Shuttle Orbiter. The Orbiter Remote Manipulator System 
ext ■ .cts the payloads from the orbiter payload bay and 
. 'ds in the docking (berthing) and connect on of the 
modules . 



I 

L 
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TABLE 2-1 
MODULE ACRONYMS 


AP 

- 

Applications Pallet 

BFM 

• 

Beam Fabrication Module 

CCM 

- 

Construction Control Module 

DM 

- 

Docking Module 

HAB 

- 

Habitability Module 

HM 

- 

Habitability Module 

LM 

- 

Logistics Module 

MM 

- 

Manipulator Module 

MSM 

- 

Mission Support Module 

PSP 

- 

Public Service Platform 

RP 

- 

Research Pallet 

SM 

- 

SubsystKU Module 

SPM 

- 

Space Processing Module 

SW 

- 

Solar Wing 

SWA 

- 

Solar Wing Array 

TA 

- 

Turret Assembly 


OlTr.lNA^. ^ ^ 


2-14 



Figures 2-2 through 2-12 identify the buildup beyond 
the initial configuration to the final Space Base. 

It should be noted that for each configuration there is 
a period of time in which the Orbiter is either docked 
or connected to the base via the RMS. For these cases 
it should be assumed that the Orbiter Reaction Control 
System (RCS) provides the stabilization and control of 
the entire structure. Prior to attachment the Space 
Base/Orbiter configuration will be maneuvered to a 
minimum momentum buildup orientation ui.aed upon the 
selected attach point. This orientation will be 
accomplished by the Base and the Orbiter separately 
prior to docking, in order to minimize Orbiter RCS 
usage . 

2.1.1 Module Definition 

The basic modules associated with the space base are 
described briefly. Acronyms for each are listed in 
Table 2-1. 

° Subsystem Module - basically, the control center 
for the spacecraft operations. It houses all life 
support, communications and other electronic equip- 
ment. Fully habitable. 

° Habitability Module - Provides crew quarters, exer- 
cise facilities, health, sleeping and dining facilities . 

° Turret Assembly - The electromechanical structure for 
mounting the solar panels. Assume it has full freedom 
about the axis and can be servoed to a sun sensor 
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null. This structure also contains a pressurized 
passageway to interconnect adjacent modules. 

° Solar Wings A&B - Solar panel assemblies that will 
provide all space base power requirements until the 
beam building equipment is utilized (configuration 
10). Assume 50% occultation will provide adequate 
power . 

^ Docking Module - Simply, a structure that will allow 
a maximum of six separate modules to be attached in 
different directions. This module contains pres- 
surized passageways for crew operations. The docking 
ports on this module shall be considered to be flex- 
ible and shall be simulated by a hinge-spring joint. 
Axial rotation shall also be considered in the con- 
trol system. 

^ Mission Support Module - A module that contains 
additional equipment such as experimental hardware, 
consumables, spare subsystem equipment, repair shop, 
test facility. 

° Logistics Module - A pressurized habitable module 
containing consumable provisions and equipment for 
furthering the space base buildup. 

° Space Processing Module - A laboratory setup for the 
experimentation of manufacturing hardware and pharma- 
ceuticals in a zero gravity environment . 


.. . TT 1 
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Application Pallet - This module is essentially a 
storage unit for equipment necessary in furthering the 
function of the Space Processing Module. The equip- 
ment is externally fastened to a structure whose 
outline is cylindrical (to be compatible with the 
Orbiter). The pallet equipment must be retrieved by 
Extra Vehicular Activity (EVA). 

Research Pallet - Similar to the Applications Pallet 
but essentially holds equipment for further space 
experimentation. It could include a gimballed or 
floating structure for deep space telescope (high 
energy, UV) experiments. 

Construction Control Module - This module is defined 
as the control center for the construction phase of 
the large antenna array system. 

Manipulator - During the main stages of actual con- 
struction of structures in space the Orbiter supplies 
the space base with two structures identified as 
Manipulators. These manipulators are viewed as mecha-. 
nical arms that have full rotational freedom about the 
Xy axis and 90° of "L-BOW" motion. The Xy rotational 
freedom is supplied by two identical turret assemblies 
like those supporting the solar wings. These Mani- 
pulators are controlled through the Construction 
Control Module. 
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I ° External Tanks - During the latter stages of base 

buildup two Orbiter booster tanks are added to the 
configuration. These tanks are empty, uninhabitable 
and serve only as strongbacks (support structures) 
for the construction of beams. 

^ Beam Fabrication Module - The Beam Fabrication Module 
contains the raw materials and machinery for the beam 
manufacturing. It is assumed that the space base crew 

* 

I operates the equipment from within the module and 

\ stores th sections externally, on the external tanks 

; until final assembly begins. 

I ° 150KW Arrays - These are large solar arrays that will 

provide the power necessary to operate the beam con- 
struction equi )ment and ultimately the final space 
base configuration. Assume the structure has unlimited 
freedom about Zy and can be sun servoed. 

■ o 

Public Service Platform - The Public Service Platform 
Is the antenna conmunication array along with other 
equipment capable of monitoring items such as hurri- 
canes, earthquakes, weather traffic control, personal 
communication, etc. 

2.1.2 Mass Properties 

The coordinate system to be used for analysis and 
definition purposes is shown in Figure 2-13. It is 
right-handed, with the origin being placed at the Sub- 
system Module/Habitability Module interface. 
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TABLE 2-2 


SPACE COMmUCTiCOK AASB HASS/OIMEnd]3N 51PMAAV 


mm gsL.m 

L tf T 


LOCATI OM /fHCHESl 

K r 


1 ' i 

4^ A 


SUSiYSIGH tIODULE 

U,440 

ll’ K 

U.5*D 


6S 

C 

0 

^JASlTAfilLm' MODOU 

16.4S7 

26' X 

14,5'D 


-156 

C 

0 

TURRer AS:»i!DIBLY 

6S0 

ir X 

8.75*D 


' 158 

C 

i 

SOLAR VINO A 

SSI 

104* X 

26* X 


198 

C 

715 

SOLAR WING B 

SSI 

104* X 

26* X 


158 

w 

-715 

OOCXIMQ MODULE 

3,S2S 

lO.S* X 

U.7* X 

U.7’ 

-377 

0 

3 

MISSION SUPPORT MODUU 

13,697 

20.8* X 

14.5*0 


389 

0 

0 

DOCKING MODULE 

5,625 

lO.S* X 

U.7* X 

U.7* 

578 

0 

0 

LOGISTICS MODULE 

20,575 

17.5* X 

14.5'D 


•346 

0 

0 

SPACE PROCESSING MODULE 

25,794 

51.5* X 

14.5*0 


932 

c 

0 

Ai>l>tICATIONS PALLET 

16,400 

25* X 

14 *D 


-377 

0 

220 

RESEARCH PALLET 

16,400 

25* X 

14 *D 


-377 

0 

-220' 

LARCH LOGISTICS MODULE 

20,375 

17.5* X 

14.5'D 


578 

175 

0 

RESZARCIi LOGISTICS NODULE 

20,575 

17.5* X 

14.5*0 


578 

-175 

0 

HABITABILITY MODULE 

16.4S7 

26* X 

14.5*0 


-377 

358 

0 

SUBY3T3IS MODULE 

16,440 

11* X 

14.3'D 


-377 

136 

0 

HASITAaXLITY MODULE 

16,487 

26* X 

14.3*0 


-377 

-358 

0 

SUBSYSTEMS MODULE 

16,440 

11’ X 

14.5*0 


-377 

-136 

0 

COHSTRUCTION CONTROL MODULE 

13,670 

17.5* X 

14.5*D 


1366 

C 

0 

DOCKING MODULE 

3,625 

U.7* X 

10.8* X 

U.7* 

1541 

0 

250 

TURRET ASSEMBLY 

650 

8.73'D 

X ll'A 


1541 

C 

125 

DOCKING MOOUT.E 

5,625 

11.7* X 

10.8* X 

U.7* 

1541 

0 

0 

TURRET ASSEMBLY 

630 

8.75'D 

X 11* 


1541 

0 

-125. 

DOCtCING MODULE 

5,625 

11.7* X 

10.8* X 

U.7* 

1541 

0 

-250 

MANIPULATOR 

SOS 

43* X 

1.5’D 


1360 

0 

125 

ilANIPULATOR 

SOS 

43* X 

l.5*D 


1723 

0 

-125 

(EXTERNAL TANK A) 

(75,000) ( 

(27.6*0 X 

134*) 


1541 

0 

-1237 

SEAM FABRICATION MODULE A 

33,020 

15*D X 

13’ 


1791 


-250 

(EXTERNAL TANK 3) 

(73,000) < 

:27.6'D X 

154) 


1541 

0 

1237 


^ I I 

It A -« 


*WK2N DIAMETER APPEARS FIRST, IT ZMPLIIS DIAMETERS LIS IN PLANE X 70 X>Y PLANE 


'%.V- 




^ It 
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TABLE 2-2 (CONTINUED) 



MASatf) 

saa.m 

L « T 


LOCAiicxr 
X Y 

z 

BEAM fASRlCATtOH MODULI B 
(“ISO Kw" arbay vmiG A) 
("150 WJ” ARRAY WING B) 
("ISC KW" SUB8Y8TBIS) 

35,020 

(8.276) 

(5.276) 
(45,613*) 

15'0 X 13' 
(82 X ISO X 
(82* X 190* X 
(5'D X 5*> 

21') 

21’) 

1791 

-1143 

-1143 

-1143 

0 

0 

0 

0 

1294 

-1294 

0 

CONSTRUCTION LOOIITICS NODUU 

20,575 

17.3' X 14.3'D 


1716 

0 

0 

(PUBLIC BBRVICI PUTPOUl A) 

(21.600) 

(103' X 440* X 

312') 

1889 

0 

0 


«B)KRGY STORAGB IN C ONSHWCTZCK COimOL MOIHILI (ASSttlPTIOlO 


TABLE 2-3 

ELEMENTS OF THE INERTIA MATRIX 
o GIVEN IN VEHICLE COORDINATES (XV, YV, ZV) 


^XY 


'YX' 


'YZ 


- I 


ZY' 


'ZX 


• I 


x:. 



WEIGHT 

ON-AXIS (DIAGONAL) TERM 

PRODUCT 

OF INERTIA TERMS 

Dl 

(LBS) 

(SLUG-FT^) 


(SLUG-FT^) 



W 



^XX 

lyY 

^ZZ 

^XY 

BBi 

^ZX 

1 



gmnin 


0 


6 

2 

80801. 

-- - 

.29815E6 

.27208E7 

.24854E7 

0 

0 

0 

3 

106595. 

. 31920E6 

.73281E7 

.70927E7 

0 

0 

0 

4 

139395. 

.68683E6 

.9S917E7 

.88138E7 

0 

0 

0 

5 

180545. 

.99243E6 

.11383E8 

.11077E8 

0 

0 

0 

6 

246399. 

.20896E7 

. 14562E8 

. 15299E8 

0 

0 

0 

7 

260069. 

.21008E7 

. 19630E8 

.20367E8 

0 

0 

0 

8 

354854. 

.2 20 3 IE 8 

.75525E8 ^ 

. 56566E8 

0 ^ 

! 

0 

+.21268E8 

8A 

354854. 

.16802E8 

.88217E8 

.74487E8 

0 

0 

0 

9 1 

_| 

464874. 

.52726E8 

.13463E9 

.85231E8 

0 

0 

+.19761E7 

10 

562061. 

.61369E8 

. 19729E9 

.13928E9 

0 

0 

+.588c7E7 

11 

582636. 

.61393E8 

.20261E9 

. 14459E9 

0 

0 

+.56941E7 

12 

588637. 

. 70445E8 

. 20090E9 

.20986E9 

-.34308E7 

+.10579E7 

+.59784E7 
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Table 2-2 lists the module dimensions, weight and 
center of mass for each component. More detailed moment 
of inertia properties for each module is given in 
reference 2-2. 

Moments and products of inertia and the weight of each 
configuration are listed in '’’able 2-3. Reference 2-2 is 
the source of more detailed configuration information. 

It lists the configuration breakdown summaries, including 
radius of gyration (K), total weight, configuration CG 
with respect to the coordinate system origin, vehicle 
axis moments of inertia, principal axis moments of 
inertia, direction cosines and angles between the ve- 
hicle axes and principal axes. 

2 , 2 ORBITAL PARAMETERS 

The initial Space Base configuration will be launched 

into a 270 nautical mile orbit of either 28.5 or 55 

degree inclination. The entire mission of buildup and 

operation will be in Low Earth Orbit (LEO). Orbit decay 

may require periodic reboost but this problem shall not 

be addressed in this study. The orbit will be assumed 

to be circular. If we then assume a spherical earth, 

the orbital period will be 5668.2 seconds. The corres- 

-3 

ponding orbital rate will be 1.10850x10 rad/sec. The 
linear velocity for this altitude is approximately 
24,975 ft/sec. 
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2.3 OPERATIONAL REQUIREMENTS 

Operational Requirements of the Space Construction Base 
were non-existent in the data package compiling the 
study contract. In order to determine a feasible con- 
trol system for the SCB, certain assumptions will be 
made regarding its operation and mission accuracies. 

All of these assumptions are listed hero. As the study 
progresses, these parameters may be altered based upon 
new information or driven by other constraints. 

2.3.1 Mission Timeline 

The period of buildup and operation of the Space Base is 
shown in Figures 2-14 a and b. This timeline is hypo- 
thetical but can be used as a baseline for consideration 
in determining duration of each configuration in irbil 
and the orbital maintenance and lifetime of subsystem 
components. The lifetime of configuration 12 can be 
eventually considered open-ended for purposes of the 
study. Note that each configuration is expected to have 
a lifetime of approximately 4 months. 

Along witii the timeline presented, the crew is c.'cpected 
to grow in numbers as shown in Table 2-4. 

2.3.2 Attii-.ude Pointing and Maneuvering 

Attitude pointing and maximum slewing requirements are 
listed in Table 2-5 as a function of conf igurati... 
number. The maximu.m slewing rates are reduced for later 
configurations to bound the angular momentum require- 
ments . 
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TABLE 2-4 


CREW REQUIREMENTS 
FOR THE SCB MISSION 

CONFIGURATION CREW 


1 3 

2 3 

3 6 

4 6 

5 9 

6 9 

7 9 

8 9 

9 9 

10 12 

11 12 

12 12 
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POINTING AND MANEUVERING REQUIREMENTS FOR THE 
SPACE CONSTRUCTION BASE 
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TABLE 2-6 


GENERAL REQUIREMENTS 
FOR THE SCB MISSION 


Both turret assemblies have full freedom 
about Zy and can be servoed to the sunline. 

Assume 50% of solar panel area can be occulted 
without loss of required power to spacecraft 
functions. 

There is no requirement for periodic rotation 
of the spacecraft about any axis for the 
purpose of thermal control (no "barbecue" mode). 

Every docking module joint is to be considered 
"flexible" (simulated by a spring-hinge). 

Solar panel bending modes are estimated by 
using the results of a study of the SEPS 
Solar Array Dynamic Analysis (May 1976). 

It is most likely that docking will be refined 
to "berthing" in which tho impulse imparted 
to the spacecraft will be less than the present 
docking forces. However, we will use docking 
data for analyses since it will be worst 
case data. 



i 


TABLE 2-7 

CONFIGURATION PECULIAR REQUIREMENTS 
FOR THE SCB MISSION 

o IN ADDITION TO GENERAL REQUIREMENTS 


CONFIGURATION 

REQUIREMENTS 

1 

a) Entire configuration is prefab and delivered 
to orbit by the shuttle orbiter. Solar panels 
unfold. 

b) Turret assembly joint to spacecraft main body 
should be considered to be identical to a 
docking module joint. 

2 

a) No peculiar requirements. 

3 

a) Space processing module added; crew may 
double in size (from 3 to 6 men). 

b) Assume that electromechanical equipment can 
impart impulses to the spacecraft of 100 
times a typical crew motion (wall push-off). 

4 

a) Assume Research Pallet contains a separately 
stabilized experiment. 

b) This is the first configuration to require 
inertial pointing capability. 

c) Must consider solar panel occultation by 
Research and Application Pallets. 

5 

6 

a) First configuration with modules placed in 
YyZy plane. 

a) Additional modules in YyZy plane. 

b) Assume additional crew members (max of nine). 
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TABLE 2-7 (Cont'd) 


CONFIGURATION 


REQUIREMENTS 


a) No additional special requirements. 

a) Major imbalance in vehicle axis inertias. 

b) Tank is a strongback. 

c) Consider an alternate configuration where the 
single tank is placed along T-Xy direction 
until the second tank is added. 

d) Building manipulator arms appear. 


a) Second tank is added. 

b) Solar panel occultation may be significant 
to alter the main vehicle orientation. 


a) Large solar arrays are added. 

b) Assume same articulation as smaller solar 
arrays. 

c) Assume large solar arrays supplement smaller 
units. 

d) Assume large array connection to main space- 
craft body is a spring-hinge joint. 


a) Beam construction begins for large antenna 
structure of 12. 





TABLE 2-7 (Cont’d) 


CONFIGURATION 


REQUIREMENTS 

12 

a) 

Antenna structure is completed. 


b) 

Subsystems are added for information gathering 
and transmission. 


c) 

Assume antenna structure is earth pointing 
continuously. 


d) 

Examine a configuration where the X principal 
axis of inertia is the controlled axis. 


e) 

The combination of the antenna structure 
and subsystems module equipment form the 
"Public Service Platform" . 
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The time duration for inertial hold modes during con- 
figurations 4, 5, 6 and 7 may be limited by the 
capability of selected momentum exchange actuators and 
the desaturation system. 

2.3.3 General Requirements 

The general requirements for the SCB mission are listed 
in Table 2-6. It should be noted that these general 
requirements apply to all configurations unless other- 
wise stated. 

2.3.4 Configuration Peculiar Requirements 

In addition to the general requirements, Table 2-7 
lists requirements which are peculiar to individual 
conf igur at ions . 

2.4 REFERENCES 

2-1 Jennings, J. , MT-40,801, Definition of 
Space Base Buildup, 21 October 1977 
2-2 Space Construction Base Buildup Summary, 

Data Package Attachment to Exhibit A 
Scope of Work, Request for Quotation 
1-7-ED-07552-AP131D, Marshall Space 
Flight Center, 10 March 1977 
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SECTION 3 


ACTUATOR SIZING 

It is presumed that the SCB actuator system is composed 
of momentum exchange devices; in particular, double 
gimbal control moment gyros (DGCMG’s). One inherent 
advantage is the spherical angular ir.omentum envelope. 
Modified SKYLAB-type DGCMG's are assumed: 2300 ft-lb-sec 

each and with gimbal stops removed. The number of DGCMG 
units reouired will be a function of a momentum profile 
based on the magnitude of disturbance torques to which 
the SCB is subjected. In general, this CMG requirement 
is based upon the long term cyclic variation in angular 
momentum disturbances. The angular momentum buildup 
due to bias disturbance torques will be used to size and 
select momentum desaturation schemes in Section 4 for the 
various SCB configurations. 

Both long and short term disturbances on the SCB will be 
defined in the following paragraphs, with CMG sizing 
being primarily a function of the long term effects. 

LONG TERM DISTURBANCES 

Most of the long term angular momentum variations are 
caused by torques due to gravity gradient and aerodynamic 
effects on the SCB. These are described in references 
3-1 and 3-2, respectively. Other less important sources 
are magnetic and radiation torques, which are described 
in references 3-3 and 3-4. Each will be described and 
then combined into a composite angular momentum require- 
ment for each configuration. 


Gravity Gradient 

The largest torques, which a spacecraft is subjected at 
the altitude of interest, are ordinarily those due to 
forces from earth gravitational effects. Various space 
base orientations are also examined to determine those 
which are more practical with respect to minimizing the 
gravity gradient angular momentum requirements. These 
gravity gradient torques are obtained for all twelve 
space base configurations for the following orientations: 

1. X-POP (Xy perpendicular to the orbital plane) 

2. Y-POP 
.3. Z-POP 

4. Yy Solar Inertial (wo^st case g = 45°) 

The last listed orientation was included to obtain a 
worst case inertial pointing requirement for space base 
configurations 4, 5, 6 and 7. Pointing experiments from 
the Application Pallet or the Research Pallet will 
probably be most expedient at these times. From con- 
figurations 8 and on, the space base becomes cumbersome 
and the many large appendages may occult the field of 
view of any pointing experiments. It is assumed that 
for these experiments, the +Zy axes may be called upon 
to point anywhere in the celestial sphere; and item (4) 
above is a worst case orientation with respect to 
gravity gradient bias torques. 
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The geometry of the space base in earth orbit is shown 
in Figure 3-1; the orientation shown is Yy solar Inertial. 
Some of the assumptions used in computing gravity 
gradient torques were the following: 

1. Orbital altitude - 270 N.mi. 

2. Orbital inclinations up to 55 degrees 

3. Spherica’ earth (non-oblate) 

4. Circular orbit 

One other assumption is used for determining feasibile 
space base orientations: 

Solar wings and the 150 KW Solar Arrays are rotatable 
about the Zy axis. 

Space base mass data is available in Section 2. It is 
presented in the following form: 

1. Moments of inertia about the three vehicle axes 
(Xy, Yy^ Zy). 

2. Principal moments of inertia. 

3. Transformation pCy from vehicle to principal axis. 

Vehicle moment of inertia matrix ly is obtained by 
using the following similarity transformation: 

ly - p'^v’’ 'p p=V 

where each term is a 3x3 matrix, and Ip is a diagonal 
matrix of the principal moments of inertia. The resulting 
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moment of inertia matrix wi‘‘.l be: 



- 




^XX 

^XY 

^XZ 

'v - 

^YX 

Iyy 

^YZ 


^ZX 

^ZY 

^ZZ 

where 

^YX 

= 

^XY 


^ZY 

= 

^YZ 


^XZ 


^ZX 


The gravity gradient torques about the vehicle axes 
were calculated using the following expressio. s from 
reference 3-5: 


gx 


gy 


SI".. - 


I )a a +I (a^ -*a^) + I a a -I a a*] 

\y\T ' \7 rr ry f v'y V xr vxr V <2| 


yy' y z yz^ y z' xz x y xy x z 


-^^*‘yxVz-'yzVx| (3- 


j * '^1(1 “I )a a +I (a^ -a^) + I a a 
gz L yy XX ' X y xy ' X y’ zy z x 


y“x] 

-'zx«z“y] 


where 

16 3 2 

k is the gravitational constant * 1.407654x10 ft /sec 
R is the distance from the planet mass center to the 
spacecraft mass center = 


I. . with i,j * X, y, z are the components of the inertial 
dyadic 





.H. 


1 1 
Hi'- * 


> t 

: I 

m i 





.■ I 

. t 


with i = x,y,z are the direction cosines of the 
vector R with respect to the X^, Y^, coordinate 
frame (see Figure 3-1) 


and 

= earth's radius = 3440.0756 N.mi (mean) 
h = 270 N.mi. 


The magnitude of R is thus obtained, 

R = R +h = 3710.076 N.mi = 2 . 254285xl0^ft 
o 

The coefficient outside the bracket of (3-3)can also be 
defined as: 


^ - 
r3 " 



(3-4) 


where 


= orbital rate = 1.10850x10 rad/sec 


I 

* a 


The orbital period T is also obtained and calculated 

o 

3S 5668.2 seconds. 







The equations of (3-3) were implemented on the digital 
computer for X-POP, Y-POP and Z-POP orientations and were 
also integrated as a function of time to obtain angular 
momentum changes due to gravity gradient torques. The 
angular momentum was then expressed in cyclic and bias 
buildup terms. 

An additional configuration was added by eliminating 
the asymmetry of the single external tank of configuration 
8. The tank was temporarily repositioned along the +X^ 
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TABUS 3-1 



L I ANGULAR MOMENTUM REQUIREMENTS 

* DUE TO GRAVITY GRADIENT TORQUES 

I ^ 

ii O ORIENTATION: LOCAL VERTICAL WITH 0 = 45° 

I o ALTITUDE: 270 N. Mi. 

I s 
« '4 

I I o AT EACH INSTANT OF TIME, 

* ^ MAXIMUM MAGNITUDE IS PRESENTED 


BUILDUP H CYCLIC H 

(FT-LB-SEC/ORBIT) (FT-LB-SEC) 


1 

1 

«z 

«T 

1 

1 

sm 

«T 

3017. 

0 

0 

3017. 

240.0 

4778. 

5117. 

7001. 

1598. 

0 

0 

1598 ■ 

127.2 

5929. 

6109. 

8513. 

3852. 

0 

0 

3852. 

306.6 

7766. 

7332. 

10680. 

3852. 

0 

1 

0 

1 

3852. 

306.5 

i 

10735. 

10305. 

14884. 


tm ft , 














axis by using the middle docking module. The mass 
characteristics were recomputed and the symmetric structure 
is designated as configuration 8A. All of these results 
are tabulated in reference 3-1 and the X-POP orientation 
was selected for configuration 1 through 11 . It is 
presently planned that configuration 12 will operate in 
the X^ local vertical orientation with the microwave 
antennas pointing towards earth. 

The angular momentum magnitudes for the selected 
orientations are summarized in 3.1.6. Angular momentum 
magnitudes are shown in Table 3-1 for a solar inertial 
orientation for configurations 4, 5, 6 and 7. The 3 
angle was a worst case of 45°. This, of course, is not 
a preferred orientation. If the sizing of the CMC 
system is based upon a preferred, low torque orientation, 
its capacity may not be quite enough to handle this 
worst case solar orientation - - except for a very 
limited time. 

3.1.2 Aerodynamic 

Aerodynamic forces, torques about the c.m. and angular 
momentum are obtained for each for all twelve orbital 
configurations. This effort was concentrated on those 
vehicle modules having large surface areas at some dis- 
tance from the overall center of mass (c.m.). The 
portions of the vehicle in this classitication are the 
following: the solar wings, the solar arrays and the 

external tanks. The tank(s) are included in configurations 
10, 11 and 12. 
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Worst case aerodynamic torques are obtained using the 
"free molecular flow" theory which is recommended fcr 
altitudes above 65 N.mi. (120 KM). Normal forces on 
flat areas (solar wings and arrays) are then obtained. 
Angles of attack on flat surfaces, and "shading" effects 
are much more complex and are not included in this report. 

In an atmospheric environment, the aerodynamic force 
would be 

F = q Cp A (3-5) 

where 

F = the aerodynamic force 

C 

D = drag coefficient for the body of the interest 

A = projected area of the space base element normal to 

the incident flow. 

2 

q = i pV = dynamic pressure 
p = mass density of the atmosphere 

V = space base velocity, or relative velocity of the 

space base element relative to the local atmosphere 

Reference 3-6 includes a table for determining the Cj^ 
for several basic body shapes in a hyperthermal free 
molecular flow condition. is also a function of the 

speed of reemitted molecules and whether a diffused 
reemission or a specular reflection occurs. The worst 
case was taken and is conservatively specified as 
being 2.6 for a flat plate and 2.5 for a cylinder. 
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The flat plate areas for the two solar wings and the two 

2 2 

solar arrays are 5,408 ft and 31,160 ft , respectively. 

2 

The cylinder side area per external tank is 4250.4 ft 

each. For a 270 N.mi (5Co Km) orbit, the space base 

velocity is 24,975 ft /sec, and the worst case mass 

-14 3 

density (reference 3-7) is 1.53 x 10 slugs/ft . A 
listing of mass densities as a function of night and day 
and also high/low solar activity is given in Table 3-2. 

Using (3-5) for the worst case atmospheric density, 

F * (4.772 X 10"® Ib/ft^) Cj^A 

For flat plates, 

F = (1.241 X 10"® Ib/ft^) A 

For cylinders, 

F = (1.193 X 10"®lb/ft^)A 

Using these expressions, 

F = 0.0671 lb (solar wings) 

F = 0.3866 lb (solar arrays) 

F = 0.0507 lb (each external tank) 

Both the solar wings and solar arrays were assumed to be 
rotatable around the Zy axis. With the panels facing 
towards the forward velocity and the initial Yy axis, 
the moments would be as follows: 

Ljjy = 0 (because of balancing aero forces, except for 
configuration 8) 


‘T V 
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TABLE 3-2 

ATMOSPHERIC MASS DENSITY 

• ALTITUDE OF 270 Naut. Hi (500 KM) 

• FROM NASA 8P-8021, MODEL OF THE EARTH'S ATMOSPHERE 


CONDITIONS 

ATMOSPHERIC 

MASS DENSITY 

SLUGS/ FT^ 

NORMALIZED 
TO MAXIMUM 

HIGH 

NIGHT 

1.53 X 10 

1 

SOLAR ACTIVITY 

__DAY 

9.70 X 10"^® 

1/1.58 

LOW 

NIGHT 

2.44 X 10“^® 

1/62.7 

SOLAR ACTIVITY 

DAY 

6.40 X 10“^^ 

1/239 




LyV - 0 

Lzv “ very high 


* 4i 


With all solar panels facing towards the original +Xy 
axis with edges into the relative "wind", 


4v-o (except for configuration 8) 


The effort here w^ll therefore concentrate on the former 
case where becomes very high. In general, 


-L^v ■ X (X^p - Xj^jj) 


(3-6) 


where 

Fyy = each vector force parallel to the Yy axis 
^CP = each aerodynamic center of pressure location in 
the Xy direction 

XcM ~ vehicle center of mass location along the Xy 


The angular momentum can then be computed as follows: 


®ZV “ ^ ^ZV 


(3-7) 


For a solar inertial condition, the L^y torque will be 
cyclic and 


Hzy(PEAK) =(T^/4 tt) Lzy(PEAK) 


(3-8) 
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where 


= orbital period = 5668.2 sec 

For an Xy local vertical condition the momentum on the 
solar wings and arrays will be as in (3-8); but for the 
external tanks: 

H^y (BUILD UP PER ORBIT) = T^L2y(BIAS) (3-9) 

The assumptions for these aerodynamic characteristics 
are reiterated as follows: 

(a) Only space base modules with large surface areas 
at a large moment arm from the c.m. are considered 

(b) No "shading" of the airstream by one section of 
the vehicle on another is considered 

(c) Aerodynamic force computations are based on a 
space base velocity along its Yy axis (assumes the 
sun is also along this axis) 

(d) Surface drag coefficients used were for the worst 
case hyperthemal free molecular flow 

(e) Worst case atmospheric density was assumed (high 
solar activity at nighttime). 

With regard to this last item, it must be pointed out 
that the atmospheric mass density (p) can also be 
approximately 240 times less than that used in the 
calculations of aerodynamic force. For low solar 
activity at daytime, all computed forces, torques and 
angular momentum magnitudes would be reduced by a 
factor of about 240. 
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The (X^p - X^jj) radius arms, aerodynamic forces, Lj,y 
torques and angi’.lar momentum levels for the solar 
wings, solar arrays and the external tanks are tabulated 
in reference 3-2. The composite aerodynamic angular 
momentum magnitudes are summarized in section 3.1.6. 

The results are given for the worst case atmospheric 
density (high solar activity at night) and also when 
normal (actually minimum when low solar activity and in 
daytime) . 

3.1.3 Magnetic 

The peak cyclic magnetic torques are obtained using the 
following vector expression: 

L = ii X B (3-10) 

where 

L = torque vector (N-m) 

2 

VI = spacecraft magnetic moment vector (A-m ) 

2 

B = peak of the earth's flux density (tesla or webers/m ) 
above the magnetic poles 

Reference (3-8) presents an approach for estimating y, 
the magnitude of a spacecraft's magnetic moment, on the 
basis of two factors : 

1. The degree of control over current loops, magnetic 
materials, quality control of potential moment 
sources, test and magnetic compensation. 



\ 


i 
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2. The size of the vehicle, specifically its mass. 

Minimal control of (1) above is designated Class III. 

An estimate of the dipole magnetic moment per unit mass 

2 

for a Class III spacecraft would be approximately O.OlA-m /KG. 

A model of the earth's magnetic field (B) in orbit is 

available in reference (3-9). The magnetic flux is 

given in terms of a spherical harmonic expansion series. 

Only the first term or so is required to obtain an estimate 

of the maximum magnetic torques. At an orbit altitude 

of 270 N.mi., the magnitude over the magnetic poles is 

-5 

approximately 0.60 gauss or 6 x 10 tesla in the vertical 

direction. The magnitude at the magnetic equator is 

—5 

0.30 gauss or 3 X 10 tesla, and the direction of the 
flux is horizontal to the earth's surface towards the 
south magnetic pole. An orbit parallel to the plane 
of the magnetic equator would therefore result in a 
constant spacecraft torque for an inertially oriented 
vehicle. A higher torque is possible over the magnetic 
poles, but the result is cyclic. 

The peak magnetic torque is calculated above the magnetic 
poles for each space base configuration for the following 
assumptions : 

a. The orbit passes over the magnetic poles 

2 

b. The space base magnetic moment magnitude is O.OlA-m /KG; 
and the direction is always normal to B in the magnetic 
polar region. 
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TABLE 3-3 

MAGNETIC TORQUES AND ANGULAR MOMENTUM 
FOR THE SPACE BASE 

• SPACE BASE MAGNETIC MOMENT IS BASED ON CLASS III OF NASA 
SP-8018 (SPACECRAFT MAGNETIC TORQUES) y - .OlA-m^/KG 

• EARTH'S MAGNETIC FIELD IN MAGNETIC POLAR REGION AT 270 N.MI. 
ALTITUDE IS .60 GAUSS, OR 6 X 10“® TESLA 

• SPACE BASE IS INERTIALLY ORIENTED IN AN ORBIT OVER THE MAGNETIC 


POLES 

PEAK 

MAGNETIC TORQUE 
L = y X B 

PEAK 

ANGULAR 

MOMENTUM 
^ ^CYCLIC 

CONFIG. 

WEIGHT 

W 

MASS 

M 

MAGNETIC 
MOMENT u 

NO. 

(LBS) 

(KG) 

(A-M^) 

(N-M) 

(FT-LB) 

(FT- LB-SEC) 

1 

40904 

18554 

185.5 

.0113 

.00821 

3.70 

2 

80801 

36651 

366.5 

.0220 

.0162 

7.32 

3 

106595 

48351 

483.5 

.0290 

.0214 

9.65 

4 

139395 

63228 

632.3 

.0379 

.0280 

12.62 

5 

180545 

81894 

818.9 

.0491 

.0362 

16.35 

6 

246399 

111765 

1117.6 

.0671 

.0495 

22.31 

7 

260069 

117965 

1179.6 

.0708 

.0522 

23.55 

8 

354854 

160959 

1609.6 

.0966 

.0712 

32.13 

9 

464874 

210863 

2108.6 

.1265 

.0933 

42.09 

10 

572061 

259482 

2594 . 8 

.1557 

.1148 

51.80 

11 

582636 

264279 

2642.8 

.1586 

.1170 

52.75 

12 

.588637 

267001 

2670.0 

. 1602 

.1182 

53.30 



The calculated cyclic magnetic torques and angular 
momentum are listed in Table 3-3. 

For the unique condition of having an orbit parallel 
to the magnetic equator, the steady magnetic torque can 
be calculated by multiplying each torque item in Table 
3-3 by one-half. The magnetic angular momentum buildup 
per orbit can be obtained by multiplying each H^ycLIC 
in Table 3-3 by 2 it. This was done in summary tabulations 
given in section 3.1.6. 

The results indicate that magnetic contributions are 
negligible and will not have to be considered for 
simulation modeling or for sizing attitude control actua- 
tors and momentum exchange systems. If a magnetic 
desaturation system is used on the space base, any 
residual magnetic moment dipole will be automatically 
compensated . 

Radiation 

Sources of electromagnetic radiation that cause forces 
and possible torques to act on the space base (reference 
3-10) are the following: 

1. Direct solar photon radiation 

2. Solar radiation reflected by the earth and its 
atmosphere 

3. Radiation directly from ti>e earth and its atmosphere 

4. Radiation from portions of the space base itself. 
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The most important cause of radiation is (1) above. 

The intensity of the radiation, or solar flux, is given 
In terms of energy per unit time through a cross-sectional 
unit area: 

« 1353 w/m^ +3.3% 

The percentage is not a tolerance, but a seasonal variation 

because of the earth's orbit about the sun. The worst 

2 

case magnitude is then = 1398 w/m . The forces caused 
by other sources are usually at least an order of 
magnitude smaller. 

Radiation due to the earth's reflectance is maximum at 
the subsolar point (high noon). It reduces to zero on 
the night side. The peak magnitude when at 270 N.mi. 

(500 Km) and assuming the earth's planetary albedo of 
0.34 is 

Ig * 600 w/m^ (peak) 

The mean magnitude of I„ for a 50% daylight orbit is 

2 ^ 

approximately 191 w/m . 

Radiation directly from the earth and its atmosphere at 

2 

sea level has a mean emission intensity of 243 w/m . 

With clear skies and at the lower latitudes, the worst 
case Intensity at a 270 N.mi. orbit is 

I 2 = 150 w/m^. 
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Radiations and reflections from the space base itself, 
item (4), are the least important and will be assumed 
as negligible. 

For each of the (1) through (3) radiation sources, the 
radiation force is calculated as a function of whether 
the radiation is 

a. Completely absorbed 

b. Completely specularly reflected 

c. Completely diffusely reflected. 

The worst case turns out to be when the surface is such 
that a complete specular reflection is obtained. The 
expression defining this force (reference 3-10) is 

F = i [-(l+c^„)cosen+(l-c^^)sin8s 1 (cose)A (3-11) 

C ^ X S J/S aJ 

where 

I = radiation intensity 

Q 

c * speed of light » 2.997925 x 10 m/s 

c * coefficient of reflectivity, assumed to be 1 for 
rs 

the worst case 

0 ® angle between the surface normal and the direction 
of radiation 
n * unit normal vector 
s * unit vector along the surface 
A » area of the surface 


« 

I 


i 
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(3-12) 

(3-13) 


Fjj = 2AI/C = normal force when 0 is zero 
Fjj » 6.671281 X 10"® AI 

Using the numerical value of flux intensity estimated 
for Ig and Ig, the normal forces are 

F^ = 9.3265E-6 A 

Fg = 4.0028E-6 A (peak) 

Fg = 1.0007E-6 A(worst case) 

2 

where F is in newtons and A is in m units. 

Fj^, Fg and Fg are computed for the surfaces of three 
vehicle modules: the solar wings, solar arrays and 

external tank(s). All other module surfaces are smaller 
and/or are a., a shorter moment arm. and are assumed to 
be negligible in computing torquing moments. 

Table 3-4 is a listing of maximum F^, Fg and Fg forces 
for each of the three modules. The numerical magnitudes 
are given in pounds. The areas for the solar wings, 
solar arrays and for one external tank are A , A and 

W Or 

A^, respectively: 

A = 2 X 2704 ft^ = 5408 ft^ 
w 

A = 2 X 15580 ft^ = 31160 ft^ 
a 

2 

A (projected) = 4250.4 ft (each) 



Table 3-4 
Radiation Forces 
on 

Large Surfaces of Space Base 
•All forces given in lbs 
• Forces normal to plane surfaces 
or projection of non-planar bodies 


Radiation 

Type 

N 

Solar 

Wings 

(s) 

Solar 

Arrays 

(a) 

External 
Tank (One) 
(t) 

External 

Tanks 

(tt) 

1 

1.0534E-3 

6.0696E-3 

8 . 2793E-4 

1.6559E-3 

2 

4.5211E-4 

2.6050E-3 

3.5533E-4 

7 . 1066E-4 

3 

1.1302E-4 

6.5122E-4 

8.8831E-5 

i . 7766E-4 


Fjj * 2.03341 X 10"^ AI (pounds) 

= 1398w/m^ = 95.794 Ib/ft-sec 

1 2 = 600w/m^ = 41.113 Ib/ft-sec 

1 3 = 150w/m^ = 10.278 Ib/ft-sec 
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2 

For two external tanks, » 8500.8 ft 

Radiation torques are computed for two vehicle or lent a 
tions: XPOP inertial and X local vertical. Intensity 

of radiations causing both bias and cyclic torques are 
given in Tables 3-5 and 3-6 for the two orientations. 
In some cases, these would modify the forces given 
in Table 3-4. 

Peak torques, L^, Lg and Lg are computed using moment 
arms available from Section 2. These are listed in 
Table 3-7 for each section of the space base. The 
moment arm in each case is where 

X^p = center of radiation pressure on the module of 
interest along the vehicle's X axis 
^CM * space base center of mass location along the X 
axis 


Bias and cycl ic angular momentum components can then be 
computed from these torques. For the XPOP inertial 
orientation, all maximum torques would be about the Z 
axis. The Z axis angular momentum components would then 
be as follows: 

"bias'’’,, (3-14) 

"CYCLU; ' [- 2/2 + L 3 ] <3-15) 

where Lj^, Lg.and Lg are torques for solar wings, solar 
arrays and external tanks as each are included in a 
vehicle configuration. 
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Table 3-5 

Intensity of Radiation 
for 

Inertial XPOP Orientation 

2 

• All units in w/m 

• Results in Z axis torques 


Radiation 

Magnitude Used for Computation: 

Source 

Bias Torques 

Cyclic Torques 


1398 for panels 
1398 for tanks 


^2 

191 for panels 
191 for tanks 

+300 for panels 
+300 for tanks 

^3 


+150 for panels 
+150 for tanks 


Table 3-6 

Intensity of Radiation 
for 

X Local Vertical Orientation 

• All units in w/m2 

• Results in maximum Z axis to. '^es 


Radiat ion 
Source 

Magnitude Used for Computation: 

Bias Torques 

Cyclic Torques 

^1 

1398 for panels (S»90°] 
1398 for tanks " 

+1398 for panels(B=0'^) 
+1398 for tanks " 

^2 

191 for panels 
0 for tanks 

+300 for panels 
0 for tanks 

^3 


+150 for panels 
0 for tanks 
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Table 3-7 

Radiation Forces and Angular Momentum 

• Orbit Period (T^) is 5668.2 seconds 

o 

for 270 N.mi. orbit 





Angular Momentum (ft-lb-sec) j 

T 

+21.8 

21.4 

1.7 

12.0 

2.4 
14.0 

8.4 

22.3 

38.3 
43.7 

32.2 

35.3 
34.5 

a a 2 

00 00 

66.1 

113 

107 

117 

112 

113 

3^ 

Jik 

+28.8 

28.2 

2.3 

15.9 

3.2 

18.5 

11.1 

29.4 

50.6 

57.7 

42.6 

46.7 
45.6 

llg 

rH iH 


091 

FA8 

142 

15| 

148 

1^ 

w 

+155 

153 

12.4 
86.0 

17.4 
1,9.7 

124.3 

-159 

-273 

-312 

-250 

-252 

-246 

-5696 

-5824 

-5788 

608 

ZAt+ 

+764 

832 

797 

807 

Peak Torques (ft-lb) 



+ 

.00254 

.00021 

.00143 

.00029 

.00166 

.0010 

-.00265 

-.00455 

-.00519 

-.00383 

-.00420 

-.00410 

ill 
• • • 

1 t 1 

+.00786 

.0135 

S gi S S 


+.01035 

.0102 

.00083 

.00573 

.00116 

.00664 

.00398 

-.0106 

-.0182 

-.0208 

-.0153 

-.0168 

-.0164 

-.379 

-.388 

-.386 

+ .0314 
.0539 

g s « s 

i!) u5 . lO 

O o o O 
• •• • * 


+.0241 

.0237 

.00193 

.0134 

.0027 

.0155 

.00027 

^.0247 

-.0424 

-.0484 

-.0357 

-.0392 

-.0382 

-.884 

-.904 

-.898 

+.0733 

.1255 

+.1187 

.1292 

.1237 

.1253 

1 

s 

'd,« 

+22.90 

22.48 

1.83 

12.68 

2.57 

14.69 

8.80 

-23.44 

-40.26 

-45.95 

-33.92 

-37.18 

-36.27 

-145.67 

-148.93 

-148.02 

+88.48 

151.59 

+71.66 

78. 

74.73 

75.65 

Conflg. 

1 

2 

3 

4 

5 

6 
7 

8 

8a 

9 

10 
11 
12 

Zl 

TT 

01 

to S 

9 

10 
11 
12 

1 

Module 

Surface 

1 1 s 

SOLAR 

ARRAY'S 

(a) 

EXT.TANK 

(t) 

EXTERNAL 

TANKS 

(tt) 






5* ^ 
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Table 3-8 
Angular Momentum 
Due to 

Radiation Forces 

Space Construction Base Surfaces Considered are 
Solar Wings, Solar Arrays and External Tank(s) 



XPOP INERTIAL 

CONFIG. 

®BIAS 

(FT-LB-SEC/ORBIT) 

^CYCLIC 

(FT-LB-SEC) 

1 

+155 

+7.00 

2 

153 

6.87 

3 

12.4 

0.56 

4 

86.0 

3.88 

5 

17.4 

0.78 

6 

99.7 

4.49 

7 

124. 

2.69 

8+ 

313. 

28.4 

8a 

535. 

48.8 

9 

453. 

48.5 

10 

-5114. 

305. 

11 

-5279. 

310. 

12 

-5227 

308. 


ANGULAR MOMENTUM 


X LOCAL VERTICAL 


^lAS 

(FT-LB-SEC/ORBIT 


+155 


153 

12.4 

86.0 

17.4 
99.7 

124. 


256. 

438. 

361. 

-5214. 

-5375. 

-5324. 


Hy are also possible because of asymmetric configuration, 
For XPOP: * 288., to “ +26.1 ft-lb-sec. For X 


+28.8 

28.3 

2.30 

15.9 

3.22 

18.5 

11.1 


95.5 

164. 

165. 
1213. 
1235. 
1230. 


A w a. AS — A ww « 

Local Vertical : 


“CYCLIC - 
288., 


CYCLIC “ ±87-8 ft-lb-sec, 


, 1 . 

I [■ 

4 k 
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For an X local vertical orientation, the external tanks 
will not contribute to a Z axis torque for radiations 
from the earth (Ig and Ig). The Z axis angular momentum 
components were based on vehicle and solar panel orien- 
tations as determined in reference 3-4 . The resulting 
components, therefore, are the following: 

■•bias - [h * 4 /”] - [I' 2 /’'JtaSKS 


where the first term includes contributions from solar 
wings, solar arrays and the tanks. Also 


H, 


CYCLIC ' t L 2/2 + L 3 ] 

- <To/2’ )[V2 * 4], 


(3-17) 


'TANKS 

where the first term again includes all contributions. 


The results of equations (3-14) through (3-17) are illu- 
strated in Table 3-8 for all space construction base 
configurations. These angular momentum magnitudes are 
included in the surmnaries given in Section 3.1.6. 

Although these momentum levels are less than gravity 
gradient and worst case aerodynamic angular momentums, 
radiation effects appear to be greater than those obtained 
from magnetic torques. For a sizing estimate of a momentum 
exchange or a momentum desaturation system, torques and 
angular momentum caused by radiation may be considered 
negligible. 
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3.1.5 Vehicle Roll During X Local Vertical Orientation 
During a local vertical orientation of an orbital 
vehicle, it is necessary to rotate the vehicle about its 
vertical axis if solar panels are articulated about a 
single axis. This provides the solar panels with what 
is effectively "two-axis pointing". 

An X local vertical vehicle is illustrated in Figure 3-2, 
each coordinate space is defined as follows: 


0 orbital coordinates, defines the orbit and 
are in the orbital plane 

inertial coordinates, where Xj points tov.’ard the sun 
V vehicle coordinates, where Xy is the roll axis 
S solar panel coordinates, where Yg is the axis normal 
to plane of the panels 


The angles shown are defined as follows: 


“y vehicle position in orbit, defined as zero when Xy 
is parallel with X^ 

rotation about defining the solar angle between 
X and X, 

O 1 

vehicle roll angle about X , where zero is defined 

X V 

when Y^ is parallel to Y^ 

S solar panel angle about Z. , where zero is defined 

z V 

when Y_ is parallel to Y„ 

s V 

Angular rotations from 0 to !_ or S are indie ted on the 
figure where an interim coordinate system is defined: 

U vehicle coordinate system anywhere in the orbit 
where “ 0 
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Transformations may then be generated between the S 
and coordinate systems as a function of the angles. 
In order for the solar panels to point towards the 
sun, the following must be true: 

Yg . Xj » 1 (3-18) 

Other relationships which also must be true are 



The roll angle 4> and the solar panel angle 6 was then 
determined in reference 3-11 using the functions of 
(3-18) and (3-19): 

tan* = (3-20) 

sin6 = -cosa cos8 (3-21) 

where the notation was simplified by eliminating the 
angle subscripts: 
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Equations (3-20) and (3-21) are plotted in Figure 3-3 
as a function of the position in orbit (a) and the 
angle between the ecliptic and whe orbital plane (B). 
For a circular orbit, is directly proportional to 
time. 


Observation of Figure 3-3 for very small B angles 
Indicates a potential problem. Theoretically, the 
vehicle must roll 180® about in a very short period 
of time. If the moment of inertia of the vehicle about 
i-i appreciable, extremely high changes of angular 
momentum will be required twice per orbit - once positive 
r.nd then negative. The angular momentum swing about 
the X axis of configuration 12 is listed in Table 3-9 
as a function of the B angle. The amplitude of the 
angular momentum is defined by the following: 


**x “ ^xx '^MAX 


(3-22) 


“^MAX 

(3-20); 


is obtained by taking the time derivative of 


^MAX ' 


(3-23) 


and * orbital rate for the 270 N.mi. orbit = 1.1085x10 
o 

rad/ sec . 


Table 3-9 indicates the consequence of having a single 
axis articulation of the solar panels when the B angle 
has a small magnitude; the requirement for extremely 









Table 3-9 


ANGULAR ROLL RATE 
AND MOMENTUM REQUIREMENTS 

• X local vertical orientation 

• ■ 7.044512x10^ Ft -Lb-Sec^ for configuration 12 

• An equivalent CMG unit has an angular momentum of 
2300 Ft- Lb-Sec 

• Assumes sun line is perfectly normal to the solar 
panels for all 3 angles 


±0 

^MAX 

MAGNITUDE OF 
MAX. ROLL RATE 

EQUIVALENT 
NUMBER OF 

(DEG) 

(RAD/ SEC) 

(DEG/ SEC) 

CMG UNITS 

0 

00 

00 

00 

1 

0.0635 

3.64 

1950. 

10 

0.00629 

0.360 

192.7 

20 

0.00305 

0.1745 

93.4 

30 

0.00192 

0.110 

58.8 

45 

0.0011085 

0.0635 

34.0 

51.79 

0.0008726 

0.050 

29.4 

60 

0.000640 

0.0367 

19.6 

70 

0.000400 

0.0229 

12.35 

80 

0.000196 

0.0112 

6.0 

90 

0.0 

0.0 

0.0 
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high magnitudes of angular momentum about the X axis. 
Solutions to this problem of small 3 angles, which 
should be considered are one of the following: 

1. Double axis articulation of the solar panei.s. 

or 

2. Set a fixed <j) = + 90 degrees with the Z axis per- 
pendicular to th<* orbit (see Figure 3-2) and with 6 being 
controlled by sun sensors to keep Y aligned as close 

o 

as possible to Xj (rotating once per orb-w) for sti^all 
B angles. 

The problem withl. is the complexity of the additional 
solar panel control hardware. The disadvantage of 2. 
is the loss of some solar power: e.g., :o 86.6% of 

nominal for 6 = 30 deg or 70.7% of nominal for 6 = 45 
deg. This power reduction can be minimized by adding 
dedicated control or momentum exchange devices for 
the X^ vehicle axis (reaction wheels) such that two 
axis pointing can be resumed for [1 angles less than 30 
degrees. 


3.1.6 Summary 

All long term angular momentum requirements are combined 
to determine: 

o A total momentum envelope for prospev'.tive momentum 
exchange devices 

o Maximum momentum builduj^s due to bias torque for 
sizing momentum desatu ttion devices in Section 4. 


V 


L. 
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Angular momentum requirements for gravity gradients, 
aerodynamic, magnetic and radiation torques were determined 
in sections 3.1.1 through 3.1.4, respectively. These 
are combined on a vehicle-axis basis to determine 
maximum components of total angular momentum vectors. 

A summary of cyclic angular momentum requirements are 
given in Table 3-10 for all twelve space construction 
base configurations. The recommended attitude orienta- 
tion for each c >nfiguration is given. Alternate, 
modified orientations are also listed for configura- 
tions 8 through 11 where 8M, 9M, lOM and IIM are rotated 
slightly about the YV axis (from the nominal XPOP 
attitude) until the XP axis (principal X axis) is 
normal to the orbital plane. 12M is an alternate to 
configuration 12 (XV local veicical) wheie the XP axis 
is always along the local vertical . Alphabetical symbols 
after nu”'erical magnitudes identify the source of the 
angular i.x>ment\im contribution: 

A Aerodynamic torque 
G Gravity gradient torque 
M Magnetic torque 
R Radiation torque 

Maximum components along each of the ’ehicle's axes are 
labelled as H^, Hy and H^. The magnetic angular momentum, 
which is computed as a direct .unction of the vehicle 
mass, is given as a spherical component since the 

ma.netic moment of the vehicle has no pro 'erred orienta- 
tion. 
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' * Table 3_io 

{ r Summary of Angular Momentum 

Due to Long Term Cyclic Torques 



• Max 

Total is 

for 

Period 

Of High Solar 

Activity 



ANGULAR MDMEMDM 


(FT-IB-SBC) 



MAX OCMPONEMS 

TOTALS 

CXUFIG. 

ATTnUDE 


i. 

«Z 

— 

_5xxz 

MAX 

NORMAL 

1 

XPCP 

194. G 

0 

693. A 

4.M 

730. 

198. 





7.R 




2 

XPCP . 

196. G 

0 

680. A 

7.M 

721. 

203. 





7.R 




3 

XPCP 

196. G 

0 

56. A 

lO.M 

215. 

206. 





l.R 




4 

XPQP 

480.G 

0 

384. A 

13. M 

630. 

493. 





4.R 




5 

XPCP 

254. G 

0 

78. A 

16. M 

282. 

270. 





l.R 




6 

XPCP 

613.G 

0 

445. A 

22. M 

782. 

635. 





5.R 




7 

XPCP 

613.G 

0 

266. A 

23.M 

692. 

636. 





3.R 




8 

XPCP 

(t)A 

0 

1310. A 

25000. G 

26371. 

25065. 





28. R 

32.M 



f?’ 

XPCP 

7340. G 

0 

1310. A 

32. M 

7494. 

7372. 

U) 


(t)A 


28. R 




8A 

XPCP 

11415. G 

0 

2350. A 

32.M 

11696. 

11^47. 





49. R 




9 

XPCP 

41100. G 

0 

1890. A 

42. M 

41188. 

41142. 





48. R 




9M 

XPCP 

40970.G 

0 

1890. A 

42. M 

41053. 

41012. 

(II) 




48. R 





(Continued) 




Table 3-10 (Continued) 


ANGULAR MGMENTOM 


XP-LV 


MAX OGUPCMOnS 


lOCUnG. I ATTITUDE I Hjj 


(HI) 


484% .G 0 22900.A 

305. R 

47860.G 0 22900.A 

305.R 

48450.G 0 23700.A 

310.R 

47910. G 0 23700. A 

310. R 

"o 0 23400. A 

(IV) 1230. R 

0 0 23400. A 

(IV) 1230. R 


(FT-IB-SBC) 


TOTALS 


MAX 


52. M 53790. 48524. 

52. M» 53241. 47914. 

I 


53. M| 54126 . 48505. 

53643. 47965. 


24630. 1328. 

(IV) (IV) 
24630. 1328. 


(t) a minor trace of angular momentum exists 

(I) jZ>Y = 25° to put the XP axis normal to the oibit plane 

(II) jjy <5° to put the XP axis noniial to the orbit plane 

(III) to put the XP axis along the local vertical 

(IV) This entry does NOT Include an Hjj ccmponent required for tnt leuvering 
about the X axis. 












TABl£ 3-11 

SUMMARY OF ANGUUR MCMENIVM 
DOE TO LCRG TERM BIAS TCTHQUES 
MAX TOTAL IS JCR PERIOD OF HIGH SOLAR ACTIVITY 


ANGULAR MOMENTUM (FT-IB-SBC/ORBIT) 


OONnG. ATTITODE 


MAX OCMPOKENrS 


HX «Y «Z 


200.AB 
155. R 


215.AB 
153. R 


8.AB 

2.R 


125. AB 
86. P, 
25. AB 
17. R 


140. AB 
100. R 


85. AB 
124. R 


(t)AB 220000.G 415. AB 

313. R 

(t)AB 0 415. AB 

313. R 

0 0 715. AB 

535. R 


20650. G 600. AB 

453. R 

0 600. AB 

453. R 


XITAIB 


380. 

181 

413. 

199. 

90. 

72. 

291. 

167. 

142. 

117. 

385. 

246. 

359. 

274. 

220200. 

220200. 

930. 

517. 

1450. 

738. 

20937. 

20915. 

1313. 

716. 


(Continued) 
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tabu: 3-11 (Continued) 


OCNFIG. 

ATTITODE 

10 

XPCP 

lOM 

(11) 

tl 

11 

tt 

IIM 

(II) 

n 

12 

XLV 

12M 

(III) 

XP-LV 


ANGULAR tOfETTrUM (FT-IB-SEC/GRBIT) 


MAX OCMPCNENTS 


% BZ 


61520.G 


11 " - 59490.G 7540. AB 340.M 61195, 

5279. R 

IIM " ■ 0 7540.AB 340.M 13159. 

(II) 5279. R 

12 jEv 124900.G 71680.G 350.M 

43500. A 
7450.AB 
5324. R 

12M XP-LV - 0 4'3500.A 350.M 

(III) 7450. AB 

I 1 5324. R 

(t) a minor trace of angular inxnentum exists 
(T) jSy = 25° to put the XP axis normal to the orbit plane 

(II) JJy P^^ ^ noimal to orbit plane 

(III) <10° to put the XP axis along the local vertical 

AB caused by half-wave rectification of aerodynamic cyclic torques 


124900.G 


TOTALS 




315.M '2719 . 5459. 


340.M 61195. 60067. 


340.M 13159. 5650. 


350.M 179158. 147191. 


5114. R 

0 7290.AB 315.M 

5114. R 


350.M 56624. 


5886. 
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The magnitudes listed under components are for maximum 
aerodynamic torques obtained for the rare, periodic 
occasions of high solar activity which results in max- 
imum atmospheric densities (at nighttime). Totals are 
listed for this maximum and also the norma] , where 
aerodynamic torques are at their normal levels. 

Angular momentvun buildups due to bias torques are 
presented in a similar manner in Table 3-11 . Here 

AB Represents aerodynamic bias torques cause by 

rectification of cyclic torques, where the vehicle 
is in daytime or nighttime for one-half of the orbit, 
the difference in air density causing an apparent 
half-wave rectification of aero torques. 

The resulting angular momentum listings given in Tables 
3-9 and 3-10 will be used to select actuator groupings 
in 3.3, while those listed in Table 3-11 will be used 
for determining a momentum desaturation system in 
Section 4. 

SHORT TERM DISTURBANCES 

Short term disturbances were defined mainly as forcing 
functions for the computer simulation of the SCB. 
Disturbances being considered are crew disturbances 
and Shuttle Orbiter docking. These disturbances can 
also be used in sizing the torqu: capability of the 
attitude actuator system. 



other disturbances which cannot and will not be considered 
at this time are: 

a. Moving parts of equipment (such as cameras, manipulator 
arms) 

b. Rotational wheel unbalances 

c. Attitude control system noise from sensors, quantiza- 
tion, etc. 

Locations of disturbances to be applied to the SCB simu- 
lations are listed in Table 3-12. All possible docking 
and/or crew disturbance force directions are noted on 
the table. Module locations and section numbers were 
defined in Section 2. 

. 1 Crew Disturbances 

Crew disturbance torques are basically caused by 
forces generated at a radius arm from the vehicle body 
center of mass. In vector form, 

L = R X F + M (3-24) 

where ' 

L = resultant torque disturbance 

R = radius arm from the center of mass to the point of 
application of the crew disturbance force 
F = crew disturbance force 

M = any disturbance moment generated as a couple. 
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TABLE 3-12 

DISTURBANCE LOCATIONS ON THE SCB 


MODULE 

BODY 

SECTION 

DISTURBANCES 

FORCE 

AXES 

CONFIGURATIONS 
FOR DOCKING 

CREW 

DOCKING 

FROM 

TO 

DM1 

1 

/ 

/ 

X 

1 

1 



/ 

/ 

y 

1 

5 



/ 

/ 

Z 

1 

3 

DM2 

6 

/ 

/ 

X 

2 

2 



/ 

/ 

y 

2 

4 



/ 

/ 

Z 

2 

12 

DM3 

7 

/ 

/ 

X 

8 

10 



/ 

/ 

y 

8 

12 


1 

/ 


z 

- 

- 

DM4 

10 

/ 

/ 

X 

8 

12 



/ 

/ 

y 

8 

12 



/ 

/ 

z 

8 

8 

HM3 

1 

/ 


X ! 

1 

(6)* 

(12)* 

(at the 


/ 


y 1 

(6)* 

(12)* 

end in the 

/ 


Z i 

(6)* 

(12)* 

-y direction) 



1 




♦crew disturbances only 




A time history of a "wall push off" crew disturbance 
force is shown in Figure 3-4; this was obtained from 
reference 3-12. The force can be restricted as being 
along the vehicle X, Y or Z axis, and the resulting 
scalar equations are 


Lx “ RyF^-R^FY+Mjj 


(3-25) 


Radius arm components (Rj^jRy.R^) were obtained frcii 
Section 2, and are listed in Tables 3-13 and 3-14. The 
conditions in Table 3-13 are for an ideally rigid SCB 
while Table 3-14 represents a flexible vehicle where 
modules are divided into bodies as in reference 3-13. 

Each body normally contains one or more modules. Table 
3-14 also lists the body number, SCB configuration number, 
body mass and the moment of inertia about its center of 
mass. Each body moment of inertia was computed about 
its i-th axis as follows: 

■'ll “ ” [ 
where 

m = the body mass 

= radius of gyration about the i-th axis 
r., r = module location along j-th and k-th axes with 

J K 

respect to the origin of SCB axis definition 
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TABLE 3-13 

PHYSICAL CHARACTERISTICS OF SCB VEHICLE 
DISTURBANCE LOCATIONS 

o RIGID BODY 
o i ■ X, y, z AXES 


DISTURBANCE 

LOCATION 

VEHICLE 
CENTER 
OF MASS 
(IN) 

RADIUS 

ARM 

(FT) 

WEIGHT 

(LBS) 

‘'ii 

MOMENT OF 
INERTIA 
(SLUGS) 

VEHICLE 

CONFIG. 

APPLICABLE 

MODULE 

POSITION 

(IN) 

DM1 

-377 

-76.8 

-25.0 

40,904 

2.6620E5 

1 


0 

0 

0 


5.3219E5 



0 

0 

0 


2.9685E5 


DM2 

578 

833.2 

-4.60 

588,637 

7.0445E7 

12 


0 

29.7 

-2.475 


2.0090E8 



0 

-48.1 

4.0083 


2.0986E8 


DM3 

1541 

ft 

75.65 

It 

ft 

12 


0 


-2.475 





0 

1 


4.0083 




DM4 

1541 

tf 

75.65 

Tf 


12 


0 

i 

-2.475 : 





250 


24.842 




^iM3 

-377 

ft 

-84.183 

H 

M 

12 


-514.2 


-45.325 

1 

i 




0 


4.0083 

1 

i 
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TABLE 3-14 

PHYSICAL CHARACTERISTICS OF DrSTURBANCE LOCATIONS 

FOR SCB BODIES 


o FLEXIBLE BODY 
o i « X, y, z AXES 


DISTURBANCE 

LOCATION 

BODY 
CENTER 
OF MASS 
(IN) 

RADIUS 

ARM 

(FT) 

gm^ 

WEIGHT 

(LBS) 

*^ii 

MOMFNT OF 
INERTIA 
INERTIA 
(SLUG-FT^) 

BODY 

CONFIG. 

APPLICABLE 

MODULE 

BODY 

POSITION 

(IN) 

FROM 

TO 

DM1 

1 

-377 

-88.742 

-24.011 

39,202 

3 . 1021E4 

1 

1 



0 

0 

0 


2.6466E5 





0 

0 

0 


2.6466E5 



DM2 

6 

578 

535.191 

3.567 

60,472 

3.3632E5 

5 

12 



0 

0 

0 


1.3981E5 





0 

0 

0 


4.2742E5 



DM3 

7 

1541 

1150.995 

32.5 

45,089 

3.5866E4 

8 

12 



0 

0 

0 


7.4545E5 





0 

0 

0 

1 

) 

•; "^73£5 



DM4 

10 

1541 

1541 

0 

82,080 

i42F' 

9 

12 

i 


0 

0 

0 


/ . ^ ' 1 





250 

1145.65 

-74.971 

1 


1 

I'. 

1 



HM3 

1 

-377 

-321.92 

i 

-4.59 

143,031 

1 , 3691E6 j 

6 

12 


1 

-514.2 

0 

-42 . 85 


1.5339E6 





0 

25.403 

-2.1x7 


2.1758E6 
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Cj,Ck =* locations of the center of masf of the tody 
along the j-th and k-th axes with respect to 
the origin of the SCB axis definition 


Also, 


n^l 

c « m 


(3-27) 


and m 





n 


(3-28) 


where n = module number 

N = total number of modules forming the body of 
interest 


3.2.2 Docking 

Docking disturbances to SCB will be as an initial 
condition o> angular rate, W(o). Equations for the 
solution of W(o) are generated in reference 3-14 utilizing 
Figure 3-5. A generalized scalar equation to be 
simulated for each case follows: 


**<'» 

W/o) 


R V -R V 

I 

J /m +R ^+R ^ 
xx' a y z 


R V -R V 
Z X X z 

J ,/m +R ^+R^ 
yy' a z x 


W^(o) 


R V -R V 

X y y x 

J /m +R ^+R ^ 
zz' ax y 


(3-29) 
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where 



4 


Vv. V„, V =■ linear relative velocity of the docking 
.. y z 

vehicle along the SCB X, Y, Z axes 

m mass of the docking vehicle 
a 

J ,J ,J * ^•'e moments of inertia of the space base 
XX yy zz 

about the vehicle or section of interest, 
about the center of mass of the vehicle 
'>r its section of interest 

Equation 3-29 assumes that the impact of docking is 
completely inelastic due to ,he operation of the docking 
latching mechanism. Also, the docking vehicle is 
assumed to be a point mass (or particle) once it is 
docked to the space base. A more complete derivation 
of thes ;uations would have included inertia dyadics 
for both bodies and the inertial translational and 
rotational velocities at the instant of impact. The 
resulting equations, however, would then be extremely 
complex. The docking condition would be assumed for a 
fully loaded Shuttle Orbiter at a relative velocity of 
0.5 ft/sec. The orbiter mass, m , will be assumed to be 

a 

7200 slugs (equivalent to a weight of approximately 
23^,000 lbs.). The magnitudes of Rj and can be 

found in Tablr., 3-13 or 3-14, dependent upon whether 
a rigid or flexible SCB is being simulated. Based upon 
conditions shown in Table 3-12, only one component, 
either Vj or Vj^, can exist during dockir.g. This can 
further simplify equation (3-29) by eliminating one 
of the terms in the numerator. 
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If translational docking disturbances on a body are of 
interest, the initial condition of a SCB linear velocity 
can be defined as: 

[ m 

Components of the Vg vector are scaled to V in a 
proportionate manner as shown in (3-30) 

where 

= mass of the space base 
T 

V = fv V V 1 
- L X y z J 

3.3 MOMENTUM EXCHANGE ACTUATORS 

On the basis of data presented in Sections 3.1.5 and 
3.1.6, the number and size of momentum exchange 
actuators can now be determined. An ideal momentum 
exchange actuator for a vehicle and mission of this type 
would be a double-gimbal control moment gyro (DG CMG). 
Each CMG unit would have a spherical momentum envelope, 
with modifications previously recommended to remove 
all gimbal stops. Assuming the same type as used for 
Skylab, the angular momentum for each DG CMG is 2300 
ft-lb-sec. The number of DG CMG units required for 
each configuration is given in Table 3-15. The largest 
number of DG CMG units in the Table is 24 for configura- 
tions 10 and 11. All 24 of these CMG units would also 




(3-30) 
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TABLE 3 -15 

NUMBER OF DGCMG UNITS REQUIRED FOR ANGULAR MOMENTUM ENVELOPE 
% CMC H 2300 FT-LB-SEC EACH 

• ASSL^S NO REDUNDANCY 

• SPHERICAL CMG ENVELOPE 




DGCMG UNITS 1 

CONFIG. 

ATTITUDE 

MAXIMUM 

NORMAL 

1 

XPOP 

2 * 

2* 

2 

ft 

2 * 

2 * 

3 

tt 

2 * 

2 * 

4 

t1 

2 * 

2 * 

5 


2 * 

2 * 

6 


2 * 

2 * 

7 

ft 

2 * 

2* 

8 

If 

12 

11 

8M 

XP-POP 

4 

4 

8A 

XPOP 

6 

5 

9 

ft 

18 

18 

9M 

XP-POP 

18 

18 

10 

XPOP 

24 

22 

lOM 

XP-POP 

24 

21 

11 

XPOP 

24 

22 

IIH 

XP-POP 

24 

21 

12 

1 XLV 

11# 

2*# 

12M 

XP-LV 

11# 

2*# 


♦MINIMUM REQUIRED FOR CONTROL OF VEHICLE 

# DOES NOT INCLUDE ANGULAR MOMENTUM REQUIREMENTS FOR ROLL AXIS 
MANEUVERING WHICH IS NECESSARY FOR POINTING SOLAR WINGS AND 
ARRAYS DIRECTLY AT THE SUN 
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be available for configuration 12, although only 11 
units would be required for the worst ease conditions. 

ROLL AXIS MANEUVERING 

The additional 13 units would be available for assistance 
in cycling the roll axis of the vehicle about the vertical 
axis for accurate pointing of solar panels throughout 
each orbit, if desirable. 


Reference to Table 3-9, however, indicates that a total 
of 59 or 34 equivalent CMG units would be required for 
vehicle roll about the local vertical with 3 angles of 
30 and 45 degrees, respectively. For 3=45 degrees, the 
additional angular momentuin about the X axis can be in 
the form of 10 reaction wheels (2300 ft-lb-sec), 10 
single gimbal (SG) CMG units or 10 additional DG CMG 
units. The reaction wheels (RW*s) would extend the 
momentum envelope of the system to the equivalent of 
34 CMG's in the + X direction, while 24 actual CMG's 
would still be available for control about the Y and 
Z axes. 

If 10 single gimbal CMG units were used, needless 
momentum would be available for the Y and Z axes (beyond 
the 24 DG CMG units which are already redundant). The 
10 SG CMG units would weigh more than the 10 RW's because 
of the gimbals and gimbal servos. Also, additional 
software would be required for the RW steering control 
laws. 

If 10 additional DG CMG units were added to the system, 
they would weigh more than the RW's or the SG CMG units 
because of the gimbal] ing and controls. The steering 



control law software, however, is easily integrated 
into the software existing (reference 3-15) for the 
original 24 DG CMC units. The full spherical momentum 
envelope of 34 DG CMG units would then be available 
and no additional steering control laws would be 
required. On that basis, the addition of the 10 DG CMG 
units appears to be the most attractive approach. 

The fixed roll angle approach (with the Y vehicle axis 
IN the Orbital Plane and with the Z axes NORMAL to the 
Orbital Plane) would be used for B angles less than 
45 degrees. For a Bangle just slightly less than 45 
degrees, only 71% of the solar power would be available 
at the solar panels. For smaller B angles, the efficiency 
would improve up to 100% when B = 0* For B angles 
larger than + 45 degrees, the roll about vertical 
maneuvering would be required and 100% of the solar 
panel electrical power would be available. 

Li the roll control transition is to take place at B=30 
degrees, then the minimum solar power efficiency would 
be improved to 86.6%; but 35 additional DG CMG units 
(or SG CMG's or RW's) would be required for configuration 
12. Figure 3-6 summarizes the trade off between solar 
power efficiency and the number of additional DG CMG 
units as a function of the B angle. The horizontal 
straight line is when the vehicle maneuvers in roll 
about the local vertical. The curved portion (actually 
a cosine B function) is for when the roll angle is fixed 
and the Z, or solar drive, axis is normal to the orbital 
plane. 
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FIGURE 3-6 


SOLAR POWER EFFICIENCY FOR X LOCAL VERTICAL ATTITUDE 
AS A FUNCTION OF B ANGLE AND NUMBER OF ADDITIONAL 
CMC UNITS REQUIRED 
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ALTERNATE NON-ROLLING SOLAR POINTING APPROACH 
An alternate approach would completely eliminate roll 
maneuvers around the local vertical for all B angles — 
with a slight increase in the average solar power 
efficiency. For B angles with magnitudes less than 36 
degrees, the same technique as previously described 
would be used: the vehicle roll angle would be set at 

90 degrees (relative to that shown in Figure 3-2) and 
the solar panels would be controlled to point towards 
the axis on the diagram. The efficiency would then 
be equal to cosB . 


For B angles larger than 36.1 degrees, the vehicle roll 
angle would be set to 0 degrees (Y axis normal to the 
orbit plane) and the solar panel drive would attempt to 
point towards the sun (maximize Yg«Xj). This can be accom- 
plished by using solar aspect sensors on the solar panels 
to set Xg«Xj=0. Ordinarily the solar power efficiency would 
vary in a cyclic manner as a function of orbit position 
angle a. The resulting instantaneous solar power 
efficiency would be 


- Yg.X 


I * t - 


(cosB sina)' 


(3-31) 


Also , 


tan6 - ^ (3-32) 

where 6 = solar panel position angle. 


An expression can also be generated for the average solar 

power efficiency as a function of the B angle using (3-31): 
tt/2 

n “ ^ J n da (3-33) 

o 


o 


! 



n 


AVERAGE 
EFFICIENCY 
(PER CENT) 


^ ANGLE MANEUVERS 



X LOCAL VERTICAL 


0 15 30 45 60 ?5 90 

6 ANGLE MAGNITUDE (DEGREES) 


A 

B: (fi 
C: 4) 


4) = 90" 


FIGURE 3-7 


(ZPOP), SOLAR DRIVE OPERATIONAL (POINTING TO Xq) 
(YPOP), SOLAR DRIVE FIXED (6 « 0 or 180°) 

(YPOP), SOLAR DRIVE OPERATIONAL (Xg • Xj * 0) 

SOLAR POWER AVERAGE EFFICIENCY vs. BETA ANGLE 
FOR THE ALTERNATE NON-ROLLING APPROACH 
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The integral portion of (3-33) (tabular data was obtained 
in reference 3-16) is an elliptic integral which results 
in curve C of Figure 3-7. Curve A is for the case when 
the solar panels point towards (when <(> = 90° or the 
Z axis of the vehicle is perpendicular to the orbit 
plane). A composite of curves A and C would then be 
utilized; they cross at 8=36.1° where the average 
solar power efficiency is approximately 80.8%. This is 
an apparent improvement over the method described in 
3.3.1 since no continuous roll axis maneuvering or 
additional CMC units are required. This should be the 
preferred approach. The CMC units indicated in Table 
3-15 would then apply. 


3-55 



REFERENCES 


3-1 Kaczynski, R. F. , Space Base Angular Momentum 
Requirements for Gravity Gradient Torques, 

MT 40,802, 24 October 1977 

3-2 Kaczynski, R. F. , Space Base Angular Momentum 

Requirements for Aerodynamic Torques, MT 40,803, 

17 November 1977 

3-3 Kaczynski, R. F., Space Base Angular Momentum 
Requirements for Magnetic Torques, MT 40,804, 

22 November 1977 

3-4 Kaczynski, R. F., Space Base Angular Momentum 
Requirements for Radiation Torques, MT 40,805, 
February 1, 1978 

3-5 NASA SP-8024, Spacecraft Gravitational Torques, 

NASA Space Vehicle Design Criteria (Guidance 
and Control), May 1969 

3-6 NASA SP-8058, Spacecraft Aerodynamic Torques, Anon., 
January, 1971 

3-7 NASA SP-8021, Model of the Earth's Atmosphere, 

Anon. , May, 1969 

3-8 NASA SP-8018, Spacecraft Magnetic Torques, Anon., 
March 1969. 

3-9 NASA SP-8017, Magnetic Fields - Earth and Extra- 
terrestial, Anon., March 1969 

3-10 NASA SP-8027, Spacecraft Radiation Torques, 

October 1969 

3-11 Kaczynski, R. F. , Solar Panel Position and Vehicle 
Roll Orientation During X Local Vertical Mode, 

MT 40,806, February 1978 


3-56 


3-12 Fritz, C. G.; Howell, J. T. ; Nicaise, P. D. and 
Parker, J. R. , NASA TM X-64972, A Miniaturized 
Pointing Mount for Spacelab Missions, November 
25, 1975 

3-13 Cornell, G. A., Space Base Mathematical Model, 
Bendix Research Laboratories Report, Program 
BRL-SB-1, December 14, 1977 

3-14 Kaczynski, R. F., Short Term Disturbances on the 

Attitude of the Space Construction Base, MT 40,809, 
May 8, 1978 

3-15 Kennel, H. F., Steering Law for Parallel Mounted 
Double-Gimballed Control Moment Gyros, NASA 
TM X-64930, February 1975. 

3-16 CRC Standard Mathematical Tables, 18th Edition, 

The Chemical Rubber Co., 1970 (see page 537). 


3-57 


I <1WtgR*«|j,i 'J 


SECTION 4 


4.0 MOMENTUM DESATURATION 

4.1 DESATURATION REQUIREMENTS 

Desaturation requirements are based upon the longterm bias 
torques determined in the preferred orientations ot jao; ‘^CB 
configuration. The requirements for momentum des.. curat io'. are 
listed in Table 4-1. The worst case average torque for one 
orbital period is also tabulated. Almost all of the disturbance 
torque is about the space construction base (SCB) Z axis for 
all twelve configurations. All of these are summarized in 
reference 4-1, 

Magnetic torquing, reaction control (RCS) and gravity gradient 
(GG) desaturation were evaluated for the XPOP orientations and 
also for the X local vertical alignment for configuration 12M. 
Characteristics of these approaches were independently deter- 
mined in reference 4-2. They will be reviewed in the following 
paragraphs and then compared at the end of this section. 

4.2 MAGNETIC DESATURATION 

A counteracting magnetic torque is generated on the vehicle 
by M, a torque coll system magnetic moment, interacting with 
the earth's magnetic field B according to the physical law 
given by: 


L * M X B (4-1) 

— — 

where each vector in (4-1) can be put in vehicle coordinates. 

2 

SI units for L^, M and B are N-m, A-m and tesla, respectively. 
The average torque vector per orbit should be equal and 
opposite to the average disturbance torque of the previous 
orbit, as measured from the momentum exchange control system. 


4-1 



SCB REQUIREMENTS FOR MOMEN*fUM DESATURATION 



! 

t 



CO 



M 


I: 

M 



> 



M 

CO 


E3 

h 


< 




o 


< 

HI 


O 

S3 

> 


00 







A 


5 



s 



HI 



X 

< 

I 


a 

w 

< 


OS 



p 



Si 




< 


M 



o: 


so 

< 

5 



a 

s 

» 

S 

o 

p 

a 

o 

a 

£-t 

H 

a: 




§ 

00 


< 

(0 

o 


CD 

S 


40 



CO 

s 


HI 

D 

» 




Q 


o 

O 

pi3 

s 

M 



a: 

O 


M 

S 


Oi 


< 


s 


p 

p 


< 

s 


H 

HI 

o 

H4 

X 

s 

OQ 

< 


cxi 


< 

O 

0 

0 

0 



pmmm 

o> 

H 


lO 

(0 


CO 

05 






P 

tH 

40 

w 

o 

O 

CO 

00 

H 

p 

CO 

N 



< 

CO 

CO 

H 

c>] 

CCI 



P 

05 

p 

O 


ca 

a 

o 

o 

o 

o 

o 

o 

o 

O 

H 


O 

o 

p 

S 

• 

• 

• 

• 

• 

« 

• 

• 

• 

• 

• 

• 

Of 

o 












H 

8 

CO 

< ^ 
sa, 

!2C 


























» Eh 

a hi 

a 

o 

0> 

m 

CO 

H 

0> 

CO 






o: 

p 

N 

C3l 

40 

H 

40 

N 

CO 


05 




M 

a 

CD 


H 

40 

CQ 

p 

to 

to 

CO 


05 

05 

> 

HI 

O 

o 

o 

o 

o 

o 

o 

ri 

05 

05 

CO 

05 

< 

X 

• 

• 

• 

• 

• 

• 

9 

* 

« 

• 

• 

• 


< 










05 

05 

05 


a 














< 

H 

CJ) 




to 

rf 


to 

05 

o 

CD 


a 

00 

0> 


p 

H 


0* 

H 

H 

40 

40 

00 


cc 

n 

H 


rW 

tH 

05 

05 

40 

0* 


CO 

00 


o 










40 

40 

40 

HS 

? 3 

2 K 

go 

a 
























— 
















o o 














X M 














50 














b: I 














< ca 

a 













•j ^ 


o 

CO 

o 

H 

05 

40 

05 

O 

CO 

05 

05 


D 1 

a' 

00 


0) 

cr 


00 

lO 

'0 

H 

H 

40 

05 

O H 

HI 

CO 



<N 

iH 

CO 

CO 

05 

CO 

N 

ri 

CO 

2 

X 









H 

05 

CO 

CD 

< W 

< 










H 

H 

40 


a 













" 1 














2 

Q 








Dk 





O 

M 

Pi 

H 

0< 







2 




> 


NH i 

r: 

i 

z 

r 

z 

• 

z 

r 



r 

- 

1 

& 

P 













X 

o 












■MHWMI 


M 









a 

a 

a 

a 

a 

b 


H 


CO 


40 

to 

0- 

00 

o> 

o 

rH 

05 

§ 

• 










iH 

tH 

f-l 

o 

o 













o 

- 1 

ss; 

-J 
















4-2 






^ j m m m »» 


A basic orientation of four magnetic torques relative to vehicle 
axes is given in Figure 4-1. If X is set at 35.26 degrees, each 
component of the magnetic momemt vector m « mg m^ 

has direction cosines of equal magnitude with the X, Y and Z axes. 
A block diagram for the typical implementation of a magnetic 
torquing system is shown in Figure 4-2. Both the minimum energy 
(ME) and cross product control laws are described. The cross 
product law is easier to compiu-hend and implement, but the ME law 
is expected to use less electrical energy and could be a candidate 
for the SBC application. 


An explanation of some of the symbols follows: 


B * predicted magnetic field of the earth at the instantaneous 

location of the SCB in orbit, in vehic*;? coordinates (tesla) 

B * actual B (tesla) 

* measured B (tesla) 

—fn 

— AVG * average bias torque (previous or presently predicted orbit) 
in vehicle coordinates (f.-m) 

L * generated magnetic torque in vehicle coordinates (N-m) 

m,m„ “ actual and commanded magnetic moments, respectively, in 

° 2 
coil coordinates (A-m ) 

M,M = actual and commanded magnetic moments, respectively, in 

— 2 
vehicle coordinates (A-m ) 

T = time period of ME control law solution (sec) 

AH »• angular momentum vector change (or commanded dump) over the 
period T (N-m-sec) 
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BLOCK DIAGRAM OF MAGNETIC TORQUING 
FIGURE 4-2 





i 


' ! 


I 

' J 


X 

X 


Note that the diagram assiimes m * The first two blocks • i | 

on the diagram must be implemented in the computer, while . . \ 

the second two blocks are inherent to the individual orient- ^ ^ 

at Ions and the operation of the torquer colls. 1 

l ; 

Table 4-2 is a summary of the magnetic torqulng system char- i 

i 

acterlstlcs. An average earth's magnetic field flux magnitude t ; 

—4 

of 0.1 X 10 tesla is assumed to be available on either the ' ‘ ; 

vehicle X or Y axis. This is a fairly conservative estimate 

I 

for all orbits. ■ • ' 


From the block diagram we see that 

Lm 2 “ 0.577 + mg + m^ + m^) B - 2.3094mB (4-2) 

gives the Z-axls magnetic torque if the directions of the mag- 
netic moments are properly obtained from the control law. 

Utilizing the torque requirements from Table 4-1, the required 
magnetic moment per coll (m) was calculated. It was also re- 
calculated for the case when one coil (as shown in Figure 4-1) 
has malfunctioned. A magnetic torquer size of 3300 A-m^ was 
selected for configurations 1 through 7. 

Each coll is estimated to be 84 inches long and 3.1 inches in 
diameter. The weight of each magnetic torquer will be 90 pounds 
and will require 22.5 watts for a maximum magnetic moment. 

The magnetic torquing system was tripled for configurations 
8M and 9M, and multiplied by another tenfold for configurations 
lOM and IIM. If magnetic torquing were to be used for con- 
figuration 12M (X local vertical) without any additional gravity 
gradient or RCS assistance, the system would have to be multi- 
plied 130-fold when compared to the original system for the earlier 
configurations, the reason for this being the large solar panel 
surfaces being located so far from the center of mass. 
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♦Thruster at - X end only, C.M. is at -76.8 in. for Configuration 



RCS Desaturation 

An RCS desaturation system would be operated on demand, l.e. 
whenever the momentum envelope of the CMG system Is approached 
or whenever the opportunity exists. RCS thruster locations 
are assumed to be available at the largest radius arms along 
the X-axls and to operate as a couple (for all except con- 
figuration 1). A specific Impulse Igp was assumed to be 300 
seconds. The time period At for each configuration was as- 
sumed to be four months each, also the worst case average 
torque was assumed to exist over the complete A t period of 
operation. The propellant weights per configuration and the 
total accumulation were then computed. All of these char- 
acteristics were tabulated In Table 4-3. Although the propel- 
lant weights do not Include tankage, valving, piping and 
thrusters, they appear to be very competitive with magnetic 
torque desaturation. The weight numbers are quite low con- 
sidering that the fuel Is consumed over a four year period. 

Gravity Gradient Desaturation 

Whenever the SCB Is not In a precise pointing period during 
Its mission. It can be tilted slightly from Its nominal 
orientation to generate a gravity gradient torque which is 
equal and opposite to the cumulative bias effects of aero- 
dynamic, radiation, magnetic and also gravity gradient dis- 
turbances. For example, in the XPOP mode a slight tilt about 
the Y or Z vehicle axes will generate desired gravity gradient 
bias torques about the Y and Z axes, respectively. Torques 
may also be generated about the X axes through a combination 
of Y and Z axes tilts; this effect, however, will be much 
weaker because of the actual inertia properties of the SCB. 

The approximate (small angle) gravity gradient bias torque 
equations, which assume a circular orbit and neglect product 
of inertia terms, follows: 
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SIN 0 SIN 'i' 

A 


Ly . Ky SIN 0 


(4.3) 


L = K SIN 
z z 


where = 3W^ (I^^ - lyy) /2 


^ (I,^ - I^,) /2 


<'yy - 'xx> /2 


I„„, I„„, I^_ = moments of inertia about vehicle axes 
xa yy zz 

“ orbital rate (rad/sec) 

0, 'F = small tilt angles about the vehicles Y and Z axes, 
respectively. 

A similar relationship also exists for the X local vertical 
attitude during configuration 12. 


Ordinarily, the tilt angles will be constant over the whole 
orbit and will be updated every orbit or so. But when an 
inertial orientation is required for celestial viewing, there 
are some problems in using gravity gradient desaturation; 

(a) Large angular momentum build-ups may require large tilt 
angles during occulted portions of the orbit. 
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(b) An additional momentum capability would be required 
for slewing through large angle maneuvers. 


(c) Most of the disturbance bias torques will be due to 
gravity gradient and aerodynamic effects; if the aero- 
dynamic bias torques change appreciably for large 
tilt angles, it is possible that an angular momentum 
buildup may not be bounded by simple gravity gradient 
control. 

(d) More complex software. 

For the reasons stated above, it is recommended that another 
desaturation scheme (magnetic) be available for SCB configu- 
rations requiring orientations other than XPOP or X local 
vertical. Gravity gradient desaturation, however, does have 
some advantages for these preferred orientations: 

(a) It is virtually weightless and requires no power (no 
additional hardware). 

(b) Software for continuous tilting is fairly simple 

(c) No contaminants are released in the immediate vicinity 
of the vehicle. 
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4 . 5 Sununary 

The three desaturation schemes which were discussed may be 
compared by listing some of the disadvantages of each as 
listed in Table 4-4. The magnetic torquing system appears 
to be the most flexible approach. For the last few con- 
figurations (lOM, IIM, 12M), however, the magnetic system 
gets extremely heavy in weight. For configurations lOM and 
IIM RCS thrusters maybe mounted on the CCM and LMl modules. 

The RCS will then be the primary desaturation approach for 
1011 and 1111, with some assistance from the magnetic torquing 
system. It would actually be recommended that the RCS system 
be Installed prior to configuration 8M or 9M and be available 
as a backup or for an emergency desaturation condition. 

No additional desaturation equipment would be required for 
configuration 1211, since a "fixed" gravity gradient tilt 
angle from the local vertical can be used to counteract a 
long term angular momentum buildup. One additional compli- 
cation with 1211, however, is the continual cyclic rolling about 
the local vertical for 3 angles greater than 30 or 40 degrees. 
The CMC angular momentum buildup should be separable from the 
cyclic momentum for controlling the roll angle about the local 
vertical: this buildup would be an input signal to the gravity 
gradient moment management control law. 

A summary of recommendations for SCB momentum management is 
shown in Table 4-5. A magnetic torquing system can be used 
for configurations 1 through 9M, with the RCS system used as 
a backup for 8M and 9M. The logistics of propellant fuel re- 
plenishment and thruster relocation for an RCS desaturation 
system ^.s thus circumvented for the first seven configurations. 
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Rather than a tenfold increase of the magnetic torquers for 
configurations lOM and IIM, the RCS system can become the 
primary desaturation system at that point. Progressing to 
configuration 12M would ordinarily require a larger increase 
in propellant weight. But fortunately configuration 12M is 
in an X local vertical orientation. A fixed small tilt angle 
from local vertical can be utilized such that a fixed gravity 
gradient torque cancels the large Z axis aerodynamic torque. 

The twelve magnetic torquers and the existing RCS thrusters 
can be used for backup and for short term assistance in de- 
saturation. 
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SECTION 5 


5.0 OONIBX SYSTEM APEfOACH 

5.1 INIKmTTION 

This section Is part of • study of control system concepts for the 
Space Construction Base. This effort is directed to investigate control 
during buildup, i.e., vMle the Space Construction Base is being assembled 
in low earth orbit. Study emphasis is on coupling control and control 
coordination since these are identified as problem areas peculiar to this 
application. This section presents a clarification of anticipated design 
problems and control oonc^ts. 


Within the next 10 years, earth aatellites are expected to include 
a new breed of large structures constructed in space. Each one may serve 
a wide range of functions from scientific experiments to nanufacturing 
in the zero-g environment. Others may have special purposes such as 
collecting solar energy and either retransmitting it to earth via micro- 
wave or servicing other space systems. Hany of these satellites will be 
constructed from sections transported by the Space Shuttle. 


These large structures introduce new challenges to control system 
design. The changing physical structure during construction affects the 
optimum distribution of sensors, the control system structure, and con- 
trol parameters. The large size and requirements for minimizing weight 
can bring structural modes within the control systeai bandwidth - result- 
ing in stability problems and coupling effects. 
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The Spece Conetructlon Base provides a conceptual baseline for 
various ayotea studies related to these large space structures. It is of 
"Intersedlate site" %rlth outer dimensions of 320 by 420 tt. Ic Is 
conceived to be constructed in a sequence of twelve configurations over 
a period of four years or longer. 


The objective of the study described In this s&' i ionis to Investi- 
gate Space Construction Base control problems throughov buildup and 
during any dynamic testing while in low earth orbit (XX.0) . The control 
system must maintain overall stabillty» include vernier control for at 
least one element or module, and provide mmneuvering capability. 


The emi^asle in this work is on those control problems which are 
peculiar to vehicles such as the Space Construction Base. This Includes 
appendage and module stabilising control in which the coupling between 

modules and appendages, shape of Individual modules or appendages, and 
the damping of structural modes are controlled. This also Includes the 
vehicle control coordinator idiich must vary control parameters and 
select sensors according to configuration and performance changes. Spe- 
cial problems also include model uncertainty and reliability. 
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Sl2 0MK)L PROBUEUS and APRCACH 


Eventually, it is desired to define a control system concept which 
will stabilize the Space Construction Base attitude, in whole and in part, 
and provide vernier, maneuvering, and momentum management capability. 

Such a system must be designed in light of model uncertainty, adapt to 
conf: guration changes, function in the presence of strong intermodule 
interactions and low structural damping, and have a high degree of system 
integrity with respect to possible component failures. 

The Space Construction Base concept introduces a special set of 
control problems. The control system must accommodate itself to changes 
in the structure during buildup. It must be relatively insensitive to 
uncertainties in the structural dynamics since they will not be known 
precisely. The structure must be regarded as being flexible due to its 
large size and due to the need to minimize its weight. It must be recog- 
nized that the Space Construction Base will be used for a wide variety 
of purposes. 

The purpose of this section is to outline a control system concept 
for the Space Construction Base. Emphasis is given to those elements 
which are peculiar to this application. The categories of control and 
control requlreir.ents are discussed in Section 5.3. This is followed 
in Section 5.4 by a description of the overall control system concept. 
Important problems special to the Space Construction Base are outlined 
in Section 5.5. Sections b.6 and 5.7 outline approaches for two control 
functions of importance for large spacecraft. 

5.3 CONTROL FUNCTIONS 

As a convenience, the control of the Space Construction Base may 
be divided into several functional areas: 
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A. Appendage and module stabilizing control includes coupling 
control, shape control, and artificial damping of structural 
modes. . Coupling control ranges from complete decoupling to 
artificial stiffening in which one module is slaved to 
another. These functions must be accomplished while the 
Space Shuttle is docking, station-keeping or maneuvering. This 
functional area is somewhat analogous to th. stability augmen- 
tation function of an aircraft. 

B. Attitude control holds a given overall orientation of the 
spacecraft. This would be referenced to an inertial frame 
or local vertical. Typical accuracy requirements would be 
on the order of 0.5 deg. 

C. Vernier control holds an instrument, element, or module to a 
precise attitude. Typical accuracy requirements are 1 to 

10 arc-s. 

D. Maneuvering control changes the overall orientation of the 
Space Construction Base such as required for docking or for 
minimizing gravity gradient torque after a module has been added. 

E. Momentum desaturation is required with moment exchange devices 
such as reaction wheels (RMs) and control moment gyros (CMGs). 
Thus, secondary actuators are used along with CMGs and RWs. 

The functional structure of the control system design reflects 
the philosophy and complexity of the control system. It must express 
the delegation of computer authority, a definition of control loops, 
and guidelines for sensor and actuator distribution. Further definition 
of the functional design would include control loop parameters such as 
gain factors, filter constants and sensor and actuator dynamic require- 
ments (bandwidth, position and rate limits, etc.) 

The delegation of computer authority is tied to the use of central- 
ized or decentralized control. The viewpoints of hierarchical and multi- 
level control may be adapted to this application. Decentralization has 



bfe'one more attractive for digital control systems recently vlth the 
growth of microcomputer technology. Intervention from ground control 
should be sn Integral part of the concept. 

Tlie extent of coupling among control loops Is an Important part 
of the functional design. Artificial stiffening Is an attempt to force 
the space vehicle to act as a single rigid body* l.e., control action 
substitutes for a very h«isvy rigid structure. Complete decoupling is an 
attempt to cancel Interactions between the control loops. Partial de- 
coupling limits the interactions to assume sufficient stability margins. 
Artificial stiffening between modules would be desirable during maneu- 
vering so as to rotate the v^icle as a unit. Decoupling control, how- 
ever, may be more desirable while holding a given attitude so as to 
minimize the coupling of disturbances. Notch filters serve to decrease 
the interactions between control loops and the structural modes but make 
the design sensitive to model uncertainties. 

Guidelines should be established for sensor and actuator distri- 
bution. These will establish the relationship of location selection 
to mode shapes, function of the given nodule, and center of gravity. 

The number of sensors and actuators on a given module also depends on 
these considerations. The possibility of Including more sensor locations 
than actuators Should also be considered since mSny sensors are of lighter 
welgd^t and lower cost than the actuators. Also, optimum sensor loca- 
tions may then be selected according to changes in structural im>des. 

The mechanization concept would Include many hardware selections. 

For exanq>le, there would be some consideration of using digital or analog 
computers at each level of control. Sensor choices may Include different 
classes of gyros and accelerometers, as well as strain gages. Actuators 
may Include CMGs, torquers /motors, and/or reaction Jets. Possibly the 
earth's magnetic field combined vlth local electromagnets may be used 
for momentum management. 
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5.4 SYSTEM CONCEPT 


Figure illustrates a functional concept for the Space Ccnatruc- 
tion Base control system. It is intended to fulfill all of the categorle> 
of control outlined in Section 5.3. Each block vlll be changed as the 
Space Construction Base is assenfcled. 

The Vehicle Control Coordinator %n>uld change control parameters 
as the configuration changes. The selected parameters may be predeter- 
mined nominal values corresponding to each configuration, values based 
on sensed behavior of the system, and/or values remotely commanded by 
a ground station. It would also be capable of selecting or mlghtlng 
sensor outputs based on estimated modal shapes. The coordinator may 
include an identifier for estimating the dynamics of the vehicle. The 
Identifier may simulate an adjustable model which is updated based on 
observations of the actuator commands and sensor outputs. Sone of the 
coordinator functions may be performed in a ground computer. The 
vehicle control coordinator is discussed further in Section 5.7. 

The Vehicle Attitude Reference Unit may be located in the Habita- 
bility Module or the Subsystems Module since these are Included in all 
configurations. It would include gyros, a fixed reference system such 
as a star tracker, and a computer Cor computer segment) for strapdotm 
guidance computations. 

The Vehicle Manueverlng and Attitude Control unit seeks to maintain 
a desired overall velilcle attitude or !^lew rate. Primary control will 
likely be achieved through use of the double ginhal control moment gyros 
(DGCMCs). 

The priamry actuators for momentum management are assumed to be elec- 
tromagnetic. Thrusters should provide backup. A magnetometer is likely 
needed to determine ambient magnetic field so that electromagnetic re- 
quirements aay be computed to achieve a given moment. 

Vernier controls are separate functions. A separate set of sensors, 
control logic and CMGs is required. These would be localzed to vidiicle 
elements such as the research pallet. 
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figure 5-1 - Space Base Control Concept 









The control concepts described in the above discussion are slinilar 
to those foimd on many existing spacecraft. The coordinator functions 
are probably more important and complicated than In many other space 
vehicles. However, the newest feature is the Appendage and Module 
Stabilizing Control (or just stabilizing control). This feature is 
discussed in Section 5.6 with supporting simulation and analysis results 
given in Section 6. 

The Force and Torque Sensing and Processing Units provide primary 
inputs to the stabilizing controls. The applicable sensors are strain 
gages from which coupling torques and relative angles are computed. 

The coupling torques are at the interconnections between modules and 
between modules and appendages such as the solar panels, external tanks, 
and antenna support structure. The strain gages will also be used to 
provide vibration modal data for the Vehicle Control Coordinator and 
Ground Control. The nuai>er of sensor locations may be limited by such 
considerations as computer capacity and calibration requirements. Al- 
ternate sensors such as optical devices may also be considered for stabi- 
lizing control. 

The primary actuators for the coupling control function of the 
stabilizing controls are appendage to..~ue motors. Shape control and 
artificial damping may also be accomplished with torque motors. For 
larger appendages such as the large solar panels, CMGs may be used for 
actuation providing there is another type of actuator for momentum 
desaturation and for balancing large constant torques. 

5.5 SPECIAL PROBLEM AREAS 

Because of its large size and its buildup in orbit, the Space 
Construction Base presents special control problems. These must be 
considered in any trade studies and in design development. 

One of the primary considerations is model uncertainty. It 
is not practical to conduct meaningful preflight dynamic tests of the 
entire structure. In fact, the final configuration may be unknown 
when the first sections are placed in orbit. Modes of multiple-connected 
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Table 5-1 - Control System Design Trades and Considerations 


Functional 

Elements 

Design Feature 
Trades 

Considerations 

Vdilcle Control 
Coordinator 

Functional assignments 
to ground control or 
vehicle 

which control parameters 
should be adjustable 
which sensors should 
be selectable 

Space computer requirements 
Response time 
Flexibility in adapting 
to configuration and 
operational changes 
Reliability 

Design sensitivity to model 
uncertainties 
Capability to Influence 
dynamic performance 
Cost 

Actuators 

Application of rotary 
actuators versus 
CMGs to appendage 
and module stabili- 
zing 

Actuators for desatu- 
ration of CMGs 
(passive devices 
versus mass expul- 
sion thrusters) 

Capability to Improve dynamic 
behavior 

Availability of appropriate 
actuators 

Power requirements 
Size and weight 
Reliability 
Cost 

Desaturation requirements 
Resupply requirements 

Sensors 

Strain gages versus 
optical devices for 
relative angle data 

Size and weight 
Accuracy 

Sensitivity to model uncer- 
tainties (including modal 
shapes) 

Dynamic response 

Linearity 

Range 

Reliability 

Cost 

Calibration requirements 
Multiple usage capability 
Power distribution requirements 

Appendage and 
Module Stabili- 
zing Controls 

1 

1 

where to insert 
coupling controls 
where to Insert 
shape controls 
where to insert 
artificial damping 

Capability to lnq>rove dynamic 
behavior 
Reliability 
Cost 

Actuator requirements 
Power requirements 
Sensor requirements 
Sensitivity to SK>del 
uncertainties 
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large etructures cannot be accurately predicted analytically. Alao, 
the analytical nodela uaed in control aystem design are necessarily 
aliq>lified. The control systen must be designed to satisfy performance 
requirements in spite of these uncertainties. 

Another consideration is that the ci^ntrol configuration will change 
as nodules are added. Control of a given module may require sensor data 
from other modules leading to greater complexity as others are added. 
Actuator requirements on a given module m'.y depend on the number and 
size of other modules, as well as local weight and wel^t distribution. 

The control of a g. an module must adapt to disturbances propagated 
from other modules. As modules are added, new structural modes are 
added and others shifted. Fixed sensor and actuator locations on a 
given module may not be optimum for all configurations. 

Even if the configuration of the Space Construction Ease were 
constant, there would still be special control problems. Strong inter- 
action among the modules may lead to a need for somu degree of decoupling. 
Since a complex control system can be highly susceptible to component 
failures, the system integrity must allow satisfactory performance of 
many functions in spite of some loop failures. Control system computa- 
tional requirements must also be considered. There are restrictions 
on locating actuators. For example, the small solar panels will un- 
likely support the mass of control moment gyros (CMGs) . 

Table 5—1 is a summary of aone control system design considerations 
for the Space Construction Ease. 

5.6 APPENDAGE AND MODULE STABILIZING CONTROLS 

The controls under discussion in this section include coupling 
controls, shape control and artificial damping. They serve to stabi- 
lize the Space Construction Base by controlling either module inter- 
actions or individual modules directly. 
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Potential actuators include rotary torquera, control aioment gyros 
(CMGs) • and notor driven cable winches. The rotary torquers would be 
placed at the Joints between modules and used primarily for coupling 
control. Some rotary actuators will be free to rotate to large angles - 
e.g. , a solar wing must rotate to maintain pointing toward the sun. 

CMGs may be used for any of the above functions but must be supplemented 
with other actuators for momentum desaturation. Cable winches may be 
used in coupling or shape control if the dynamic response requirements 
are minimal. For example, such a situation would exist if there is a 
large constant torque required to maintain a desired shape (to overcome 
a spring restoring torque) . 

Coupling control may accomplish either decoupling or stiffening at the 
points between modules and between module and appendages. Both charac- 
teristics may be obtained through control mode switching. Decoupling 
control would be applied %rhen it is desired to minimize the transfer of 
disturbances. On the other hand, artificial stiffening may be applied 
during maneuvers so that the structure will rotate as a whole. 

Figure 5-2 Illustrates how a rotary actuator may be interfaced to 
an appendage hinge joint with the main body. Double Gimbaled Control 
Momemt Gyros (DGCMGs) are assumed to be controlling the main body. As 
shown, the rotary actuator will apply a torque to both the main body and 
appendage. The reaction torque on the main body may be cancelled with 
a command to the DGCMGs. Alternately, the attitude control on the main 
body could be allowed to make the necessary corrections. 


Decoupling control is shown in Figure S-3. The main body and 
appendage each has its own control loop (via and Kg^) ^or stabiliza- 
tion. The decoupling serves to cancel or partially cancel the coupling 
terms shown in Figure 5-2. The extent of decoupling is governed by the 
gain factors, and (0 S ^^ci — 
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Figure 5-2 - Concept for Main Body and Appendage Actuation 
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Figure 5-3 — Concept for Decoupling Control 
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Artificial stiffening control Is Illustrated in Figure 5-4 In 
this case, the appendage attitude connand is set to the aain body angle. 
There is no decoupling and the control parameters will generally hav? 
different values than with decoupling co. .rol. The objective is to 
slave the appendage angle to that of the main body. This would be 
desirable during any maneuver. 

Both the decoupling control and artificial stiffening incl'jde an 
appendage control loop with the parariters &Qd K’^. Th:'s configura- 
tion has been found convenient frr the study of decoupling and artificial 
stiffening controls, but it is not necessarily the best for the damping 
of structural modes. Therefore, alternate configurations such as that 
of Figure 5-5 should be Investigated. The concept shown Includes a band- 
pass filter to Isolate the mode being damped and has a higher rain to 
be responsive to the controlled mode. 

The use of coupling control will certainly be selective. That is, 
for a given appendage and axis of rotation, the coupling may lert 
uncontrolled. For example, the couplings from the solar wings to the 
main body about the Z-axis (wing torsion) will probably be uncontrolled 
" 0 in Figure 5-3 . This is because a motion of a solar wing about 
this axis will have little effect on Che main body due to the large 
difference in moments of inertia. However, coupling from the main body 
to each solar ying about the Z-axls will probably be controlled (K ^2 ^ 
since any disturbance on the main tody will otherwise propagate to the 
solar wing. 

5.7 VEHICLE CONTROL COORDINATOR 

Due to its large size and the buildup process, the Space Construc- 
tion Ease control system will require some means to adjust control parame- 
ters and select sensors. These functions are achieved through the Vehicle 
Control Coordinator with some interaction with Ground Control. 
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TWo concepts for asking these sdjustaents ere discussed briefly 
In this section. These represent two extreaes In the level of coaplexlty 
of the spscebome control system. The ectual configuration asy include 
features of each concept but vlll probably be closer to the sis^ler 
approach of Figure 5-6. 

Figure 5-6 Illustrates a concept with a minimum level of on-board 
capability. All of the central parameter and sensor selections are 
determined on the ground through the mission support function. The ground 
simulation is updated based on performance data received from the space- 
craft and is used for determining desirable control parameters and sensor 
locations. The central ccitrol tmit functions primarily as memory so as 
to minimize radio communication requirements. There is some capability 
for manual intervention by space personnel. Those variables which have 
errors characterized as stochastic are estimated with a Kalman Filter or 
Extended Kalman Filter. Ibose which are characterized as deterministic 
are computed with an Observer. 

Figure 5-7 shows a concept with extensive on-board capability. 

The control parameter and sensor selections are determined primarily 
from on-board computations in the Vehicle Control Coordinator. The 
groiind simulation serves to Interpret performance data transmitted from 
the spacecraft. Ground Control provides overall instructions such as ma- 
neuvering and docking commands. Ground Control may also provide a backup 
to space computations. The Control Configuration Controller* Identifier, 
on-board dynamic model simulation, and estimatcr-observer provide adap- 
tive control capability. The identifier coi>q>ares the system performance 
with the dynamic model and adjusts the model accordingly. The compari- 
sons may be made using test signal inputs or be based on the response 
to normal disturbances. 

5.8 pistri&ution of sensors and actuators 

A basic part of the control system design will be the selection 
of actuator and sensor types and locations. This will be based partly 
on the determination of torque requirements for countering disturbances 
such as gravity gradient and aerodynamic drag. 
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Tabic 5-2 liata aoaa of tha trades associated vlth soxsor and 
actuator distribution. The extent to which performance such as stability 
and pointing accuracies can be improved over uncontrolled behavior is of 
fundamental importance. Certain locations such as the small solar wings 
are not feasible for relatively large and heavy actuators such as CMGs 
and the accompanying momentum desaturation actuators. Structural modal 
shapes must be considered so as to avoid strong actuator Interactions 
and large sensor errors. 

Table 5-2 - Considerations for Sensor and Actuator Distribution 


Item Under 
Consideration 

Alternatives 

Considerations 

1. Actuator Locations 

Hodules on main 
body 

On appendages 
At appendage hinge 
Joints 

At module intercon- 
nections 

Performance require- 
ments 

Effectiveness in ob- 
taining desired 
moment or force 
Size and weight 
Structure modal shapes 
Safety 

Access for maintenance 
Power distribution 
requirements 
Momentum desaturation 
requirements 

2. Sensor Locations 

Modules on main 
body 

On appendages 

Size and weight 
Sensitivity to desired 
measurements 
Structural menial shapes 
Power distribution 
requirements 
Access for maintenance 











5.9 MULTIVARIABLE DESIGN ANALYSIS 

Although the priaary tool for control systea design Is coaputer 
slaulstlon, there ere soae analytical approaches that are useful in the 
design derivation. Most of the single-loop aethods are well known and 
do not require review here. 

Multivariable or aultiple-loop control » however, is an area of 
continuous research - warranting acre discussion. Mjltivariable 
analytical approaches include: 

A. Multilevel Control - This is a viewpoint of systeas as con- 
sisting of layers of decision levels. An overall systea is 
divided (or decomposed) Into a hierarchy of goal-seeking 
subsystems or decision problems. 

B. Linear Quadratic Regulator (LQR) - This approach, also known 
as Linear Optimal control, views the systea as a %ihole with 
the performance specified in terms of a single index. 

C. Multivariable Nyquist Array (MNA) - This is an extension of 
frequency response aethods to multivariable control systeas 
and includes the concept of partial decoupling. 

D. Characteristic Loci - This is another extension of frequency 
response aethods to multivariable control systeas. It is 
analogous to the root locus technique and does not require 
partial decoupling as a separate step. 

E. Other Algebraic Techniques - Pole/zero placement and aatrlx 
transfer functions may be analyzed from a variety of alternate 
viewpoints which have appeared recently in the literature. 

Figure 5-8 shows a typical application of linear optlul control. 
The design procedure is conceptually straightforward using atvallable 
computer prograas. The gain aatrlx is often taken as a steady-state 
constant aatrlx, but it aust be changed as modules are added to space 
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Figure 5-9 ~ Digital Control System 







base. An estinatlon Is required if the measureDent noise and disturbances 
are SK>deled as randou processes. An observer is applicable to detemlnls- 
tic systeas which have some states that cannot be measured directly. 

In practice, the design procedure is iterative. First, the choice 
of an approximate linear model must be made. Second, a performance index 
with desired weighting factors on each state and control variable must be 
selected. For stochastic representations, the statistics must be deter- 
mined, although they are generally unkno%m. After a design has been 
computed, the result should be evaluated on a more exact, possible non- 
linear simulation. Since that result is unlikely to satisfy the designer, 
the procedure must be repeated with different performance weighting or 
a different linear model. 

Figure 5-9 shows a digital control representation of the overall 
space construction base control system - linearized about some reference 
condition. The control parameters for the "digital controller" nay be 
derived using the LQR method. The first step would be to derive a con- 
tinuous model for ecch configuration of the control actuators (e.g. , 

CMC servo loops), the overall ^pace base dynamics, and the sensors. 

Each model must be linearized and simplified for representation of per- 
turbations about some reference point. Each model will have a separate 
set of control parameters, l.e., an LQR design. 

Each continuous model Bust be converted to a discrete model. This 

may be accomplished once a sample rate has been selected. The state 

transition matrix is obtained by summing matrices obtained from a trun- 

AT 

cated exponential series expansion (of e where A ■ continuous plant 
matrix and T ■ sample period). For disturbances and noise, the covari- 
ance smtrlces are converted by dividing by the sample period. 

The objective of the control system design is to minimize a per- 
formance criterion: 



i%(i) F(i) X (1) + u’^d - 1) G(1 - 1) u(l - 1)> 
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where 


X(l) • column state vector at i'tli sample time (n components) 

u(i - 1) •* column control vector at (1 - l)st sample time 
(r components) 

F(l) ■ n X n symmetric positive senl-deflnlte matrix (plant 
state variable weightings) 

G(i) r X r symmetric positive semi-definite matrix (cortrol 
variable weightings) 

E(Y) ■ expected value of Y 


In order to solve the optimal control problem, it is first necessary to 
select the performance weighting matrices, F and G. A convenient way 
for doing this is to start with the identity matrices of state and con- 
trol orders. Each diagonal element is normalized by dividing each by 
the square of the maximum expected value of the corresponding state, 
error, or control variable. These weighting matrices are normally 
diagonal matrices since any off-diagonal terms would lack physical 
meaning. The normalized wel^tlng matrix for the state variables is 
multiplied by one constant and the weighting matrix for the control is 
multiplied by another constant. The ratio of these constants then 
determines the ratio of control activity to state variable variation. 

Once the solution is obtained for several sets of performance 
index parameters, it is necessary to test the resulting system design 
in a nonlinear simulation. The performance of the simulated system 
is compared as the performance index parameters are varied. Thus the 
optimal control approach reduces the problem from trial and error vari- 
ation of control parameters to variation of the ratio between control 
and state variable weighting. The linear model dimensions may also be 
varied for evaluation through simulation. 
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The Multivariable Nyquiat Array (MNA) aathod reaulta in a ayaten 
auch aa ahovn in Figure 5-10. The control would aerve to provide auffi- 
cient decoupling between aodulea to provide a given atablllty wargln. 
Thla function la labeled "partial decoupling" and la derived to achieve 
a property called diagonal dominance. The actuator commanda Into the 
"Space Baae dsmanlca" laay each InfluOiice the orienta.'lon of acre than 
one module. On the other hand, the vector componenta Into the "Partial 
Decoupling" each has primary Influence on a alngle module. Tfuia, once 
partial decoupling is Inaerted, the "Individual Loop Controllera" may 
be each deaigned separately as Independent control loops. This could 
be applicable to the analysis of stablllelng control for the Space 
Construction Base. 


Thus, In the MNA analysis, the selection of control parameters 
is divided into two problems (see Figure 5-11). The first is the par- 
tial decoupling of otherwise strongly interacting loops. This is 
accomplished through the proper selection of the parameters, 1^^. Once 
this is accomplished, each loop is analyzed individually to select the 
parameters, f^^. 

The control inputs to the plant, u^ and u^, each affect both 
outputs, 6^ and 62. Partial decoupling through the proper selection 
of 1^2* ^21’ ^22* links between u^ and *nd 

U2 and That is, uj^ mainly affects 6^ and u^ mainly affects 62. 

The extent of decoupling defined in Figure 5-11 for the 2-loop case, 
is based on a theorem by Gershgorln. 


5-26 



ivnxDv 







loor i 
COCTMt 


oicouptwc mnox 
I or COK7KOL 


co«m 


• - c(S) L r # 


'll«> 

•l,(l)\ 



lll<S) 

Hi($)\ 

21^»> 

«22(S)y 

L OOMINASCt or CL 
SUCH THAT 


Ujl(S) 


U,2(S) 


Hi Hi 


Hi Hi 


Vfttrt: q 


Hi ' ‘u *X2 Hi 
Hi • *11 Hi * *ii Hi 


Hi * *11 Hi * • 
Hi “ *ii Hi * •: 


22 Hi 


22 Hi 


IHi‘*>I * i<22**^l 


(low OOMIJ«AJ.Ct) (COIWC, P0W«A.NC£) 

ro» I . j« ovtR TM mQi/wcY w.vc£ or xictwsT. 

Figure 5-11 - Separation of Decoupling Operation 


REPRODUCIBlLnt OF TOf 
ORIGINAL PAGE IS POOR 









SECTION 6 


6.0 DIGITAL SIMULATION 

Tvo Dodeling Approaches have been adopted for this study. A sinple 
one ditnensional simulation is used for preliminary design development 
and initial prediction of performance. A more detailed, three dimensional 
n-body simulation has also been developed for futher evaluation of ex- 
pected performance. Selected eituations have been simulated cm both 
proRrams and ctxipared. 

The purpose of this section is to describe the one dimensional model 
and some preliminary results from its application. 

6.1 DESOlIPTIC»J OF ONE-DIMEWSIQNAL M3DE1S 

The n-body, one dimensional system is illustrated in Figure 6-1. 

It consists of masses (J, , J-, ... J ) connected by massless rods having 

i. 4 r. 

spring constants (K^, K^, ... ‘nd viscous damping (o^, o^, ... 

Torques are applied to the masses (T^, Tj, ... T^) , which are each a sum 
of disturbance and control torques. This model may be easily adapted to 
represent as many masses as desired. 

nr OF THE 
X IS KK)R 


RFPROttUCIBn. 



Figure 6-1 N Body, One Dimensional Model 
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The equations of notion are: 


J 

n - 


used. 


. -Kj 4 Kj ej - Oj 


J, 6 ^ - Kj 6 j - (Kj 4 K^) 4 Kj 4 6 ^ - (o^ 4 Oj) 6 , 


4 02 f 3 4 T2 


(6-1) 


e ,»K ,6 ,-(K «4K ,)e ,4K ,e 

n-l n-2 n-2 n-2 n-1 n-1 n-1 n 


4 a 


«6 - — (q - 4 o ,)6 -I ^ 4T . 

n-2 n-2 n-2 n-1 n-1 n-1 n n-1 


J 6 - K ,f ,-K ,e 40 ,6 ,-c ,6 4T 

n n n-1 n-1 n-1 n n-1 n-1 n-1 n n 


For this study, a three body one dimensional nodel was extensively 
With n ■ 3 in equations (3-1): 


■^1 h ■ "^1 ®1 *^1 ^2 ■ “l ®1 “l ®2 ^1 


Jj ej “ ®1 ’ ^*^1 **■ ®2 *^2 ®3 ®1 ®1 ‘ ®2^ ®2 


4 02 63 4 Tj 


(6-2) 


J3 63 - K2 62 - K2 63 «2 ®2 - “2 ®3 * ^3 
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A corresponding state variable representation is; 


. Xj, Xj . 


Xj - iTj - Oj (Xj - Xj) - Kj(X - Xj))/Jj 




\ - iTj + a^(X2 - X^) 4 X3) - o^CX^ - X^) ( 3 -: 


- KjCX^ - X3)]/J^ 


3 • S* ^5 “ \ 


X^ - IT 3 4 a^CX^ - X^) 4 K 2 (X 3 - y.^)]/3^ 


In vector fore, equations (3>3) nay be expressed: 


JC*t- 





, d.l 1. iuwtr.t.d in klnck dla- 

one diioensional ■od«* order transfer 

The three bod> one o* ^ second order 

eon. in n.uie ..«•« — 

„e .ncv. co*ln.tion of .ys»n 

.lon.1 nnlal ,p.c. ,.pT.»nt.tlon. for 

. v«l.tv of vl«P=in»- contrnl. A 

le may he used in the app f--„„l8ted, which in vector f 

:r"ui.l coninol .,n.llnn n« n..n lornU.1. . 
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%>here 


■ disturbance torque 

■ control torque 

d.. ■ rate feedback sain 

13 

k.. ■ position feedback gain 



Depending on the control gainsCk^^, d^^) selected, a variety of 
control loop configurations may be realized. 


Figure 6-3 Illustrates the vibrational modes of the three dimen- 
sional model. Dynamic rotation about each of the three axes may be 
studied individually, using the one dimensional model. The two modal 
frequencies associated with each axis were chosen to be representative 
of dynamics of Configuration 1 of the Space Construction Base. 

In order to relate the three dimensional model parameters tn the 
one dimensional model, it is necessary to compute the effective moment 
of inertia of each body or appendage. These inertias are referred to 
the actuation point. The spring and damping constants were calculated 
by deriving expressions for the eigenvalues of of the one dimen- 

sional model, and equating these to the desired values. 


One dimensional model parameters used to represent Configuration I 
are summarized in Table 6-1. 

A five body, one dimensional model has also been developed for 
this study. Figure 6-4 shows various applications of the one dimen- 
sional model. The three body one dimensional model has been used to 
represent each axis of the main body and one solar wing of Configura- 
tion I. The wing is divided into two bodies to represent two modal 
frequencies Four modal frequencies of the solar wing may be 

represented by using the five body version of the one dimensional model 
(6-4b). Also, the five body model may be used to represent the main 
body with both solar wings present (6-4c). 
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Figure 6-4 - AppHcetions of the One Dimensional Model 


Table 6-1 - One Dimensional Model Parameters for Configuration 









These one ditnensional nodels ney also be used to represent the 
large solar panels of Configurations XI, and XII of the Space Construe- 
tion Base. In this case, one or both of the solar panels «ay be repre- 
sented, with the remainder of the structure lumped into the 'main body.' 


OFT*® 
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6.2 SOLAR WING DEXXXJPLING OONTOOL STUDY USING THE ONE-DIMENSIONAL MODEL 


During periods when the Space Construction Base (SCB) Is being 
held at constant attitude, the decoupling control decreases the propa* 
gation of disturbances between the solar wings and main body. The thret 
body one dimensional model has been used for this portion of the study. 
Figure 6-5 Illustrates the decoupling control concept. Momenumi ex- 
change devices In the main body of the SCB are represented by torque 
Tcj. Torque motors located at the hinge point between the main body and 
solar wing are represented by torque T^ 2 * Actuators cannot be placed on 
the outer section of the solar wing (represented by mass Jj), hence 
T ^2 • 0. has been omitted from the Illustration. A constant atti- 
tude is held when and 620 constant. With • Kj >2 “ 1*0. 
complete decoupling of the bodies occurs; however, due to Imperfect 
sensors and actuators, this may not be possible to do in practice. 


The rate and position feedback gains (K‘j, K* 2 » 
calculated assuming that the controlled body is completely decoupled 
from the rest of the structure. A control loop bandwidth of 0.3 
radians/s and a damping of 0.7 were chosen for this control. The gains 
were calculated from: 

Si “ “nl^ *^1 ^■^“‘•B/RADIAN 

51 “ ^'“nl *^1 PT-LB/RADIAN/S 

52 * "N2^ *^2 PT-LB/RADIAN 

K *2 • 2;«„2 *^2 ^->*B/RADIAN/S 


The resulting gains for the torsional, normal, and lateral axis controls 
are summarized in Table 6-2. 
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TABLE 6-2 DECOUPLING CONTROL GAINS 


AXIS 

Si 

Si 

S2 

S'. 

TORSIONAL 

B6449.5 

403431. 

68.152. 

318.045 

NORMAL 

2794.05 

13038.9 

1602.0 

7476.0 

LATERAL 

86449.5 

403431. 

1674.0 

7812.0 


For compdrison purposes, a performance index has been defined for 
each angle of the one dimensional model. 


(Pl)j ■ (®-S) 

0 

where 

(PI)j « performance index for the j*th angle 
a^ ■ value of the j'th angle, in radians 
T ■ duration of computer run, in seconds 
• initial condition of a selected angle, in radians 

For this part of the study, perfect actuators and sensors are as- 
sumed. In further work, realistic models for these will be incorporated 
into the simulation. Disturbances are represented by an initial condi- 
tion on one of the masses (a^^ in equation 6-5} for tht decoupling 
study. 
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Figure 6-8 - Decoupling Performance, Lateral Axis 
of Solar Wing 



1 


•V 

k t 


i 


: t 


In the case of the torsional axis, Jj»J2 so ^^at 

movements of O2 kittle effect on Jj. For this reason, it was de- 
cided that not be decoupled from J2 ■ 0). Decoupling of distur- 
bances propagating from to J2 is provided for, the amount of decoup- 
ling determined by gain K^2‘ ^ series of computer simulations were 
done, varying K^2 0 to 2.0, the results of vd^ich are shown in 

Figure 6-6. 

Performance indices PI2 and PI^ take on their largest value when 
Kq 2 * 0. This is the case where no decoupling torque in applied and the 
greatest disturbance transfer occurs. 

The indices increase linearly with K^2* range of gains 

represents the case where too much decoupling torque is provided. The 
indices decrease linearly until Kq2 • 1»0, where PI2 • PI3 ■ 0, 
corresponding to complete decoupling. As K^.2 is increased from 1.0 to 
2.0, the indices increase linearly with K^2* * range of gains 

represents the case where too much decoupling torque i: provided. 

A similar control was applied to the normal axis. Due to the fact 
that Jj, J2, and are of comparable size, there is a strong coupling 
of disturbances between them. Therefore, it was decided to decouple J2 
from as well as to decouple from J2. A series of computer simu- 

lations were done, varying Kqj and Kq 2 over a range of 0 to 1.0. Tht 
results are shown in Figure 6-7. Good performance is achieved for gains 
from 0.5 to 1.0, indicating that complete decoupling is not necessary 
for this axis. 

In the lateral axis, masses J2, and again of comparable 

size, being the largest. This presents a particularly difficult 
control problem, since J3 cannot be controlled. Simulation of the 
and J2 controls for the lateral axis were done, again varying and 
Kq 2 over a range of 0 to 1.0, the results of which are shown in 
Figure 6-8. Variations in these control gains have little effect on Pl^ 
and PI2. Index PI3 shows only slight changes up to the point where 
Kq 2 • where it sharply increases. It is not understood at this 
time why this is so. The effect of the outer control loops (rate t 
position feedback) on the decoupling control has not been determined as 
yet, and this may explain this result. This work will be done in the 
future. 
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6.3 STODIES WITH 3-DIMENSIONAL I03ELS 


A simulation Is In the development process for representing each 
configuration of the Space Construction Base (Figure 6-9) . It has been 
organized to limit the maximum number of rigid bodies to 20. Except 
for the Research Pallet, the bodies are assumed to be connected by hinges. 
Each hinge has three dimensional freedom in angular rotation. The Research 
Pallet has the added feature of being free to translate Independently in 
three directions. 

The digital computer simulation has been developed for configura- 
tion I and it is in process of final verflciation. Each hinge point has 
been tested separately in the torsional, lateral and normal vibrational 
modes. Applying appropriate initial conditions allowed one to neglect 
the effect of all bodies except the two on either side of any hinge. The 
relative inertias of the two remianlng bodies were forced to differ by 
several orders of magnitude so that one body could be deflected and allowed 
to oscillate about the other. Table 6-3 shows selected analytical results 
from one such oscillation along with the corresponding simulation results 
for a typical hinge check. 

Since there exists a direct relationship between the one dimensional 
and three dimensional models in the torsional mode, comparisons were also 
made to check higher order oscillations in the latter model. A typical 
comparison is shown in Table 6-4. 
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Table 6-3 - Typical Hinge Point Check 


Time 

(Seconds) 

Analytical Solution 
(Radians) 

Simulation Result 
(Radians) 

0.00 

0.0100000 

0.0100000 

0.15 

0.0051247 

0.0051246 

0.30 

-0.0039420 

-0.0039420 

0.45 

-0.0079993 

-0.0079993 

0.60 

-0.0038890 

-0.0038890 

0.75 

0.0033448 

0.0033449 

0.90 

0.0063922 

0.0063921 

1.05 

0.0029400 

0.0029399 

1.20 

-0.0028254 

-0.0028255 

1.35 

-0.0051025 

-0.0051025 

1.50 

-0.0022131 

-0.0022129 

1.65 

0.0023772 

0.0023773 

1.80 

0.0040688 

0.0040688 

1.95 

0.0016578 

0.0016577 

2.10 

-0.0019929 

-0.0019930 




Table 6-4 - Typical Comparison of Results Between the Three 
Dimensional and One Dimensional Models for the 
Torsional Mode 
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SECTION 7 


7.0 


7.1 


APPLICATION OF MULTILEVEL CONTROL TECHNIQUES TO SPACE 
CONSTRUCTION BASE 


INTRODUCTION 

This section begins with a presentation of a state variable 
form of the sixty-six degree of freedom mathematical model, 
comprised of discrete rigid bodies, of the Space Construction 
Base described by Cornell (7-1 )♦ of Bendix Research Lab- 
oratories in his memorandum of December 14, 1977. This state 
variable model was decomposed into a set of decoupled first 
order scalar differential equations to render it more 
amenable to the application of hierarchical multilevel control 
techniques. 

In a prior memorandum by Chichester (7-2) of Bendix Guidance 
Systems Division, multilevel techniques were demonstrated by 
applying them to the optimal control of a single axis torsional 
model to which control had been applied previously by Porcelli 
(7-3) using another method. The overall multilevel approach 
was described in terms of the following sequence of steps. 


1. Express mathematical model of plant, (system to be 
controlled), in state variable form. 

2. Decompose mathematical model of plant into set of 
decoupled equations. 

3. Construct performance index. 

4. Form Hamiltonian. 

5. Develop costate equations with associated coordination 
equations . 

6. Develop control algorithm. 

7. Construct subproblem hierarchy. 

8. Discretize equations of each subproblem in the 
hierarchy. 


♦ These numbers refer to references listed at the end of 
Section 7. 
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This section summarizes the application of the first 
seven steps of this sequence, and an additional step, to the 
twelve configurations of the mathematical model of the Space 
Construction Base documented in the memorandum by Cornell 
(7-1). The additional step cited here is the incorporation 
of special necessary conditions for optimality in the costate 
and control equations due to the non-separability of the 
performance index required for multilevel local vertical 
attitude control of the Space Construction Base. 

7.2 GENERAL DECOMPOSED MODEL 


7.2.1 


background 

Figure 1 and Tables 7-1 and 7-2 are reproduced from Cornell is 
memorandum (7-1). Figure 7-1 presents a topological tree that 
shows how the rigid bodies comprising the mathematical model 
are connected. A single line represents a three degree of 
freedom spring hinge suspension while a double line rep- 
resents a six degree of freedom suspension. Table 7-1 lists 
the number of rigid bodies associated with each configuration 
of the model. Table 7-2 lis';s the numbers of degrees of 
freedom associated with each of the twelve configurations. 

As indicated by Figure 7-2, the equations comprising the 
overall control problem in state variable form may be assembled 
into two related subproblem hierarchies, the translational 
hierarchy and the rotational hierarchy. 


TABLE 7-1 


CONFIGURATION 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


NUMBER OF RIGID BODIES 

5 

6 

7 

8 
8 
8 
8 

10 

10 

19 

20 
20 
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Figure 7-1 . Topological Tree of Space Base Rigid Bodies 




TABLE 7-2 - Degrees of Freedom 


Configuration 

Degrees c 
Translational 

>f Freedom 
Rotational 

Total 

1 

3 

15 

18 

2 

3 

18 

21 

3 

3 

21 

24 

4 

6 

24 

30 

5 

6 

24 

30 

6 

6 

24 

30 

7 

6 

24 

30 

8 

6 

30 1 

36 

9 

6 

30 

36 

10 

6 

57 

63 

11 

6 

60 i 

66 

12 

6 

60 i 

66 


Decomposed Translational State Equations 

The decomposed translational equations for configurations 1 
through 12 may be written in the form of equations (1) 
through (12) with the aid of equations (13) through (43). It 
should be noted that the state variables, and are 

scalar components along the Space Base axes of the translational 
displacement vector of the ith rigid body of the body and the 
state variables ^Tiy ^Tiz corresponding 

scaler components of the translational velocity vector where, in 
this instance, i = 1,8. 


Translational 

0) 

Rotational 

Subproblem 

< 

Subproblem 

Hierarchy 


Hierarchy 


0) = (“i, “2, “20^*^ 


Figure 7-2. Overall Control Problem Structure 
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Decomposed Translational State Equations 
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Scalar State Equations 
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(7) 

( 8 ) 

(9) 
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( 11 ) 
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mm 
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(14) 


SF 


n 

• IT 
i-1 
1^8 


ix 


ZF 


n 

' ZF 
ij»8 


ly 


n 

ZF « ZF 

i-1^^ 

IjtS 


= (^ix. ^iy. ^iz> 


(15) 


(16) 


(17) 


The upper limit of the summations in equations (13) through 
(17), n, varies between 5 and 20 according to the specific 
configuration of the space base being modeled as shown in 
Table 7-1. It also may be seen in Table 7-1, that for 
configurations 1, 2 and 3, the translational equations of 
rigid body 8 are omitted. 


The following coordination equations are defined for each of 
the decomposed translational state equations. 
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(18) 

(19) 

( 20 ) 
( 21 ) 


( 22 ) 

(23) 


®T8x “ ^^T8x 


(24) 
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(25) 
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(26) 


(27) 

(28) 
(29) 


where : 

T 

®18 * ^®18x ^18y ®18z^ vecter representing the 

suspension force exerted on body 1 by body 8. 

T * 

h = (L_ L L ) is the aggregation of terms linear in o). 
X , y , z 1 

resulting from the summation of the translational equations 
of all of the rigid bodies of the model except body 8, 


N N N are components of terms along the space base axes 
X , y , z 

each of which is quadratic in and 


20 

^x “ ^^Txi “"i ^20) 

i=l 


where : 


• 

“i “ 

<"ix, “ly, “lz> 

(“^l) 

^Txi 

T 

“ ^Nxix’ ^Txiy’ ^Txiz^ 

(32) 


20 



Z bm . w. 
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where : 


I 


■ t 

_ f 


T 

bTyi " (b^yix* ^T/iy' Vyiz^ 
20 

^ '^Tzi “i 

i«l 

i ><8 


(34) 

(35) 


where : 
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h V. 

^ “Tziy' °Tziz^ 



(36) 

N = 

X 

20 

£ b) . 

i=l ^ 

i^8 

T 1 

^Txi “i 



(37) 

“i = 

<"ix, "ly 



(38) 



0 
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^Txlyy 

0 

(39) 


1 

^^'"'xixz 
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(40) 
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i^Tyiyz 

(41) 
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( 42 ) 


20 ^ 

Nz = I A^zi “i 

i=l 

ii<8 


Tzi 


®^Tzixx 

0 

^^Tzixz 

0 

^Tziyy 

^^Tziyz 

^®"Tzixz 

^^Tziyz 

0 


(43) 


The subproblem hierarchy corresponding to the decomposed 
t? anslationa! state equations, (1) through (12), and thf 
translational state coordination equations, (19) through 
(29) is shown in Figures 7-3 and 7-4. It should be noted that 
this subproblem hierarchy corresponds to the decomposed 
model without the application of control. 



Figure 7-3. Translational Subproblem Hierarchy 
Without Control 


7-9 




Body 1 Coordination Equations 



k 

Figure 7-4. Body i Translational 

Subproblem Hieraichy for i=l,8 




7,2.3 Decomposed Rotational State Equations 

The decomposed rotational state equations for configurations 
(1) through (12) may be written in the form of equations (44) 
. through (49) with the aid of rotational state coordination 

equations (50) through (55) and equations (58) through (66). 
Here it should be noted that the state variables 0^^ and 
are the Euler angles of the ith rigid body with respect 
to the space base coordinate axes and the state variables 
“iy ^^iz scaler components with respect to 

the space base axes of the angular rates of the ith rigid 

body where li: this instance 1=1,2, ,n, the value of n 

depending upon the particular configuration being modeled 
as shown in Table 7-1. 
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0 

^^Ix, jx.Jy 
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^®ix, jx, jz 
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A ' 

^ix, j 

s 

*^lx, jx, jy 

^ix, jy, jy 

*^ix, jy, jz 

t , 



*®’lx, jx, jz 

^^ix. jy, jz 

®ix, jz, jz 

K 



®^iy, jx, jx 

^®"iy, Jx, jy 

^^iy.jx. jz 

: 

^iy, j 

:s 

^^iy.dx, jy 

0 

^“^iy.Jy, jz 




^^iy, jx, j7 

^^ly, jy, jz 

^iy, jz, jz 




a . 

iz,3x, -*x 

^®^iz, jx, jy 

ia. 

IZ, JX, JZ 


a: . 

iz, J 

= 

^^iz, jx, Jy 

^iz, jy.jy 

4a. . . 

iz.jy.jz 

i 



^^Iz, jx, jz 

^^iz, jy, jz 

0 




~ 

^Ix, jx, jx 

^“ix, jx, jy 

jx, jz 

. 

^ix, j 

= 

*a. . . 

ix,jx,jy 

^ix, jy. jy 

^^ix, jy, jz 




ia . 

1X,JX,J2 

L 

ia . 

ix,iy,iz 

^ix, jz, jz 




^iy,jx,jx 

^^iy, jx, jy 

^^iy, jx, jz 


^iy, j 


^^iy,jx,jy 

^iy, jy, jy 

^^iy .jy, jz 




ia . 

iy,jx,jz 

^^iy.jyjz 

^iy, jz, jz 




^iz, jx, jx 

^^iz, jxjy 

^^iz, jx, jz 


^Iz, j 

= 

jx, jy 

^iz, jy, jy 

4a . 

iz,jy,jz 




ia. . 

1Z,JX,JZ 

ia . 

iz,jy,.iz 

^iz, jz, jz 

S'- 

V 
► 't 




1‘RECEDING page blank not filmed 


(61) 


(62) 


(63) 


(64) 


(65) 


( 66 ) 
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The subproblem hierarchy corresponding to the decomposed 
rotational state equations, (44) through (49) and the 
rotational state coordination equations, (50) through (66) 
is shown in Figures T-S and 7-6. As was the case with the 
subproblem hierarchies in Figures 7-3 and 7-4 for the 
translational equation, these hierarchies also correspond 
to the uncontrolled case. 

Inspection of the state equations (45), (47) and (49) reveals 
that they are quadratic in the state variables and 

^iz Since the traditional development of multilevel control 
techniques has been concentrated on the control of linear 
systems, this form of the state equat lore poses a special 
problem in such application. Pontryagin's maximum principle, 
upon which multilevel optimization is based, does not 
preclude the control of nonlinear systems and some recent 
papers by Hassan and Singh (7-4), (7-5), (7-6) and (7-7) describe 
some potentially useful approaches for effecting such control 
by multilevel techniques. 
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Hierarchy 



(«> 


n 





Body n 

Subproblem 
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Figure 7-5. 


Rotational Subproblem Hierarchy Without Control 








7.3 


CONSTRUCTION OF PERFORMANCE INDEX 


7.3.1 General Discussion 

After the state equations for the uncontrolled system 
(plant) have been decomposed, i.e., decoupled, the next 
step in the application of hierarchical multilevel optimal 
control is the formulation of a performance index functional 
that provides a measure of the departure of the system's 
operation from the optimal condition. The most distinctive 
characteristic of a performance index utilized in the 
application of multilevel control techniques is that it 
must be either decomposed at the outset in the same manner 
as the state equations of the plant or amenable to such 
decomposition . 

The specific form of the performance index and the variables 
utilized in it depend upon the type of control to be effected 
and the relative importance of such factors as the polarity 
of errors and minimization of the energy expended in effecting 
control in comparison with minimization of error. In the 
present section the formulation of a performance index for 
local vertical attitude stabilization involves only the 
rotational state variables c r those rigid bodies of the 
mathematical model to which torque is applied. 

7.3.2 Performance Index for Local Vertical Attitude Stabilization 
Since the performance index of this section pertains only 
to attitude control, it is not a function of any of the 
translcitional state variables. The following assumptions 
were utilized in constructing this performance index. 

1) The differences (errors) between the actual and 

♦ ★ 

specified Euler angles of body 1, 

and I respectively, are to be minimized. (The 

convention adopted hei v and later in this section is 

♦ 

that X is the specified or optimal value of x.) 
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2 ) 


The differences between the actual and specified 

body rate components of body 1, - “tv** “ 1 ,,““* 

^ 1 X IX X y 

and are also to be minimized. 


ly 


* 


3) The Euler angle errors (j>g- (j>g*, 0g- 0g* and ipg- 

and bod^ rate errors, Wg^^- Wg^*, Wgy- “gy* 

Wg^-Ug^*, are similarly to be minimized if body 8 is 
present in the configuration being controlled. 


4) The relative Euler angles and body rate components 

with respect to body 1 of the remaining bodies comprising 
the mathematical model are to be minimizei*. 


5) Departures from specified or optimal values of the 
state variables are equally important for either 
polarit, of the same magnitude. 

6) Control is effected only by the application of torque 
to each of the rigid bodies of the model, 

7) Control actuator torque is applied to each solar wing 
only at its interface with the rest of the space base. 

8) Control energy is proportional to the square of the 
torque applied to each rigid body of ihe model. 

9) ^,ach error term and each control energy term is 

multiplied by a constant coefficient, W. . representing 

1 > J > 

the rela ive significance of that term in the performance 
index. 

10) The time interval over which performance index is 
to be optimized extends from t to t. 
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The performance index resulting from the application of 
the above assumptions is the following: 


i! 


^ i 


•V 


■/ 




I 


i=2,3,6,ll 


"»1.7 (V**! 

[*i , 1 <♦ i-*l , 2 . 3\?x 

, 4 ( e 1-9 , )2*Wi ^ 5(c-iy-a.iy )'«i , gT, ?y 

«i y(»i-*i ^ _ gT^i^ ] 

*’'7,i(*7-*6>^*"7,2(“7x-“1x>H, 3 
+*7,4(07-«6’H,5<"7y-“ly>^«7,6 ^a^y 
«7.7(1'7-'l'6)^*»7,8(“7z-“ez>^«7,9 

"*8,l<*8-*8 )H,2<“8 x-“5x>H, 3 ^aL 
«8,4‘ )^**8,5<“8y-"8PH,6 ^Jy 

«8 , 7<1'8-1’S )^^*8 , 8<“8z-“8^"^*8 , 9^alz 
* ^ [’*i,l<*i-*ll>H,2<“ix-“llx>H,3TaL 

1-12,13 * »i,4(ei-hl)H,5'“ly-“uy)^**i.6’'a?y 

*»l,7(*r1’ll)°"’'l,8<“iz-“llz>H,9’''alz] 
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( 67 ) 


«l 0T 


I 


i-9.10.20 , 


(67) 
Cont 'd. 


*»1.7<»l-1'7>H.8("lz-“7z)^**i,9’'fiz' 


> dt 


A similar performance index could be constructed for 
inertially referenced optimal attitude control. It would 
differ from that shown in equation (67) in that ((>|=0|=O and 
“lx*“iy*“lz*®* Neither of these performance indexes is 
separable in tiiis form due to the generation of cross 
products of the state variables by the error terms. Hence, 
neither of them can be utilized directly in the application 
of multilevel control techniques. A paper by Singh and 
Hassan (7-7), however, presents an approach to multilevel 
optimization that utilizes non-separable performance index 
or cost functions that are transformed to quasi-separable 
forms. The corresponding form of the local vertical attitude 
control performance index of equation (67) is presented in 
equations (68) through (118). 

/ f f20 r *2 ♦ 2 

j 4lx<faix-^alx> *“ly<’'aiy-\ly> 

to 1*4,5,14,15,16,17,18,19 
"“iz<’'alz-’'Jlz>"*JjPl,j] 

*®«!,«I,*!,“lx,“Iy,“!z 

The terms and coefficients introduced in equation (68) 
are defined in equations (69) through (118) 


i 



i 



1- 


* ♦ ♦ V 

1*1,2,3,6,11 


f. ^ 


+W 


★ ♦ 


♦ ♦ 




1.5 “lAy *''l, 7 'l'l 1 'i **l, 8 "lz“lz' 


♦ ♦ 


♦ ♦ 


* 4c 


+w? i(|)g (|>7 +W 7 ^ 2“6x‘^7x'^^7 , 4 ® 6 ® , S“6y“7y 
** 7 , 7*^7**7 , 8“6z“7z-^*8 . 1 ' ^^*8 , 2 ' “8x > ^ 

■^*8,4*®8^ ■''*8,5*“8y' **8,7**8^ ** 8 , 8 *“Sz*^ 


♦ ♦ 


♦ ♦ 


♦ 4c 

’1 


* 2 <*l,l*7*i**i,2“7x“lx"*'l,4878: 

i-9,10,20 

■"*1 , 0“7y“iy-"*i ,7*7>''i*'*i“s“7z“lz> 


★ ♦ 


^(W. ^ ^ ow-i -, „ 0),; ^+w, ^ e, ;e: 


ouiin (0* ^Wi x6i^0‘ 

2 llx IX 1,4 11 1 


+ 

1=12,13 

, 5‘^lly‘*‘iy'^^i , 7’^ll’^i‘^^i , 8“llz“iz 




l, 2 =^, 2 “l^'^l, 2 <^^"^l, 3 Tlx" 


Pl, 3 =^ 4 ®l ^^. 4<®1 


Pl, 4 "^l, 5 “ly'^''^ 3 , 5 ^“ly^ '^'’^ 1 , 6 ’^ly 


2 , 


Pl. 5'^,7 > 


Pi, 6 =^ 1 , 8 “iX, 8 (‘^ 1 z>'^^ 1 , 9 Tiz 


(69) 

(70) 

(71) 

(72) 

(73) 

(74) 

(75) 


For 1=2,3,9.10,12,13,20: 


k 

i 
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V 


'’l.2'*l,2f^?x**i,3Tlx 




Pi, s'*!, 7* t' 


Pi,6'*i,8“i>i,9’'iz 


For i« 6 , 7 ,ll: 


(77) 

(78) 

(79) 

(80) 
(81) 
(82) 


(83) 

(84) 


Pl, 3 **i, 48 i" 


(85) 

Pi,4-*l,5“lH,6’'iy 


(86) 



(87) 

Pi. 6 ** 1 , 8 “iK, 9 Ti^ 


(88) 

r ^ * 

! 

(89) 

P 8 , 2 -* 8 , 2 [“ 8 ^«^x> 1 ** 8 , 3 T 8 x 

(90) 

P 8 . 3 -* 8 , 4 [® 8"*<<>1 


(91) 

P8,4-*8,5 [“8^<“8>' 

1**8,6’'8y 

(92) 

P8,5‘*8,7[^'8^*<''’^^] 


(93) 

P8,6“*8,8 [“8z**“8z^^] 

2 

+V/ T 

1 ”8,7 8z 

(94) 


t ; 

i 
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I 


ii 


: 1 


1 

i 

f ; : 


(95) 

! f * 

» 


(96) 

\ 

^1 , 4*^1 , 4*^2 ,4*^,Z, 4^^e , 4‘*‘"'ll , 4 

(97) 

► * 

W.5*^1.5"^2,6"^3.5"^6.5"'^ll,5 

(98) 

• * 

*1 , 7*^1 , 7‘^’'^2 , , 7‘*'^6 , 7 "'^11 , 7 

(99) 


<.8*^.8^’'2.3"^3.8*^6.8"^1,8 

(100) 



(101) 


^6,2“^6,2^^7,2 

(102) 


^6,4*^6,4'*‘'^7,4 

(103) 


^6,5‘^6,5'^^7,5 

(104) 


^6,7“^6,7'^^7,7 

(105) 


oo 

+ 

00 

II 

00 

"sP 

‘ ( 106 ) 


^7, l“^7, 1'^^IO, l‘^^20, 1 

(107) 


^7.2'^7.2^«9.2^^1G,2^^20.2 

(108) 


''^7 , 4 “''^7 , 4 ■^''^9 , 4 0 , 4 ■*■^20 . 4 

(109) 


^7 , 5“^7 , 5"^^9 , S'^^IO , 5‘^^20 , 5 

(110) 

^ 

^ •- 

^7 . 7“^7 . 7'*’'’^9 , 7 ’■"^1 0 , 7'^''^20 , 7 

nil) 

^ Ii 

^7 , 8'“^7 , , 8'*‘^10 . 8'^'*^20 . 8 

(112) 


i* 


iflSi?'' 





U, 1*^11, 

(113) 

11,2“^11 ,2'^^12,2‘^^13,2 

(114) 

11;4''^11,4'*'^12,4'^^13,4 

(115) 

11,5‘'^11,5‘*‘^12,5'^^13,5 

(116) 

1 i , 7*^1 1 , 7‘*'''^1 2 , 3 , 7 

(11/) 

11,8 “^ll,8'^'*^12,b '^IS.e 

(138) 


-t 
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7.4 


FORMATION OF THE HAMILTONIAN 


7.4.1 General Discussion 

In order to utilize Portryagin's minimum principle in the 
application of optimal control to a plant, it is nec- 
essary to construct the Hamiltonian of the plant from the 
plant's state equations. If the equations descrloing 
the plant are decomposed (decoupled), the Hamiltonian is 
cons j. acted from the terms of the performance index, the 
right, hand side of the state equation and the coordination 
equation associated with each of the subplants defined 
by the particular decomposition used. 

7.4.2 Hamiltonian Corresponding To Local Vertical Stabilization 
Control . 


20 6 

H= I T. p. j +G(<j.| 6* J,* 

i=l j=i 

i?*4, 5, 14, 15, 16, 17, 18, 19 

20 f 

+ V <u. (T . -T * )^+u. (T . -T ! 

^ 1 ix' aix aix ly aiy aiy iz aiz aiz 

i=l L 

iM, 5, 14, 15, 16, 17, 18, 19 




a 




■i| 


■^^i,3®Riy‘^^i,4 


The subscript 
equilibrium. 


. b) . ) u . 

IX, IX, IX ix'o IX 

eix,^ix, ®Rix^| 
iy,iy,iy ly'c ly 

”o" denotes evaluation at a point of 
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’‘iiiy'^Rly^'''ely/ty.®Rly’ 


Hi I 

JK 

i, 5 ®Riz 


0) . ) 0}j 


+(14*) T«4 F., 


1^33 *aiz ”Riz'‘*eiz, ‘iz, Riz^ 


*‘’l.l<“ix*®l"iz-®Rlx> 


20 T 


*‘’l, 2 [“I*ix,i“i* /.i<“j*ix.j“j^'’lx.j“j^ 

Jl'i 


'Ii> 1 

iij ii;i •' 


*‘^l,3<“iy-*i“iz-®Rly> 


[ rp » 20 T~ * 

0 ). A. .(i5.+ Z (oi. A.„ .o).+b.„ 4 W.) 

1 iy,i 1 ' j iy,j j ly.j J 

j=l 

j=/i 


I I 


*fl, 5 (*i“.y*“lz -®Riz> 


20 


‘’l.eKX'z.iV ^ <“j\y,j“j'"'’iy.j“j> 

j = l 


j^i 


*[hha hx » “i l2^ Tiy*8Ri2,-Sgia] 
li'l nil 


.T'iT. 
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Z V.; 


u 

III 


I 




!.] 


U1 

h 


r i 

? * 

i ; 


■^'’l.l*’'aix"’'aix^*'’i.2*’'aly'’'aiy^*''i,3^’'aiz"’'alz^ 


(H9) 


where j and j ,j * 1,2,..., 6 are Lagrange multipliers 
introduced to ensure satisfaction of the following conditions: 


“lx°??x 

X X^ 

“iy““iy 

(A) . =(*).♦ 

IZ xz 

T =T * 
iax iax 

T. =T.* 
lay lay 

T. =T.* 
laz laz 


( 120 ) 

( 121 ) 

(122) 

(123) 

(124) 

(125) 

(126) 

(127) 

(128) 


and X. .. , j = l,2,...,6 are the costate variables 

X , j ) 


The decoupled rotational state equations, (44) through (49), 
may be obtained from the following necessary conditions for 
optimality applied to the Hamiltonian. 


. ! 

‘ * i 
{ I • 







ax 

i.2 


0.= 3H 

^ ax 

i,3 


7-'>7 


Q29) 

(130) 


( 131 ) 



( 132 ) 


(133) 


(134) 


The rotational state coordination equations, (50) through 
(55), may be derived from the application of the following 
necessary optimality conditions to the Hamiltonian. 


^ * 0, i»l,2, ... ,6 


(135) 


Development of Costate Equations 

The rotational costate equations, which are equal in number 
to the rotational state equations, may be obtained by 
application of the following necessary conditions for 
optimality 


'i.r 


(136) 


(137) 


(138) 
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3(A), 


X. c=-3H 




( 139 ) 


(140) 


(141) 


From equations (70) through (94) and (119) the above listed 
neressary optimality conditions yield the following: 

For i = 1,6,7,11: 


l‘ -2*i, , 3-“iyPi , 5 - 81 , 1 


(142) 


For i = 2,3,8,9,10,12,13,20: 


*1.1' 5-81,1 


(143) 


Xi_i(tf) = 0 


(144) 


For i = 1,6,7,11: 


^i ,2 ~^^i, 2 *^ix ^^^ix, ix, ix^ix^o"^ i , 2 ~^i , 1 


-fa. . .0). +a . . . (a) . ) 0, « 

' ix,ix,iy ly ix,ix,iz iz'^o^i,2 


‘^^^iy, ix, ix“ ix^^iy, ix, iy^iy^^iy, ix, iz‘^iz^o^i,4 


~^^®^iz, ix, ix*^ix^^iz, ix, iy ^iy^^iz, ix, iz***iz^o *^i ,6 


(145) 
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For i - 2,3,8,9,10,12,13,20 


, 2 " , 2‘^lx“^^ ^Ix, ix, ix“ix^o^i , 2” , 1 

“^^Ix, ix, iy"ly*^x, ix, iz“iz^o ^i, 2 

ix, ix^^ix^^iy, ix, iy^'^iy^^iy, ix, iz^iz^o *^i,4 


■^^®^iz, ix, ix“ix'*’®’iz, ix, iy“iy"*’®’iz, ix, iz“iz^o ^i,6 




For i = 1,6,7,11 


^i,3' “^''^i,4®i"“izPi,l~^i,3 


For i ’* 2,3,8,9,10,12,13,20 


^i,3” -2’^i,4®i-‘^izPi,l“^i,3 




For i = 1,6,7,11 : 


^i,4* ■^^i,5"iy"^^^iy,iy,iy“iy^o^i,4 


■^^ix, ix, iy^lx'^^^ix, iy. iy^iy'^^ix, iy, iz^iz^o ^i ,2 


*^i, 3“^^iy, ix, iy**^ix^^iy, iy, iz“iz^o ^i,4 
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( 151 ) 


ix, iy, ly^iy'^^iz, iy, iz“iz^o^i,6 

For i = 2,3,8,9,10,12,13,20; 

, 4 “ , 5“ly“^^^iy , iy, ly“iy ^o^l,4 


-(a. j . u). +2a. . . w. +a. j j u)- ) f*. o 

' ix,ix,iy ix ix,iy,iy iy ix,iy,iz iz‘'c x,2 


ix, iy“ix^®^iy, iy , iz“iz^o^i,4 


~‘*'i^i,5"^i,4 


-fa. j . w. +2a . . . u) . +a. . . u). ) P. 

' iz,ix,iy IX iz,iy,iy ly iz,iy,iz iz‘^o 1,6 

'i,4^^f> = 0 

For i = 1 , 6 , 7 , 11 : 

^i,s“ ■^^l,7'*'i~^i,5 

For i = 2,3,8,9,10,12,13,20 

Ai,g=* 


^,5 ° 


(152) 

(153) 

(154) 

(155) 


(156) 



For i • 1,6,7,11 : 


^i,6* “^^i,8“l2"^^^iz,iz,iz“lz^o^i,6“®i^i,l 

”^*lx, ix, Iz^lx'^^ix, iy, lz“ly‘^^^lx, iz, iz“iz^o^i, 2 


■^®iy , ix, iz“ix‘^®^iy , iy , iz“iy‘^^^iy, iz, iz“iz ^0^1,4 


■ ^l,5"^i,6 

“^®’iz, lx, iz‘^ix^^lz,iy, iz“iy^o^i,6 
For i = 2,3,8,9,10,12,13,20: 


,6* , 8“iz~^^^iz, iz, iz“iz^o^i, 6"®i Pi , 1 


~^®^ix, ix, iz“ix^®^ix,iy, iz“iy^^®^ix, iz, iz“iz^c^l,2 


~^^iy, lx, iz^ix^^iy, iy , iz“iy^^®^iy , iz, lz“iz^o^i,4 
"Pi,5"^i,6 

”^^lz, ix, iz^ix^^^iz, iy, iz“iv^oPi,6 


(157) 


(158 

( 159 ) 
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The rotational costate coordination equations, which 
are equal In number to the rotational state coordination 
equations, may be derived from the following necessary 
optimality conditions 


38 = 0 
Rlx 


3H “ 0 
»Sr1x 


38-0 


M = 0 
*®Rly 


M • 0 
^ Riz 


M “ 0 

'SrIz 


Pi,l “ =^i.l 

^ = X 

1,2 '^1,2 

*^1,3 " ^1,3 

^1,4 “ ^i,4 

Pi, 5 “ ^i,5 

Pi, 6 “ ^i,6 


( 160 ) 

( 161 ) 

( 162 ) 

( 163 ) 

( 164 ) 

( 165 ) 
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DEVELOPMENT OF CONTROL ALGORITHM 

The optimal control equations resulting from a gradient 
approach to the local vertical attitude stabilization 
control problem yields the following equations: 




( 166 ) 




(167) 


- «iz <1!— >r 
^ ^ aiz 


(168) 


where r is the iteration index and q^^, qz„ and q. are 

xx X y X <6 

constants to be chosen on the basis of the rate of 
approach of the control system to optimal operation. 


an 


3T 


aix 


=2W. oT. +2u. (T. -T*. ) + X, o 

1,3 IX IX ' aix aix'^ | 


. ('i’>12 , ('i’>13 

|ii' I ^1.4 TvT 


Pi, 6 ^i,l 


(169) 


^ * ^22 

IIt:..' 2Wi,6Tly*2“ly (^aiy-^liy) * 717^- 'i,4 


iay 


jT^r' ^ 


i,2 Tlfn ^i.6 '^i,2 


(170) 
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D 

U 

u 

\ 1 


3H 

3T 


(I. ') 


- (Talz-’'Jlz) * TTTI^ 'i. 


laz 


11 


^32 


^31 

* ‘>1,2 * TT7r »i,4*^i,3 


i 7.7 ADDITIONAL NECESSARY OPTIMALITY CONDITIONS 


3H 


36 


i.l 




3H 


36 


i.2 


CO 


lx 


(0 . 
IX 


3H 




1,3 


= 0 -*• 0* * 


3H 


36 


i,4 


(0 . * CO . 

iy ly 


3H 

■yl 


i.5 


= Vi 


3H 

36 


— = 0 

1,6 


“iz “ “iz 


I 

■% * 


3H 


3v 


1,1 


T * T 
aix alx 


( 171 ) 

( 172 ) 

( 173 ) 
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7.8 CONSTRUCTION OF SUBPROBLEM HIERARCHIES 

In order to effect optimal control using multilevel 
hierarchical techniques selected subsets of the equations 
presented thus far are assembled into subproblems which, 
in turn, are assembled into a hierarchical configuration 
similar in form to that depicted in Figure 7-5. The 
subproblem hierarchy corresponding to optimal control of 
the system incorporates two additional classes of sub- 
problems to be solved on the lower level, the costate 
subproblems and the control subproblems. The addition 
of the costate subproblems to the lower level of the 
hierarchy increases the number of coordination equations 
appearing in the coordination subproblem at the apex of 
the hierarchy because each of the costate variables 
has a corresponding costate coordination variable associated 
with it. To simplify subsequent discussion of this 
hierarchy all of the subproblems of a particular type will 
be regarded as being grouped into a single subproblem, that 
is, all state subproblems into a single state subproblem, 
all costate subproblems into a single costate subproblem 
and all control subproblems into a single control problem. 
Then each subproblem of a particular type will consist 
of all equations of that type with their associated 
boundary conditions. 

In particular, for optimal local vertical attitude stab- 
ilization, the state subproblem consists of equations (44) 
through (49) which are the decomposed rotational state 
equations of the Space Construction Base. The decomposed 
translational state equations, equations (1) through (12) 
also may be included in this subproblem although they are 
not essential for the feedback control ..quations in this 
case. 
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The costate subprobletn for optimal local vertical 
attitude stabilization consists of equations (142) 
through 159) and the control subproblem consists of 
equations (169) through (171). 

In addition to the coordination equations for the decomposed 
rotational state equations, equations (50) through (55) 
and the coordination equations for the decomposed trans- 
lational stable equations, if they are included in the 
state subproblem, the coordination subproblem at the apex 
of the hierarchy also includes the costate coordination 
equations, equations (160) through (165) in order to 
effect optimal local vertical attitude control. 

It should be noted that the costate equations, (142) through 
(159), incorporate the coordination variables, 3. ., 

1 > J 

and the control equations, (169) through (171) incorporate 
the coordination variables 2 '^i 3’ 

additional equations required to define these variables 
are provided by equations (69), (182) through (19 > and 
equations (200), (202) and (204). These additional 
equations may be incorporated in the coordination sub- 
problem at the apex of the subproblem hierarchy. This 
increases the overall dimension of the coordination sub- 
problem to a considerable extent and also greatly 
increases the number of coordination variables to be 
transmitted between the subproblems of the hierarchy 
with control compared with the number required for the 
hierarchy without control. The resulting subproblem 
hierarchy for the controlled system is portrayed in 
Figure 7-7 in which the following relationships are utilized. 



STATE 

EQUATIONS 



CONTROL 

EQUATIONS ^ 


COSTATE 

EQUATIONS 



Te 


Subproblem Hierarchy For Multilevel Local 
Vertical Attitude Stabilization Control 
Without Translational Equations 


FIGURE 7-7 
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(212) 
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(213) 
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(214) 
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(215) 
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(216) 

3='^^- ft ft 
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(217) 
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(218) 

'’ ' f''l’ '’2 Vjq)’' 

(219) 
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( 221 ) 




e20^ 


Ti = 


*^ei ~ ^"^eix, *^eiy, *^eiz^ 


( 222 ) 

(223) 


It should be noted that represents the external 
disturbance torque applied to the ith rigid body of the 
model and Tg,i represents the control actuator torque 
applied to the same body. 

If the translational state equations, (1) through (12), 
are added to the state subproblem and the translational 
coordination equations, (19) through (30), are added to 
the coordination subproblem of the hierarchy of Figure 
7-7, it issumes the form shown in Figure 7-8. For this 
hierarchy, the following additional relationships are 
required. 


iji ** ( ' ^T8 ^ 


(224) 

Ti " ^®ix’ ^Tix’ 

®ly ^Tly. Hlz' 

(225) 

A* 

iji “* ( > ^T8 ^ 


(226) 

/N 

Ti “ ^®Tix’ ®Tiy’ 

^Tiz^ 

(227) 

Ip ^ ( ®»p j 9 ^T8 ^ 


(228) 

Ti “ ^®Tix’ ®Tiy’ 

StIz)^ 

(229) 
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Subproblem Hierarchy for Multilevel 
Local Vertical Attitude Stabilization 
Control With Translational Equations 


FIGURE 7-3 
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The coordination variables B and v and their corresponding 
equations may be eliminated from the subproblem hierarchy 
shown in Figures 7-7 and 7-8 by expressing each of them 
in terms of the remaining coordination variables. Each 
of the resulting modified subproblem hierarchies for 
the controlled system contains a coordination subproblem 
of substantially smaller dimension and requires 
transmission of many fewer coordination variables than do 
each of the original hierarchies. This approach to the 
reduction of the subproblem hierarchy consists of several 
steps. 

1. Expansion of the partial derivatives of G appearing 
in equations (182) through (198) defining Bj^ j by 
utilizing equation (69) which defines 

G(PJ. wjy, i|)|, 

2. Substitution of the results of Step 1 in equations 
(182) through (198). 

3. Substitution of the results of step 2 in the costate 
equations (142) through (159). 

4. Substitution of the equations defining (200), 

(202) and (204) into the control equations, (169) 
through (171). 

The costate and control equations resulting from applica- 
tion of the steps listed above are presented in Appendix C. 
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Substitution of the equations of Appendix C for the 
corresponding costate and control equations in the sub- 
problems of the hierarchy depicted in Figure 7-8 leads 
to construction of the reduced subproblem hierarchy 
shown in Figure 7-9. The definitions of the variables 
given in equations (205) through (223), except for the 
variables that have been eliminated, also apply to the 
reduced hierarchy. 

For the most general Space Construction Base configuration, 
12, the range of the rigid body index, 1, is given by: 
1=1,2,3,6,7,8,9,10,11,12,13,20. The range of i is 
reduced correspondingly for the 11 other configurations. 
Table 7-3 lists the numbers of scalar state, costate, 
coordination and control equations to be solved for 
each configuration in order to effect multilevel local 
vertical stabilization attitude control utilizing the 
reduced subproblem hierarchy of Figure 7-9. 

Inspection of the state equations, (1)-(12) and (44) 
through (49) and the costate equations of Appendix C, 

(C-1) through (C-48), reveals that they are ordinary 
first order differential equations. In general, initial 
conditions are known at the outset for the state 
variables, while final conditions are known for the co- 
state variables. The remaining equations in the sub- 
problem hierarchy, the state coordination equations, (19)- 
(30) and (50) through (55), the costp.te coordination 
equations, (160) through (165) and the control equations, 
(C-49) chrough (C-51) are algebraic. 




Reduced Subproblem Hierarchy for 
Multilevel Local Vertical Attitude 
Stabilization Control With Translational 
Equations 


FIGURE 7-9 
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TABLE 7-3 


Numbers of Scalar Equations To Be Solved For 
Local Vertical Attitude Stabilization Control 


CONFIGURATION 

STATE 

CO-STATE 

COORDINATION 

CONTROL 

TOTAL 

1 

36 

24 

60 

9 

129 

2 

42 

30 

72 

12 

156 

3 

48 

36 

84 

15 

183 

4 

60 

48 

108 

18 

234 

5 

60 

48 

108 

18 

234 

6 

60 

48 

108 

18 

234 

7 

60 

48 

108 

18 

234 

8 

72 

60 

132 

24 

288 

9 

72 

60 

132 

24 

288 

10 

90 

66 

156 

33 

345 

11 

96 

72 

168 

36 

372 

12 

96 

72 

168 

36 

372 
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Due to decomposition, the coupling between the equations 
associated with the ith rigid body and those associated 
with the remaining rigid bodies in the mathematical 
model is suppressed. The state and costate equations 
associated with the ith body, however, constitute a 
two-point boundary value (TPBV) problem to be solved 
because the initial values are known for the state 
variables and the final values are known for the costate 
variable at the outset. Such a set of equations usually 
is discretized temporally to approximate each ordinary 
differential equation with a set of finite-difference 
equations that is more amenable to numerical solution on 
a digital computer. Once this step is taken, all of the 
equations in the subproblem hierarchy are then algebraic. 
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Section 8.0 


8.0 Structural Analysis Task - A Method for Determining 
Overall Flexible Body Characteristics for a Series 
of Connected Substructures* 

8.1 INTRODUCTION 

This section summarizes two analytical methods by which 
a series of flexible bodies can be connected to each 
c^-her. Both of the techniques presei'ed assume that 
the flexible body characteristics of the indi\idual 
sub-bodies comprising the total structure are known 
in terms of their free- free modes, assuming a discrete 
coordinate formulation. The first technique assumes 
that a single connection point existr between two ad- 
jacent bodies thus forming a topological tree configura- 
tion without any closed loops. This formulation does 
not require detailed modeling of the dynamics at the 
connection point and is particularly useful in modeling 
systems that have large relative angular degrees of 
freedom between adjacent sub-bodies (i.e., articulated 
antennas, etc.). The second technique allows multiple 
connection points between adjacent bodies and hence 
would also be applicable to a vehicle topological tree 
configuration having closed loops. In this formulation, 
detailed modeling of various connection points is required 
in order to specify the interface forces and torques 
between contiguous bodies at these points. Both modeling 
techniques are presented in terms of an illustrative 

♦Surarotrized from S. C. Rybak, ”A Method For Determining 
Overalx Flexible Body Characteristics For A Series of 
Connected Substructures," AAS 78-104, March 10-13, 1978. 
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example of a structure composed of five sub-bodies. 

With proper parameter valu?s the configuration used can 
represent vehicles that have already flown (i.e., 

Saturn, Skylab) for which the characteristics of the 
overall structure and the individual sub-bodies com- 
prising the structure are known. These data can then 
be used to determine the fidelity with which each of the 
two techniques presented models the specific vehicles 
consi ;ered. 


BACKGROUND 

To meet the requirements of large orbiting space stations 
and solar power satellites, NASA is considering a 
generation of satellites that are many times larger than 
anything considered or flown to date. Attitude control 
of such large satellites presents a unique and challenging 
problem to the control engineer, no doubt requiring the 
development of new and novel control techniques to 
achieve overall satisfactory satellite control. However, 
before control system design can proceed, a model of the 
vehicle to be controlled must be generated in order to 
allow the formulation of satisfactory control tech- 
niques. It is virtually a certainty that for the size 
of satellite being considered for a space base or solar 
power station, structural flexibility and its interaction 
with the on-board cont^'ol system will be a prime considera- 
tion in its design. Due to the vehicle size overall 
vibration testing of the total structure in order to 
obtain its flexibility characteristics is not feasible. 
Modeling the total vehicle using a normal modal-coordinate 
approach, although possible, has some serious drawbacks. 


although the conunents made about the truncation possibility 
of the technique presented below are plausible, at the 
present state of development, they are still conjectural 
and additional work is required to establish their 
validity. 

Equations of motion assuming single and multiple connections 
between contiguous bodies are presented. Single connection 
points between bodies do net require the modeling of 
individual connection points and hence the equations 
generated will apply equally to a series of connected 
rigid bodies or rigid bodies connected to flexible bodies. 
However, the system of equations is restricted to des- 
cribing a topological tree configuration without any 
closed loops. The equations developed for multiple 
connection points between bodies will be applicable to 
topological tree configurations containing closed loops. 
However, in the manner in which they are formulated, 
they would not be applicable to the accommodation of two 
or more rigid bodies connected to each other within the 
connection chain. 

Equations of motion for single and multiple connections 
between bodies are developed for a vehicle consisting 
of five substructures. With proper parameter values 
the structural configurations used could represent 
vehicles such as Saturn or Skylab for which detailed 
structural data are available. The analysis techniques 
presented can be applied to these data thus determining 
the degree of modeling fidelity achieved. 

MGE bunk now Htv 
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8.3 EQUATIONS OF MOTION FOR SINGLE CONNECTION POINT 

BETWEEN BODIES 

Equations of motion were generated for Figure 8-1 shown 
below. The general procedure that was followed is de- 
scribed in tabular form and the resulting equations are 
identified. The detailed derivation and accompanying 
assumptions may be found in the complete paper described 
in the reference. 

Translational Equations of Motio n 

a) Generate equations describing the sum of forces on 
each body. 

b) Solve for the constraint forces between bodies. 

Angular Momentum of Flexible Body 

a) Derive the expression for the total angular momentum 
of the configuration shown in Figure 8-1. 

Rotational Equations of Motion 

a) Derive the expressions for the rotational equations 
of motion using the rate of change of angular mo- 
mentum expressions for each of the five bodies in 
the figure. 

b) Solve for the constraint torques between bodies. 

c) Substitute the constraint torque equation expressions 
in the rate of change of angular momentur. expression 
for body 2 and derive an overall rotational equation 
of motion for the vehicle system. 

d) Express all terms of the overall rotational equation 
of motion in body 2 coordinates. 
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b) Modal Equation Set 


. ITN p . . IRN . .ITN, . . IRN. l.„ 




(8-9) 


^ 2 f ^ 2 e ^ 2 fT 2 e ‘ <^^^(-^^^012 ^ '^^^<^ 3)^028 


+ ^2™(r2)P^24 " 


2^ cl 2 


+ rf.2RN, 2.„ . ,2RN. 2.„ . .2RN,2.„ 

^ ^^3^^c23 ^ ^ ^^8^^c24 ^ ^^6^^c25 


if2^2 , p2»2 , jf2 2 


(8-10) 


8-9 


tV ' 


I 


^3TNp . .3RN-, .3TK/ ^3.« .3RN. ^3x„ 

he he ^ he he “ (~h^h23 “ ♦ ("''4>**c23 “ 


“n^N ^ ^ 


( 8 - 11 ) 


4TNp . .4RN„ .4TN 4 „ 4RN. _4.„ 

^4e ^4e ^ '**4e he " ^ '‘"^Q>^r24 ' 


9^*c24 


9"“c24 


1*4-4 ^ ^4*4 ^ ^4 4 
^N% 


( 8 - 12 ) 


.5TNp . ,5RN„ .5TN, 5_ ,5RN. 

'^5e ^5e ^ ***5e ^5e " ^ ~ ^ 


7''*c25 


7^"c25 


^ * 44 


(8-13) 


c) Transformation equation set: 
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d) Overall rotational equation of motion: 
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e) Modal equations with substitutions for constraint forces 
and torques. 
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The necessary set of equations is: 

a. Rate of change of angular momentum for the system. 
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b. The transformation equations as defined in (8-14) 
through (8-19). 

c. The interface force equation in "j" rigid body 
coordinates. 
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d. Interface torque equation in "j" rigid body 
coordinates. 
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e. System rotational equations of motion. 
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f. System modal equations. 
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NOMENCLATURE 


8.5 


R 


r 


2 

2 


r 


2 

3 


r 


3 

4 


r 


2 

6 


r 


5 

7 


» Vector from Inertial point to instantaneous center 
(C.M.) of body j j - 1,2,. ..5 

* Vector from inertial point to instantaneous C.M. of 
vehicle system 

a Vector from center of mass of total vehicle system to 
the instantaneous C.M. of body j j = 1,2,... 5 

a Vector from rigid body C.M. of body 1 to C.M. of the 
connection subelement between bodies 1 and 2 contained 
in body 1 when it is in an undeforraed state 

a Vector from the C.M. of the connection subelement 

between bodies 1 and 2 contained in body 2 to the C.M. 
of body 2 when it is in an undeformed state 

a Vector from the C.M. of body 2 to the C.M. of the 

connection subelement between bodies 2 and 3 contained 
in body 2 when in an undeformed state 

= Vector from the C.M. of the connection subelement 
between bodies 2 and 3 contained in body 3 to the 
C.M. of body 3 when it is in an undeformed state 

a Vector from the C.M. of body 2 to the connection sub- 
element between bodies 2 «nd 5 contained in body 2 
when it is in an undeformed state 


from the C.M. of the connection subelement 
^^between bodies 2 and 5 contained in body 5 to the C.M. 
of body 5 when it is in undeformed state 
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Vector from the C.M. of body 2 to the C.M, of the 
connection subelement between bodies 2 and 4 contained 
in body 2 when it is in an undeformed state 

Vector from the C.M. of the connection subelement 
between bodies 2 and 4 contained in body 4 to the C.M. 
of body 4 when it is in its undeformed state 

Linear deformation of body "j" at the rj^ location in 
body "j" referenced to the C.M. of body "j" when in 
an undeformed state 

Vector from the C.M. of body j when in an undeformed 
state to the instantaneous C.M. of body "j" 

Vector from the C.M. of the connection subelement 
between bodies "j" and "1" contained in body "j" to 
the actual connection point between the two bodies 

Angular rate due to deformation of body "j" at 
location 

This rate is expressed in rigid body coordinates of 
body "j" 

Mass of body j 

Mass of subelement of body j 

Transformation from the subelement coordinates 
to rigid body coordinates (i.e,, when body j is 
undeformed) of body j 
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je 


JRN 
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= Transformation from rigid body coordinates of body "1" 
to rigid body coordinates of body "j" 

* External force applied to body "j" 

® External torque applied to body "j" 

* Nx3 matrix of rotational modal gains of body "j" at 
location r^^ 


jTN 

J 

Mn 

j 


Nx3 matrix of translational modal gains of body "j" 

at location r, 
k 

NxN diagonal mass matrix for body "j" 

NxN diagonal damping matrix for body "j" 

NxN diagonal stiffness matrix for body "j" 
modal coordinate for body "j" 


M 

F 

D 


cjl 

cjl 

2j 


= Constraint torque between bodies "j" and "1" 

= Constraint force between bodies "j" and "1" 

= Distance from the instantaneous C.M. of body 2 to the 
instantaneous C.M. of body "j" 


u), = Inertial rate of subelement k 

k 

6^ = Vector from C.M. of subelement "k" to differential 

mass located in subelement k 
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^i(jl) 

r1 
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» Distance from inertial point to the C.M. of the 
connection subelement at the i connection point 
between bodies "j" and "1" contained in body ”j" 

- Distance from inertial point to the C.M. of the 
connection subelement at the i connection point 
between bodies "j” and "1" contained in body "1" 

» Distance from C.M. of the connection subelement con- 
tained in body "j" at the 1^ connection point between 

bodies ”j’* and ”1" to the actual point of connection 

between the two bodies 

® Distance from C.M. of the connection subelement con- 
tained in body "1" at the i'^ connection point between 

bodies "j" and "1” to the actual point of connection 

between the two bodies 

* Distance from the composite C.M. of the connection sub- 
elements contained in bodies "j" and "1" at the i 
connection point to the C.M. of the connection 
subelement contained in body "j" 

= Distance from the composite C.M. of the connection sub- 
elements contained in bodies "j" and "1" at the i^ 
connection point to the C.M. of the connection subelement 
contained in body "1" 

* Linear motion of the C.M. of the connection subelement 
contained in body "j" at the i^ coiinection point 
between bodies "j” and ”1" due to structural deformation 





yJ 


rJ 

i(jl) 


Linear motion of the C.M. of the connection subelement 

th 

contained in body "1" at the i connection point 
between bodies *'j" and "1" due to structural deformation 

Angular rotation of the connection subeleroent contained 

th 

in body ''j" at the i connection point between bodies 
"j" and *’l" due to structural deformation 

Angular rotation of the connection subelement contained 

th 

in body "1" at the i connection point between bodies 
”j" and "1" due to structural deformation 

Location of the C.M. of the connection subelement 
contained in body "j" at the l'' connection point 
between bodies ”j” and "1" with respect to the rigid 
body C.M. of body "j” 



SECTION 9 


CONCLUSIONS AND RECOMMENDATIONS 
INTRODUCTION 

A brief summary of conclusions and re'^ommendations are 
condensed in the following paragraphs. 

MISSION PROFILE 

A mission profile, defining module locations and mass 
properties, was supplied in the RFQ. Operational 
requirements of the Space Construction Base were non- 
existent in this data package and certain assumptions 
were made regarding its operation and performance 
requirements. These many requirements are tabulated 
in section 2 and were used as a baseline fcr subsequent 
studies described in sections 3 through 8. 

ACTUATOR SIZING 

Skylab-type double-gii.ibal control moment gyro (DG CMG) 
units are recommended as the primary actuators. When 
using minimum angular momentum orientations in each 
configuration, the maximum number of DG CMG's is 24 
(see Table 3-12). This is true when the alternate 
"non-rolling" approach is used for solar pointing the 
panels during the X local vertical orientation of con- 
figuration 12. Figure 3-7 indicates that the minimum 
average solar power efficiency would be 80.8%, when B ■- 
36 . 1 degrees . 
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9.4 MOMENTUM DESATURATION 

A momentum management approach is recommended which 
utilizes magnetic torquing as the primary desaturation 
system for configurations 1 through 9M, with an RCS system 
available for backup desaturation during configurations 
8M and 9M (see Table 4-5). 

In configurations lOM and IIM, the RCS system will take 
over as the primary desaturation system, and the magnetic 
torquing system previously available will be usable 
for assistance. This prevents the addition of extremely 
heavy supplementary magnetic torquing coils. 

For the local vertical oriented configuration 12M, a 
fixed tilt gravity gradient desaturation approach 
eliminates the addition of heavier desaturation equipment. 
The RCS and magnetic torquing systems used during 
configurations lOM and IIM will then be available for 
assistance and backup desaturation at this point in the 
mission. 

9.5 CONTROL SYSTEM APPROACH 

Review of control system requirements has pinpointed 
certain features which will require special attention 
during the remainder of the study. These include: 

(1) low sensitivity to model uncertainty; (2) adaptability 
to configuration changes; (3) intermodule statability; and 
(4) high system integrity with respect to component 
failures. The functions of: (1) control coordination 

and (2) intermodule and appendage control require special 
attention due to their uniqueness in this application. 
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study conclusions relative to concept studies are: 

1. Much of the control system for the Space Construction 

Base will consist of elements similar to other 
satellites. However, there are two parts which are 
quite different: module and appendage stability 

control and the vehicle control coordinator. 

2. Coupling control should include two modes: artificial 

stiffening for use during a maneuver and decoupling 
for use in stabilizing a given attitude. 

3. Coupling control should be used selectively based 
on the payoff in performance improvement and the 
impact on actuator requirements. 

4. Rotary torquers may be used for coupling control 
actuators - including locations where CMGs may not 
be used. 

5. Artificial damping of structural modes may be 
difficult to achieve, considering model uncertainty. 

If this feature is required, it may be necessary to 
include adaptive notch or bandpass filters in the 
design in order to achieve significant damping over 
the range of structural dynamics. 

6. Redundant sensors will be required due to the shifting 
of structural modes during buildup. 

7. A vehicle control coordinator will be required to 
select sensors and adjust control parameters as the 
configuration changes. 
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The recommendations are: (1) continuation of this study 

should place emphasis on the vehicle control coordinator 
and the appendage and module stabilizing control with 
particular attention to minimizing sensitivity to modeling 
errors, and (2) there should be a separate effort to 
evaluate simple sensors such as strain gages and optical 
devices for measuring relative angles and angle rates. 

DIGITAL SIMULATION 

Two models are being simulated on a digital computer. 

One is a three dimensional representation of Configuration 
1 of the Space Construction Base. The other is a much 
simpler one dimensional model. Although a hybrid coordinate 
model may be useful with uncontrolled flexible appendages, 
it would not be appropriate for coupling control studies. 

The selected models are satisfactory for the present work. 

MULTILEVEL CONTROL 
Results 

1. The 12 configurations of the Space Construction 
Base mathematical model developed by Bendix Research 
Laboratories have been recast into decomposed state 
variable forms each of which consists of a set of 
decoupled scalar first order ordinary differential 
equations and a set of algebraic equations. 

2. The decomposed state variable form of the mathematical 

model representing each of the configurations of the 
Space Construction Base without control is comprised 
of two sets of equations: translational equations 

and rotational equations. 
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4. The requirement that each of the state variables be 
either observable or reconstructed should be 
Investigated with respect to its effect on required 
computer capacity. 

5. Should the required computer capacity be excessive, 
a new set of models for the 12 configurations of 

the Space Construction Base may have to be developed. 

One of the more promising approaches appears to 
involve the combination of hybrid coordinates with 
multilevel control techniques. 

6. The multilevel control algorithms developed for the 
decomposed state variable model of the Space 
Construction Base should be transformed in such a way 
that they may be applied to the discrete rigid body 
simulation model at Bendix Research Laboratories. 

7. Expansion of the elements of the quadratic coefficient 
matrices appearing in the coordination equations involves 
large numbers of terms. More efficient algebraic 
methods for such expansions need to be developed. 

9.8 STRUCTURAL ANALYSIS OF FLEXIBLE BODY CHARACTERISTICS 

The flexible body characteristics for a series of 
interconnected bodies have been modelled in a free-free 
modal form that permits large angular rotations between 
contiguous bodies to accommodate large deflections during 
deployment and operations. The technique applied provides 
a basis for generation of a more efficient computer 
program with a higher degree of fidelity than those 
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presently being used The unique feature of this 
technique is that the overall characteristics of a large 
space vehicle may be determined by combining the measured 
(or derived) characteristics of the smaller modules 
of which it is comprised. The method developed also 
should allow truncation of substructure modes before 
assembly of the total vehicle model. 

Recommendations concerning the developed mathematical 
model using the above technique are the following: 

a. Investigate the validity of truncation of modes. 

b. Establish criteria by which satisfactory truncations 
can be accomplished. 

c. Validate the model by applying it to a known space 
vehicle such as Skylab or Saturn. 

d. Continue to develop the model in such a way that the 
application to any N-body spacecraft will be simple 
and modular. 

e. Develop a general computer program for the mode]. 
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Subject 

Space Base Mathematical Model 



ABSTRACT 

This memorandum presents a slxty-slx degree of freedom mathe- 
matical model of the space base incorporating all twelve buildup 
configurations. A computer simulation of this model is being de- 
veloped for use in designing and evaluating feasible space base 
control systems. 

1.0 INTRODUCTION 

1.1 Background 

This memorandum is written under the Bendix Research Labo- 
ratories (BRL) tactical program, "Space Base Support," Project 7648, 
Work Order Number 7318, Simulation Model Development. The memorandum 
documents the effort accomplished to date in developing a mathematical 
model of the space base. A computer simulation will be defined from 
this mathematical model to aid in the development and evaluation of 
feasible space base control systems. 

1*2 Obi ectives 

The objective of the Simulation Model Development effort 
is to define the equations for a digital computer simulation. The 
model may be used to examine the stabilization and control of the 
buildup and assembly phase. A specific objective is to structure the 
model so as to allow, by appropriate parameter inputs, any of the 
twelve interim vehicle configurations to be studied. The model must 
be capable of distributing sensors and actuators about the space bane 
and of including short-term external disturbances such as those caused 
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by docking. Vehicle flexibility must also be Incorporated in the 
model. For purposes of studying vernier control system one module 
must be loosely coupled via a six degree of freedom suspension to 
the remainder of the space base. 

1.3 Scope 

This memorandum derives the equations of motion of the 
twelve space base configurations. It is assumed that each configu- 
ration consists of a series of rigid bodies connected by a spring 
hinge suspension. The flexibility of the solar wings are included 
by this method. 

Each Space Bawe configuration is defined in terms of the 
number of rigid bodies assumed, the modules contained in each rigid 
body, and the location of the spring hinges. 

The suspension equations are defined and provisions are 
made for external forces and torques on each body. Actuator forces 
and torques on each body are a''.so provided but the control system 
equations are not included as they are yet to be developed. 

2.0 NOTATION 

The notation used in this memorandum is defined below. 


Symbol Definition 

a^^ Elements on the left side of the equations 

of motion matrix. 


b 


i 


Elements on the right side of the equations 
of motion matrix. 



Spring damping coefficient of the suspension 
torque vector of body j on body i. 
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Definition 

hi 

Actuator force vector on body i. 

hi 

External force vector on body i 

F, 

F - F , + F , 

i 

i ai ei 

\ 

Time derivative of the angular moinentuTn 
vector of body i. 


Inertia tensor of body i. 

hi* hi* hi 

Moments of inertia of body i. 

^XYi’ ^X2i* 

Products of inertia of body i. 

^XM* ^YM’ ^ZM 

Moments of inertia of the space base modules. 

‘'Sij 

Spring stiffness coefficient of the suspension 
torque vector of body i on body j. 

•^Sijx* *^Sijx’ ^Sijy 

Components of spring stiffness coefficient. 

L 

Quantity in the equations of motion that is 
a function of only m^, and 

®i 

Mass of body i. 


Mass of a space base module. 

Zin 

Total mass of space base. 

M 

Quantity in the equation of motion that is 
a function of only and 


Hinge point j. 

**jx' hy* 

Coordinates of the hinge point with re- 

spect to the space base coordinate system. 
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Symbol 

Definition 



R . 

ai 

Position vector from the center of mass of 


1 

^ •> 

body i to the point of application of the 
actuator force. 

1 

1 

R .. 

ei 

Position vector from the center of mass of 
body i to the point of application of the 



i 


! external force. 

« 


• I 




♦ 



Linear acceleration vector of body i with 
respect to inertial space. 

Position vector of hinge point j with tespect 
to the center of mass of body i. 


R • . f ^ * j t ^ j • 

ijx’ ijy* ijz 


Components of the R^^ position vector in the 
body i coordinate system. 



Suspension force vector of body j on body i. 

Actuator torque vector on body i. 

External torque vector on body i. 

Moments resulting from the external and actu- 
ator forces and moments. 


Sij 


SijM 

»j. jj. «i 


Suspension torque vector of body j on body i. 
Torque motor input . 

Coordinate transformation from body i to body j. 

Components of the posicion vector of the com- 
posite center of mass of body i measured with 
respect to the space base coordinate system. 
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Symbol 

Definition 


\ 

3 f/- 

>eiSi 


Component.s of the position vector of the 


i 

} 

f 

: t 

i 

j 


“ij 

center of mass a space base module measured 
with respect to the space base coordinate 
system. 

“i * *ij 

. - f 



'ii 

‘“ij 

X X R^j) 

Angular velocity vector of body j measured 
with respect to body i. 

Components of the vector in the body j 

coordinate system. 

t 



''♦ir '‘®ir ''♦ij 

Euler angles of body j measured with respect 
to body i. 




♦♦ij- “ij- 

Euler angle rates of body j with respect to 
body i. 




4*i» 'I'j 

Euler angles of body 1 me.asured with respect 
to Inertial space. 




♦i- *1 

Euler angle rates of body 1 measured with 
respect to inertial space. 

\ 

IS 


“i 

Angular rate vector of body 1. 

1 



“i 

Angular acceleration vector of body i. 

> 



“ix‘ ‘"iy* “iz 

Components of vector in body i coordinates. 


r 


e • • 

Ui. % (iij 9 u. 

ix iy Iz 

Components of vector in body i coordinates. 
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Synfeol 

BFMi 

CCM 

DM! 

ETl 

HMi 

KWA 

KWB 

Kv:s 

LMi 

MPi 

MSM 

Pi 

PSP 

SMi 

SPM 

SWA 

SWB 

TAi 


Definition 

Beam Fabrication Module 1. 
Construction Control Module 1. 
Docking Module 1- 
External Tank 1. 

Habitability Module i. 

150 KW Solar Wing A 
150 KW Solar Wing B. 

150 KW Subsystem. 

Logistics Module i. 

Manipulator 1. 

Mission Support Module. 

Pallet i. 

Public Service Plate. 
Subsystem Module i. 

Space Processing Module. 

Solar Wing A. 

Solar Wing B. 

Turret Assembly i. 
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3.0 DERIVATION OF EQUATIONS 

3.1 General 

The Space Base is assumed to be divided up in a 
group of rigid bodies connected by spring hinge suspensions. 

It is assembled, starting with the baseline Configuration 1, by 
adding on additional modules. As other modules are added, a 
total of twelve configurations are defined. This section defines 
for each configuration the number of rigid bodies assumed in model 
derivation, the modules comprising each rigid body, and the location 
of the spring hinges. 

The degrees of freedom are also defined in this section 
along with derivation of the equations of motion. Coordinate trans- 
formations relating vector components to the various body coordinate 
systems are defined. Euler angle rate and suspension equations are 
also presented. 

A block diagram of the overall mathematical model is 
shown in Figure 1. 

3.2 Description of the Configurations 

The Space Base buildup consists of i. series of twelve 
configurations. Each configuration is divided up into a group of 
rigid bodies. Table 1 lists the nurber of rigivi bodies assumed for 
each configuration. 

Each rigid body is a collection of modules or part of a 
module. For example, solar wing A is divided into two rigid bodies. 
Table 2 defines the breakdown of the rigid bodies and modules for 
each configuration. 
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Table 1 


Configuration Number of Rigid Bodies 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


5 

6 

7 

8 
8 
8 
8 

10 

10 

19 

20 
20 


gure 2 presents a topological tree that shows how the 
rigid bodies are connected. A single line drawn between bodies indi- 
cates a three degree of freedom spring hinge suspension. A double 
line indicates a six degree of freedom suspension. 

figure 3 shows a sketch of the space base and locates all 
of the spring hinges. 

3 . 3 Degrees of Freedom 

As shown in Figure 2, a three degree of freedom hinge is 
assumed between each body except bodies 1 and 8. These two bodies are 
assumed to be connected by a six degree of freedom suspension. As a 
result, the total degrees of freedom for each configuration are defined 
as shown in Table 3. 
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Table 3 - Degrees of Freef! 


Configuration 

Degrees of Freedom 

Translational 

Rotational 

Total 

1 

3 

15 

18 

2 

3 

18 

21 

3 

3 

21 

24 

4 

6 

24 

30 

5 

6 

24 

30 

6 

6 

24 

30 

7 

6 

24 

30 

8 

6 

30 

36 

9 

6 

30 

36 

10 

6 

57 

63 

11 

6 

60 

66 

12 

6 

60 

66 


3.4 Equations of Motion 

3.4.1 Translational Equations 

The translational equations of motion for each 
rigid body is derived by equating the mass m^ times the linear accel- 
eration vector to the sum of the force vectors (suspension forces 

S. . and external F , anu actuator F . forces; F, “ F . + F .) acting 
ij ei ai i ei ai 

on the body. Hence, 


“l *^1 ' "^12 ~ ®13 ■ ®16 ■ ^18 " ^1,11 ^1 


( 1 ) 


“2 *^2 ‘ hi ' ^24 ^2 


( 2 ) 
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^ :• i 





^^ternal 

'^s'worsndum 

K. 1977 , 

'■«8« 14 ^•^No. 


sai 


®3 Sj *• c 

13 * S^3 + 

*4 

®4 R, • c , 

4 ^24 + P, 


®5 fir • e 
^ ^ ^35 + ^5 

*• 

\ K • s 

“ ^67 ■'■ K 


-;EPROPUCIBILTry #F TOE 

PAGE*« FfW 


“7 *7 • S 


67**'^ 


79 ^ s 


* 8 hs^ F, 




8 


’‘^ ■ ^79 F„ 


^0 « 


10 * s, , 


7.10 + F. 


10 


“11 \i * S. 


1.11 ' s. 


11.12 ' S 


■^2 R 


11.13 F 


l2 • S 


11 


I1.12 ' s, 


"13 ^3 • S 


12,14 F. 


12 


11.13 ' S 


"W «14 - S, 


13,15 F. 


13 


12,14 “ S, 


"^5 


l-i.16 F. 


14 


23,15 “ S 


15,17 F. 


”ie ^ 


15 


16 • 5 


14.16 ‘ S 


“^7^7 


16.18 F. 


16 


15.17 - S 


17.19 F, 


17 


(3) 
(4) 
f5) 
C6) 

C7) 

f8) 

f9) 

flO) 

fll) 

fl2) 

(13) 

(14) 

(15) 

(16) 

(17) 
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“l8 ^8 “ ^16,18 ^18 

,s 

“l9 ^19 * ^17,19 ^19 

“20 ^20 “ ®7,20 ^20 


( 18 ) 

(19) 

( 20 ) 


The linear acceleration of body 2 (R 2 ) can be expressed 
in terms of the linear acceleration of body 1 (R^) and the angular 

a e 

accelerations of body 1 (lo^) and body 2 ( 102 ), i.e.. 


R2 “ X R^2 “ «2 * ^22 “1 * ^“1 * ^ 12 ^ “ i*>2 X (^2 x R22) 


( 21 ) 


where R. . is the vector from the center of mass of body i to the 
th 

j hinge point. The linear acceleration of the other bodies can be 
expressed in the same manner. 

Summing equations (1) through (20) excluding (8) and sub- 
stituting equations like ('2 1 . yields an expression for Rj^, i.e.. 


where 


% “ ^ - ^18 - N - L] 


20 

TM - ”i - ”8 

i-1 


20 


EF 


Eh- 


i-1 


( 22 ) 


(23) 


(24) 
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a * (i) X R 

ij 1 ij 


hi ■ “j * '“i ’ 'ij> 


Then L and N are functions of j ’ ^ij* 


L = f(m^, 


N «= f(m^, 


As a result, the linear acceleration can be vnritten 

a 

in terms of m, , R, F F S, u) . , and lo . ~ The quantities m. 

1’ Ij ei ai 18 i i i 

and are given constants, F^^^ is a given external input, F^^ is 
control system output, is a suspension force output, is the 

a a 

integral of and finally u)^ is determined in the following section 
which derives the rotational equations of motiona 


The suspension forces are needed for the rotational 
equations of motion. They are given by 


^ij “ ^^“i’ ^i* “ij* ®ij^ 


The linear acceleration of body 8 is given v 


's ■ ± '*!« * "s' 
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3.4.2 Rotational Equations 

The rotational equations of motion for each 
rigid body Is derived by equating the rate of change of angular 
momentum to the sum cf the moment vectors acting on each body. 

e 

The rate of change of angular momentum h^ Is defined by 


hj - “i + “2 * 'i "i 


(31) 


where Is the Inertia tensor 


/ 


I 

-T 

-I,.- \ 

Xi 

XYl 

XZl ^ 


^Yi 

■^zi 

-I 

-1 

i 

XZl 

YZi 

Z1 / 


(32) 


Then for each body 


* ^12 ■ ^13 * ^13 " ^,11 ®1,11 " ^16 * ^16 

' ^18 * ^18 ' '^812 ■ "^SIS " "^Sie " "^Sl,!! " ^818 ^1 


h^ - R 22 X Sj^2 " ^‘24 ^ ®24 ^S12 " ^824 "^2 


(34) 


hj « R^3 X 833 - R35 X 835 + Tg33 - Tg33 + T3 


(35) 


^4 ■ R44 X S24 + Tg 24 + T^ 


(36) 


hj - R53 X 833 + Tg33 + T3 


(37) 
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^6 “ *^66 ®16 ■ *^67 * ^67 ■*■ ’‘$16 “ '^867 *^6 


^7 “ *^77 ^67 ■ ^79 * ^79 ‘ ^7,10 * ®7,10 " ^7,20 * ^7,20 

■'■ ^S67 ” "^879 ‘ ^87,10 " ^87,20 ■*■ ^7 


^8 “ ^88 ^18 "*■ ^818 ^8 


hg * R99 X Sy^ + Tgyg + Tg 


^10 ' *^10, 10 * ^7.10 '*' '^87,10 ’’10 


^11 “ ^11,11 ^ ^1,11 ■ ^11,12 * ^11,12 ■ ^11,13 * ^11,13 
■*■ ^81,11 ■ ’’811,12 ■ '^811,13 ■*■ "^ll 


*'12 “ *^12, 12 ^ ^11,12 ' ^12,14 ^12,14 ’'sil,12 

■ ’^812,14 ^^12 


^13 * ^13,13 * ^11,13 ‘ *^13,15 * ^13,15 "*■ ’’811,13 
" ’^813,15 ’^13 


^14 " *^14,14 * ^12,14 “ ^4,16 * ^14,16 ’’812,] 


- T + T 

^814, 16 ^14 


A-18 



Data 


Internal 

Memorandum 




December 1977 Letter No. 
Page 19 


''is “ ''iS.l'i ^ ^13,15 ■ *'l5,17 * ®15,17 ‘''$13,15 



“ ‘'^S1.S,17 

•^^15 


(47) 

^16 

“ ''l6.16 

^14,16 

“ ®16,18 ^ ^16,18 ^$14,16 



■ "^$16,18 



(48) 

^17 

* *'l7,17 * 

^15,17 

■ ''l7,19 ^ ^17,19 ^$15,17 



" ‘'^$17,19 

^^17 


(49) 

e 

h, „ 

* R ^ ^ X 

s , , 

+ T • - ^ T 

(50) 

18 

18,18 

16,18 

sib, 18 U8 

e 

h,„ 

• R- ^ « X 

S. ^ 

+ T , , . + T. 

(51) 

19 

19,19 

17,19 

S17,19 19 

*^20 

■ *20.20 * 

^7,20 ' 

► T + T 

'$7,20 '20 

(52) 


where are moments resulting from the exterual and actuator forces 
and moments, i.e.. 


T. “R. xF/+R. xF.+T.+T. 
1 ei el al al al el 


(S3) 


and where T_.. is the spring hinge torque vector of body j on body i. 
oij 
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These equations can be manipulated into a scalar 
matrix of the following form; 


®11 ®12 


®21 ®22 


*57,57_| [_“20z 
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where the eleisents are functions of R 

^XYl’ ^XZI’ ^Zi* '*’i’ ®i’ quantities 0^, and 

are the Euler angles of body i. 


ij* ^xi’ hi* hi* 


The elements are functions of the same 

quantities as a^^ along with S^^g, , 

R . , and R . . 
ai* ei 

Inverting the matrix provides the solution for 

e e e 

the components of angular>acceleration *^iy* *^iz^ 

bodies except body 8. The equation for the angular acceleration of 

body 8 is given by 


^8 “8 


*^88 * ^18 ^ S 18 ”*■ ^8 * “8 * ^8 “8 


(55) 


3.5 Transformation of Coordinates 

In order to expand the vector equations of motion, derived 
in the previous section, into scaler equations all the vectors in an 
equation must be expressed in a common coordinate system. As a result, 
provisions must be made to transform vector components from one body 
coordinate system to another. This is accomplished using the notation 

Vector 
Components 
in dody j 
Coordinates 

The transformations ^Tj, jTj, gTj, gT^, ^Tg, ^T^, 

10^7* 20*^7* 14^12’ 15^13’ 16^14* 17^15’ 18^16* 19^17 

angle transformations. Hence, they are defined by 


- T 

j 1 

Vector 
Components 
in Body 1 


Coordinates 


(56) 
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i'.l'l.'Ni 


jl,- 


i .. 






'“li> 


h 

0 

V 

1 1 

0 

f 

I 1 

-Aij) 


0 

1 

-A 4 


ij 


ij 


0 

1 

0 


Li, 


ij 
1 / 

0 

1 / 




ij 


\ 0 


1 

0 


and where Euler angles 

respect to body i. Expanding 


L ** 

j i 


1 ■ 



-‘♦iJ 

1 

^♦ij 


-A6^. 

1 

i 


ml^ 




■'\^3rLmr of 

,L PAGE TS POOR 

(57) 


( 58 ) 


( 59 ) 


( 60 ) 


of body j with 
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For bodies having a large angle articulation a different 
Lon is require 
and are defined by 


transformation is required. Hence the transformations 2 Tj^, 


T 

j i 


cos A\L(^j 

sin Aiii^j 

-A6 

-sin Au/ . . 
ij 

cos Aiji . . 
ij 

A(ji 

cos Alii . . + A4i . . sin Ail . . 
ij n ij 

A6^j sin Aiji^j - Aiji^^ cos Aii-^^ 

1 


(62) 


\ 


I 


All other transformations car. be derived from those defined 

above. 

3.6 Euler Angle Rate Equations 

le Euler angle rates . ^®ij’ '.-■elated to the 

body angular rates by the vector equation 


- (A^^j)(A6^^) aV^^ + (A^^j) A*e^^ + A>^j 
Expanding into scalars 


- ~i 

Au. . 
ijx 



0 



1 

o 

Au . . 

ijy 

Au. . 

ijz 


0 

1 

W 


0 

"♦ij 


t . 1 

- 


— 


o 

o 


0 


^♦ij 

-f 

0 1 4,.. 



+ 

0 




ij 




V » / 


0 


0 

' 1 

«. — 


__ ^ 


(63) 
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Solving for the Euler angle rates yield 
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where 


441 ., •= All).. + A6.. All),, 
ij ijx ij ijz 


A0 . , * All).. - A4 . . All),, 
ij ijy ijz 


A\^. . « A4>. . Aw. . + Aw. . 

ij ij ijy 1J2 


‘“ij • “] - 3^1 "I 


For body 1 


^ = ID- + e ID. 

1 lx 1 Iz 


6- = ID, - 4', u, 

1 ly 1 Iz 


4», * 4 , ID, + ID, 

^1 ^1 ly Iz 


3.7 Suspension Equations 


The suspension torques T-., generated by each of the spring 

Olj 

hinges are assumed to be linear functions of the Euler angles > 

^®ij’ rates (40)^^^^^, 4 id^j^), i.e.. 


T_ . . ■ K- . . A4j . + C-.. Aid.. 

Sijx Sijx ij Sijx ijx 
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T 


T 


Sijy “ *^Sljy ^®ij Sijy ^“ijy 

(73) 

Sijz * ^Sijz ^"^ij Sijz ^“ijz 

(74) 


where ^Sij spring and damping constants. 

For the bodies that have large articulation angles 
about the z axis, the suspension torque is given by a torque motor, 
l.e. , 


T 


Sijz 


T 

Sijm 


(75) 


where is a torque moto’- input. 

3.8 Equations for Computing the Composite 
3.8.1 Composite Masses 

The composite masses of the rigid bodies are com- 
puted simply by summing the masses m^^ of the modules comprising the 
rigid body, i.e. , 


3.8.2 

as follows: 


“ I m^ (76) 

Inertia Tensor 

The elements of the inertia tensor are computed 


*Xl ■ ^ 'XM * ^ 




“m (^M - ^i)^ (*M ‘ *i)^] 

^ ^ “ *i) ’*■ (*M ■ *i) J 


(77) 

(78) 
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"m [(’Si - *i)^ * (>’m - n) J 

(79) 

■Si <’Si- ’‘1?%- >’i> 

(80) 

■V '“m- V 

(81) 

^ (y„ - - ».) 

(82) 


where Zj^ are the coordinates of the center of mass of module M 

measured in the space base coordinate system and x^, y^, are the 
coordinates of the composite center of mass of body i measured in the 
space base coordinates. 

3-8.3 Hinge Points 

The hinge point components ^ijy» locate 

the hinge points P. with respect to the center of mass m^ of the com- 
posite rigid body. The components are computed by 


R. . = P. - X. 

11X ]x i 


(83) 


R. . " P. “ y . 

ijy jy i 


(84) 


ijz jz i 


(85) 


3.8.4 Composite Center of Mass 

The location of the composite center of mass of 
each body must be computed. The components y^, are c*7iputed by: 
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SUMMARY 

This interim report describes some preliminary results concerning the 
application of the Multivariable Nyquist Array method to the design of an 
attitude control system for a flexible space vehicle. The system used in 
the study is the lightly damped, three body model developed by Porcelli. 

With the exception of the open loop model description, the parameters 
resulting from the MNA design were obtained independent of the Porcelli 
control configuration. For design comparative purposes, an alternate single 
loop elimination design was made using root locus and Bode methods. 

An initial application of the MNA program to the full order model of 
the space vehicle suggests a decoupling of the control objectives- This 
result was confirmed by a CSMP study of the open loop dynamics as well as by 
physical considerations. Using this result, a reduced order model was 
developed for the design of the control units for the flexible appendages. 
Application of the MNA program to the reduced order model clearly indicates 
a need for phase lead compensation in each control loop. A CSMP simulation 
of the full order model with the decoupled control system design demonstrates 
the utility of the MNA design method for flexible space vehicle systems* 
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Section 1. INTRODUCTION 


A recent paper by Porcelli [l] details a multi^loop attitude control 
system design for a flexible space vehicle using a three body model repre-* 
sentation. Using the con rpt of mode separability, Porcelli shows that when 
some of the lightly damped closed loop dynamic modes are contained within 
the control system bandwidth^ sustained oscillations may occur in the flexible 
appendages during cransient operations. To iirprove the dynamic response, a 
pair of auxiliary control loops are established to increase the structural 
damping by active means. 

Since the main control loop is designed in accordance with the attitude 
control specifications, the auxiliary control loops may be designed independ- 
ently- The design approach proposed by Porcelli initially ignores the flex- 
ible connections between the vehicle bodies. Each auxiliary loop design is 
then based upon the rigid body dynamics of the component body to which it is 
directly related. Recognizing the conservative aspects of this design 
approach, the auxiliary loop bandwidth is progressively decreased until an 
acceptable design is obtained. For the examples in [l] , auxiliary loop 
bandwidths one decade below the bandwidth of the main control loop proved 
satisfactory. 

The Porcelli method is useful when all dynrjnic inodes attributed to the 
flexible appendages are either within or exterior to the control bandwidth- 
The examples in [l] demonstrate both cases. Situations may arise, however, 
where some of the modes associated with the appendages are within the main 
control bandwidth while other modes are exterior, ia this situation^ the 
design approach may not be quite as clear as in the former cases. Here it 
may be desira. le to have an appendage control unit which utilizes sensor 
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information from each body simultaneously. This is particularly true if 
physical considerations limit the position and number of appendage control 
units. Further considerations regarding N-body models may prove to be 
untractadsle using the Porcelli method. This would certainly be the case, 
for exeunple, if the structural modes were highly interactive. Here the 
flexible members could not be ignored in favor of a rigid body analysi*! emd 
design, as required by the Porcelli method. 

As an alternative to the Porcelli method, a Multivarieible Nyquist Array 
(MNA) design was initiated for the three body models in [l] . One distinc- 
tive feature of the MNA method is the utilization of the complete system 
model at each stage of the design procedure. Any simplifying ass^'jnptions in 
the model order or system characteristics, if any, are based exclusively on 
physical considerations. Thus, in terms of the satellite model, all flexible 
appendage modes are retained and accounted for. 

The next section develops the mathematical model for the three body 
system. Section 3 briefly reviews the open* loop dynamic conditions from a 
root locus viewpoint and outlines an alternative classical control design 
procedure by successive loop reductions. The CSMP runs using the Porcelli 
control configuration are also presented and compared with the alternate 
designs. Section 4 describes the MNA design for the lightly damped case. 
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section 2. MATHEMATICAL MODEL 


Th« system model for the flexible space vehicle is presented in Figure 1 

e 

where all system components represented are considered to be ideal elements. 



Figure 1: Full Order Model for Attitude Control 

The system inputs are represented by torque drivers 't^(t), Tg(t;» and'Tj^jj(t) 
with 6^(t)> 62(t), and 8^(t) selected as the measurable set of system outputs. 
Using the systems graph approach, Figure 2, the following set of circuit and 
cutset equations are obtained: 



Figure 2: System Graph 
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Circuit Equations 

<;- » - 0 
CiJ' 0 

Cf + ^ - 0 

^10 ^3 


vr 






Cutset Equations 

'f 1 * ”^4 * - 'r. ■ ■> . 

"^3 ” "^5 “ “ "^10 “ ° 


( 1 ) 


The component equations corresponding to Figure 1 are: 


j. l!i . r 

^ dt ^ 


. 

5T - ^®i 

T - B.e. 

1 11 

T! - drivers 
1 


i - 1, 2, 3. 

i - 4, 5 
i = 6, 7 
i « 8, 9, 10 


( 2 ) 

(3) 

(4) 

(5) 


Selecting the state variables as 6^^, 6^, 6^, T^, and 6^^, the following 
state equations are obtained 


d6 'T' 
1 ^ ^ 

dt * 

li.Tk 

dt “ J 

^.h 

dt “ J 


[■^8 -fi-V -57 It. - t 

~ rt * t - n * »6'«i - 

j; [ ^0 * n * V ’ 37 tto * t * >7 - » 3 >] 


(6) 
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i!k 

dt 

!!k 

dt 

dt 


>^4 (6^ - e^) 


Kj (02 . 63) 


with output equations 




8 - e - — 

4 

'^4 '^S 

0 .0 - i 

3 1 Kj, 


In vector form the above equations become 
X = Ax + Bu 
y » Cx 


(7) 


( 8 ) 

(9) 


where 


i -(=?) 


0 

1 


-K, 


0 

-K, 


- ^ 0 


0 

0 

0 


0 

0 

0 


1 


“2 

4 0 


0 

0 

0 


1 10 ) 
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9y - -.0001 t i *01414 

• **.00035 + j .0265 

Th# periods of the lightly damped mod^s acre, respectively, 

■ .01414 ■ 444.36 secs 

• .0265 Tj » 240.46 secs 


• ( 15 ) 


(16) 


In the next section, the results of Porcelli are compared with an 


alternate synthesis appro«»ch using single loop elimination. 







Section 3. CLASSICAL DESIGN APPROACH 

For the state variable equations of Section 2, the open loop transfer 
functions of Ap^andix A were obtained from the application of Danielevsky's 
method ^2j. It is interesting to examine the transfer function relating the 
angular displacement of mass with the input to mass 1 in closer detail. 
This transfer function is repeated here as 


(s) 


T;(s) 

8 




.00083^ + .00085^ + . 16E-6S + .8E-7 
S^(S^ + .00095^ + .00095^ + .28E-6S + .14E-6) 


The poles and zeros of obtained as 


Zeros 


Poles 


2^, 2^ = -.0000586 > j .01082 


2^0 2. = -.00034 + j .0261 
3 4 


*= -.0001 ^ j •01414 
= -.00035 + j .0265 


p^ ® 0.0 


p * 0.0 
6 


With this pole-zero arrangement it is apparent that the complex zeros 
virtually cancel the con^lex poles and the transfer function could be 
reasonably approximated by 


(17) 


(18) 


6i(s) ^ 

TgCs) ' ""ll^^* = 32 


(19) 


Figure 3 contains a sketch of the root locus for the uncompensated 
system. Figure 4 provides the Nyquist diagram for the full order model. 
The large loops in Figure 4 are due to the resonant frequencies associated 
with the appendages. Note that these U quist loops do not effect the 
stability issue for the overall system. 
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Clearly to stabilize the system for high gain and to obtain the closed loop 
dynamic response desired/ it is necessary to inject phase lead compensation. 
For all practical design purposes the control design from the terminals 
can proceed using the approximation in (19). 

Following standard phase lead design procedures# a Bode plot of 


G ( s ) * 2 

s 


is made and appears in Figure 5a. 


( 20 ) 


The phase lead compensator has the form 

-“(rfUs) 

Selecting Q = .1 to inject a large amount of phase lead at the gain cross- 
over point 

Sin 0 = 7 ^ or $ = 54.9° (22) 

m 1 + a m 

Since the magnitude curve is also adjusted when the phase is modified# the 
new crossover frequency is obtained from th§ -20 log^^(l//a) point on the 


magnitude curve # i . e • # 


-20 log^^d/Ta) = -10 db 


(i) 


m 


1.779 


The zero and pole location for the compensator are obtained from 

i = /a 0) = . 562 

T m 


(23) 

(24) 


(25) 


= 5.62 
aT 


Thus 


or 


> if 1 1>779S 1 

» 1 [ ^ .1779S J 


(27) 
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G (s) 
c 


S + .562 
S + 5.62 


The compensated system is indicated in Figures 5a and 5b. Figure 6 


contains the Nyquist diagram for the compensated G(s) in (17). 
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T a 10.9 seconds 


(32) 


with similar oscillations occurring in the appendages. In each case,* an 
initial offset of .02 radians was assumed for all vehicle bodies. 

Based upon the above analysis, it can be concluded that the main control 
xinic for the satellite is only concerned with the dynamics of the large mass. 
This is further supported by the CSMP runs wherein the angul^ displacements 
associated with the appendage masses at t » 3.9 seconds remain within 1.3% 

and .01% of their initial offset:*, respectively (Appendix C) . Thus for the 
design of the appendage control units, it can be ’ssumed that the main body 
has been returned to the equilibrium point (the origin). 

Ignoring the main body is equivalently reflected in Figure 7. 



Figure 7: Reduced Order Model 

Following a development similar to the full order model (FOM) the state 
equations for the reduced order model (ROM) are obtained as 

X * Ax ♦ Bu (33) 

y « Cx (34) 
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where 



The transfer functions for the reduced order model appear in Appendix B. An 

examination of th* transfer function relating with ^ 

reveals the pole-zero configuration 

Zeros Poles 

*1' *2 * i .019998 p^, p^ ■ -.0000586 + j .010824 

(39) 

^3' ^4 * "•00'^3414 + j .02613 


B-16 



7 ^% 


«s indicated in Figure 8 with the Sode for the ROM in Figure 9 




Figure 8: Root Locus for Reduced Order Model 




For this model, pole-zero cancellations may once again be assumed with the 
result 


e^(s) 


^ (si 2 -4 -4 

* 9 ' ' S + 1.172(10 )S + 1.1716(10 ) 

With a one percent steady state error specification 


lim 1 

s -► 0 1 + G(s) 


.01 


from which the open loop gain is obtained as 

K = 1.16(10*^) 


(40) 


(41) 


(42) 


Based upon the transfer relationship in (40) it would appear that a 
bridged-T-compensator would be significantly more effective than the simple 
lead compensator used by Porcelli. For the coefficients in (4v^/, the 
bridged-T^compensator has the form 


, . * 1.172(lo“‘ *)S 1.172(10*^) 

c'®' * S(S ♦ .9888) 


(43) 
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which yields the open loop transfer function 

G (s)G(s) ■ ^ (44) 

S(S + .9888) 

The gain K in (44) could be chosen in correspondence with a steady state \ 

error specification to a ramp input or to obtain a desired set of do 
loop poles with specified damping characteristics. To modify the dynamic 
characteristics further, a lead-lag network could be cascaded with the 
bridged-T network . 

At this point in the design, a new state model would be developed for 
the space vehicle with both control loops closed. Using Danielevski * s 
method, the system transfer function from could be obtained. 

The feedback design could then be completed in a manner similar to the design 
procedures above. CSMP runs could then be used to evaluate the dynamic 
response and adjust the appropriate compensator to correspond to the practical 
requirements of the vehicle. 

From an examination of the transfer function for 6 ^(s)/T^q(s> in 
Appendix B, it is clear that the analysis above would apply directly to the 
design of the feedback control for the third body. Experience suggests, 
however, that feedback control designs obtained for each 2 op independently 
may have a deleterious effect on the dynamic behavior when all loops are 
closed simultaneously- This condition is due primarily to the effects of 
system interaction among the control loops which were not incorporated into 

the above designs. To obtain a measure of the interaction levels for the 

. * w 

control designs ^^op<^ed above as well as those due to Porcelli, Davison‘s 



The concerns cited above for control system designs obtained by ignoring 
the interaction effects are relieved when the Multivariable Nyquist Array 
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method is employed since interaction levels are used ii. the design procedure. 
The next sections present a preliminary design for the flexible space vehicle 

e 

using ths KNA dssl9n Mthod. 




Qp 'TttB 
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Section 4. MULTIVARIABLE WYQUIST ARKAY METHOD 




The fundan^ental objective of the design method is to decrease cross- 
coupled system interaction to such an extent that the closed loop system 
design reduces to a set of independent single loop design problems. Although 
simply stated, the actual reduction procedure proposed by Rosenbrock [ 4 ] 
and implemented by Munro [s] requires a high degree of designer intervention 
and is fundamentally a trial and err^x process. 

In Figure 10, G(s) is an mxm transfer matrix representing the coupling 
of m inputs and m outputs. 



Figure 10: Mua civariable System Configuration 

The pre- and post-compensator i/.atrices K -ind L, respectively, are each of 
dimension mxm. TJie feedback gain matrix, F, is assumed to be diagonal and 
of similar dimensions. Clearly, if 

Q(S) * LG(s)K (45) 

is diagonal, loop closure may proceed on an :ndividual loop basis with a 
guarantee of zero Icop interaction. It is this premise upon which the MNA 
design philosophy is based. The adherence to strict diagonalization is 
relaxed, however, w* th the substitution and exploitation of the concept of 
diagonal dominant matrices. 
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Definition 


A matrix Z(s) is diagonal dominant if either or both of the following 
conditions srs fvs3«nt for all s: 
m 

a. ][ | 2 ..(s}|/|z..(s)| <.6. <1 for all isl,2« . . . ,m (46) 

j-1 ^ ^ 

i/j 

m 

b. I |/| 2 ^^(s) I < < 1 for all i«l,2,...,m (47) 

j*l 

j/i 

Equation (46) defines row dominance while (47) defines the column dominance 
condition where 6^ is the level of dominance obtained tor the ith diagonal 
element. 

Before the design process can proceed further^ Q(s) for the DNA method 
or Q ^(s) in the INA method must be made dominant by manipulation of the 
elements of the compensator matrices. Once dominance is achieved the design 
process is completed using single loop theory to select the diagonal elements 
of F. This selection process is enhanced tfirough application of the Gershgorin 
and the Ostrowski theorems for dominant matrices. 

The Gershgorin theorem [4j states that the eigenvalues of a matrix 
(cither Q(s) or Q ^(s)) are located in the union of the bands centered about 
the diagonal elements with widths determined by the sum of the moduli of the 
off diagonal elements by row or by column. Using the envelope procedure 
developed by Crossley [js] and considering each control loop separately, a 
graphical display of open loop system interaction results. Figures 11 and 12 
indicate a typical display for an INA and DNA design formats respectively. 

The feedback gain selection for control loop i is then made in correspondence 
with the generalized Nyquist criterion amd the stability theorems of 
Rosenbrock £ 4 ] . 
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The Ostrowski theorem may be used to further shrink the Gershgorin 
bands# thereby reducing the area of uncertainty in each loop. This s^et of 
bands is frequently referred to as a set of **fuzzy*' Nyquist plots (or 
inverse Myquist plots for the INA) • Using the innermost band as a conser* 
vative estimate of the Nyquist contour in each loop# the design proceeds on 
a single loop basis. Feedback gain selection must be made exterior to the 
Gershgorin (Ostrowski) band. Thus# phase margin# gain margin and dynamic 
compensation may be used to evaluate and/or improve the loop design with a 
guarantee of low interaction from the closure of the remaining loops. 

Diagonal dominance for the Direct Nyquist Array method requires the 
selection of pre- and post-compensator matrix paralneters so that (46) is 
satisfied when (45) represents the open loop transfer matrix. For the 
Inverse Nyquist Array method Q ^(s) is used and the parameters of K ^ and 
L ^ must be selected. 

An efficient and reliable method for the evaluation of the matrix 
coefficients is described in [^7 j . The dominance algorithm uses a conjugate 
direction function minimization algorithm to adjust the parameter set until 
a performance index composed of the dominance definitions in (46) and (47) 
is minimized. For the INA method in a row dominance mode# the optimization 
problem can be separated into three independent optimization efforts; one 
for each row. Here the p>erformance index by row is 

m ^ 

J (K ) - Max I lqi.(s)l/lq. .(s) I (48) 

J w j=l J 

if'j 

A -2 

where qj^j(s) is an element of Q (s). For each i, the ith row of K is 
adjusted until minimized. In practice* the ratio in (48) is 

computed for each discrete frequency point in the range of interest. This 
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array is then scanned to identify the maximum ratio. Adjusting the elements 

of row i in X ^ yields a set of final dominance levels 

e. ■ Min J. (K, .) (49) 

1 ^ i ij 

If the dominance levels in (49) are less than unity# diagonal dominance has 
been achieved. In the event that some of the dominance levels are greater 
than one# the designer may initiate a dominance sharing search or restart 
the program using new starting values for tha unspecified compensator 
parameters. 

The concept of dominance sharing is detailed in [sj • It is funda- 
mentally a rescaling of the compensator matrices to the extent that low 
dominance levels may be intentionally increased to a point where the previous 
non-dominant levels may be shifted to a range of acceptability. This proce- 
dure has been automated in the latest version of the dominance algorithm 
and is initiated by the designer after the final set of dominance levels have 
been evaluated. » 

Once a set of coefficients have been determined for dominance of the 
open loop transfer matrix# each control loop may be treatec independently 
using single-input single-output control theory. The dominance algorithm 
briefly outlined above shifts the burden of esta±)lishing the dominance 
condition from the designer to an automated procedure. Thus designer inter- 
vention is only required during the actual design process and is no longer 
needed to establish the requisite dominance condition. Experience with the 
algorithm suggests that dominance may often be obtained within several CPU 
minutes or less depending upon the characteristics of G(s)* 
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Section 5: MNA DESIGN 


The Multivariible Nyquist Array method described in Section 4 was 
applied to the full order model of the space vehicle. Based upon an 
evaluation of these results ^ the reduced order model was used subsequent to 
the design of a control xinit for the large mass. Although the design 
described herein is of a preliminary evaluation for space vehicles, the 
results clearly demonstrate the utility of the MNA method to flexible vehicle 
control system design. 

Using the transfer matrix of Appendix A as the frequency domain repre- 
sentative for the full order satellite model, a DNA design run was initiated. 
For this DNA design the postcoi..pensator matrix was prespecified as the 
identity matrix with the precompensator matrix to be selected in accordance 
with the algorithm in [ 7 ] . Figure 11 indicates the closed loop diagram for 
the satellite with 

Q{s) - LG(s)K (50) 

representing the open loop transfer matrix. 
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Figure ll: Closed Loop Satellite Control Design 




Direct application of the DMA method to the FOM without regard to the 
ain^lifying conditions described in Section 3 produced a non-dominant condition for 
Q(s) in (50) . Using a frequency range of 

.005 < 0) < .03 

for dominance evaluations, it is evident from the computer printout that 
dominance for each column of Q(s) is lost in the frequency range near the 
resonant frequencies associated with the flexible appendages. Since this 
design effort is using the DMA column objective, loss of dominance is not 
particularly significant providing that feedback gains could be selected 
external to the corresponding Gershgorin bands. This condition, however, is 
non-existent for loop 1 as the Nyquist band is coincident with the negative 
real axis. This condition immediately relates the dominating influence of 
the two open loop poles of Q(s) at the origin (i.e. pnase angle of 180°). 

Thus a confirmation of the single loop analysis from Section 3 is obtained 
directly from the DNA diagram for column 1 of Q(s) (see Figure 4 ). With a 
Nyquist diagram of this form, phase lead compensation is required to provide 
an adequat .■ phase margin for the main control loop. 

Using the phase compensator designed in Section 3, the angular displace- 
ment of mass 3 ^, is returned to the equilibrium point (origin) before any 
significant motion of the appendages takes place. Hence the model configura- 
tion can be modified to the reduced order model previously developed. 

App' cation of the MNA program to the ROM over the frequency interval 
0 £ (I) £ .03 radians (51) 

''ields a dominant condition for each column of the ROM transfer matrix Q(s). 

Here the matrix G(s) in Appendix B provides the frequency domain description 
when mass is clamped. 
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The dominance producing compensators are 


r 1.0034 

-.94747 

[-.70949 

.66997 

r~ 


j 1.7174 

0 

L 0 

1.2128 


(52) 


(53) 


Figures 12 and 13 display the direct Nyquist diagrams for columns one and two 
respectively. The envelope of the Gershgorin bands are indicated in each 
figure. 

Since the feedback gains must be selected exterior to the Gershgorin 
bands f it is clear from the DNA diagrams that some form of phase lead 
compensation is required. This form of compensation will effectively swing 
the DNA bands down and away from the negative real axis, thus providing an 
improved phase and gain margin design. This procedure will then allow for 
an increase in the system gain space for each control loop. 

Alternatively, a bridged-T compensator could be used in each control 
loop to eliminate the lightly damped modes, in favor of two poles on the 
negative real axis. The system gain could then be selected to correspond to 
the desired degree of damping in each loop. 

Following either compensator design procedure identified above, an 
effective closed loop design for the appendage control units would be 
obtained. For purposes of illustration, each control loop was configured 
with a phase lead compensator of the form 


G (s) « a 
c 


TS ♦ 1 
OTS + 1 


(54) 


with a » .1 auid T « 30. A CSMP run for the ROM indicates that with an 

initial angular off-set of .02 radians for each irass, the system was returned 
to their respective equilibrium points within 250 seconds as indicated in 
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Figure I2i Gsrehgorin BdkiuS for 
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Appendix 0> Faster responses could be obtained by adjusting the pole- 
zero positions in (54) or by an increase in compensator gains. 

In contrast to the Porcelli control units where 't^(t) and '^^^(t) are 
only functions of 6^(t) emd 02 (t), respectively/ the MNA control configura- 
tion is of the form 

Tg(t) - B^{t) + 03(t) (55) 

t^Q(t) » C 3 02(t) + 03(t) (56) 

The main control unit, however, is in correspondence with the Porcelli 
model, i.e., 

Tg(t) = C 3 0 j^(t) (57) 

where through are constants determined from K, L and F(s). 

In view of the above preliminary results using the MNA design philosophy, 
the following comments are in order; 

1. The MNA design method is an effective design tool for flexible space 
vehicles. 

2. The MNA method applied without consideration of the physical properties 
of the system, identifies the need for system decomposition. 

3 . FOi.' highly interactive systems of the POM form, the MNA method 
clearly demonstrates the need for a design which incorporates the inter- 
action into the control configuration rather than ignoring its presence. 

4* As a frequency domain technique for multivariable systems, the MNA 
method identifies the need for compensation in the proper control loops and 
also indicates the form of compensation required. 

5. The MNA design is readily adapted to higher order models and could 
be easily applied to more sophisticated satellite design configurations. 

This 18 not necessarily true of the Porcelli method of design. 
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Section 6s CONCLUSIONS 


It was the intent of this report to demonstrate the feasibility of 
the multivariable Nyquiat array design method to the attitude control design 
of flexible space vehicles* Using a three body model for the vehicle with 
coefficients specified in the Porcelli paper, three design approaches were 
considered: 

1* Standard single loop elimination design by classical procedures. 

2* MNA design for the full order model. 

3. MNA design for the reduced order model. 

In each case, the design configuration independently confirm the results 
of Porcelli. 

The first method conforms more directly to the standard compensator 
design approach identified in the traditional control texts. It is this 
design approach which is most likely to be used in the MNA compensator design. 
From the single loop elimination design in Section 3, it is readily apparent 
that this approach compares favorably with the Porcelli method. The root 
locus for the FOM transfer function clearly delineates the control 

problem and bandwidth considerations which must be taken into account in the 
design of the main control unit. An additional feature of this single loop 
approach lies in the recognition of subsystem interactions associated with 
the appendage control design as well as the utilization of bridged-T compen* 
sation. This form of compensation does not appear in Porcelli* s paper. 

The second design approach ignored any previous information pertaining 
to the pole-zero locations in the FOM. Although system dominance was not 
obtained in the MNA form, subsequent analysis revealed the problem to be 
associated with the rigid body dynamics. This condition presented itself 
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in the Nyquist diagram as a 180^ phase for all frequencies not associated 
with the resonant frequencies of the appendages. Analysis of the remaining 
MNA diagrams clearly indicates the need for lead compensation* 

To simplify the design for the appendage control units, the first mass 
was assumed to be at the equilibrium point. This assumption is easily 
justified from the time domain analysis via CSMP. An MNA design for the 
RON was made with results supportive of the FOM run. Lead compensation is 
clearly required for the control unit design. 

t 

► 
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APPENDIX A 


Transfer Matrix for Full Order Model 

Gij^(s) ■ (S^ 4 .OOOBS^ ♦ .00083^ + .16E-6S + .8E-7)/A 

■ (.00023^ ♦ .00023^ + .16E-63 ♦ .8E-7)/A 

Gj^^(s) - (.6E-7S^ + .16E-63 + .8E-7)/A 

Gjj^Cs) « (.4E-63^ + .00023^ + .00023^ 4 .16E-6S 4 .8E-7)/A 

G22(s) ■ (2S'* 4 .0013^ 4 .0013^ 4 .16E-63 4 .8E-7)/A 

G^jCs) - (.00083^ 4 .00083^ 4 .16E-6S 4 .8E-7)/A 

Gj^(s) - (.4E-63'’ 4 .32E-9S^ 4 .8032E-7S^ 4 .16E-6S 4 .8E-7)/A 

■ (-00083^ 4 00083^ 4 .16E-6S 4 .8E-7)/A 

G^^(s) » (43^ 4 .0023^ 4 .0023^ 4 .16E-63 4 .8E-7)/A 

A = 3® 4 .00093^ 4 .0009s'’ 4 . 28E-6S^ 4 .14E-63^ 
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.APPENDIX B 


Transfer Matrix for Reduced Order Model 

G^j(i) » (2S^ ♦ .0008S + .0008)/A 
“ (.0008S + .0008) /A 
G2,(s) = 

G22(s) = (43“^ + .00163 + .0016)/A 
A » 3^* + .00083^ + .00083^ + .16(10* /3 + .8(10”"^) 
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APPENDIX C 


Full order model under main loop control 
using compensator developed in Section 3. 
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APPENDIX D 


Full Order Model with Control Units 
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